je continue mes débuts sur R.
J'essaye de dessiner une régression avec intervalles de confiance avec ggplot2.
Voici mon code. Il dessine une droite de régression à l'extérieur de l'intervalle de confiance !!

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
#charge la librairie (à ne faire qu'une fois)
install.packages("ggplot2", dependencies = TRUE)
 
#monte la librairie principale et ses dépendants
# doc sur file:///C:/Documents%20and%20Settings/planchoo/Bureau/ggplot2/html/00Index.html
library(ggplot2)
 
#ouvre le jeu de données diamonds inclu dans le package ggplot2
data("diamonds", package = "ggplot2")
 
#crée un sous-ensemble aléatoire de 100 lignes pour certains graphiques
set.seed(1410)
dsmall <- diamonds[sample(nrow(diamonds), 100),]
 
#effectue la régression
model <- lm(price ~carat, data=dsmall) 
 
# crée les données X à plotter pour le fit
carat = seq(from = min(dsmall$carat), to = max(dsmall$carat), length.out = 20)
nd = data.frame(carat)
 
# Prédictions du modèle sur les nouvelles données
err <- stats::predict(model, newdata=nd, se = TRUE)
 
# affecter la prédiction +/- 2 stdev à de nouvelles variables
nd$ucl <- err$fit + 1.96 * err$se.fit 
nd$lcl <- err$fit - 1.96 * err$se.fit 
 
#plot les points + les prédictions avec intervalle de confiance
qplot(carat, price, data=dsmall, alpha=0.5, xlim = c(0,3)) +   geom_smooth(aes(x=carat, ymin = lcl, ymax = ucl), data=nd, stat="identity")
le problème est dans la dernière ligne :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
#plot les points + les prédictions avec intervalle de confiance
qplot(carat, price, data=dsmall, alpha=0.5, xlim = c(0,3)) +   geom_smooth(aes(x=carat, ymin = lcl, ymax = ucl), data=nd, stat="identity") 
[