Barres d'erreur

[Guillaume Jubelin]

Bonjour,
j'aurais voulu savoir comment rajouter des barres d'errreurs sur mes plots
Merci d'avances

[Renaud Lancelot]

Question un peu vague. Voir la fonction segments (ou lsegments pour des graphiques trellis) qui permet de tracer des segments de droite. Exemple avec une des fonction trellis:

Code : Tout sélectionner

foo <- data.frame(
  y = rep(c("A", "B"), each = 10),
  x = rep(c(2, 5), each = 10) + rnorm(20))

library(lattice)
dotplot(y ~ x, data = foo,
  panel = function(x, y){
# trace les points
    panel.dotplot(x, y, pch = 1, cex = 1.5)
# moyenne par groupe
    m <- tapply(x, y, mean)
# écart-type de m
    se <- tapply(x, y, function(x) sqrt(var(x) / length(x)))                   
# y est un facteur: on doit le transf en nbre
    Y <- as.numeric(unique(y))
# trace la moyenne
    lpoints(m, Y, pch = 8, cex  = 2, col = "red")
# trace IC 95%
    lsegments(m - 1.96 * se, Y, m + 1.96 * se, Y, lwd = 2, col = "red")
    })


[Guillaume Jubelin]

Merci beaucoup,
je vais tester cela.
En fait, en un point donné d'un plot je veux pouvoir tracer une barre veticale représentant la différence entre la valeur estimée (celle du point) et la valeur mesurée en réalité.

Encore merci

[Renaud Lancelot]

OK, alors voir segments ou lsegments, de manière très semblable à ce que j'ai proposé dans l'exemple.

Renaud