Bonjour,
Je voudrais reproduire une ellipse pour un nuage de points mais de grosseur variable suivant l'intervalle de confiance choisi.
Pour l'instant j'ai trouvé quelque chose du genre :
x<-rnorm(100)
y<-x+10-rnorm(100)
c<-cor.test(x,y)
c<-c[[4]]
plot(x,y)
reg.line(lm(y~x))
xy<-cbind(x,y)
exy <- ellipsoidhull(xy)
lines(predict(exy))
qui me donne une ellipse incluant tous les points. Or en général, on définit plutôt un paquet un peu plus petit. Je n'arrive pas trop à comprendre ce qui est utilisé dans le package ellipse...