vocabulaire: "zone de confiance"? d'une droite de

[Stéphane Laurent]

Bonjour,

Je me demande comment on appelle les 2 courbes en pointillés ci-dessous, ou bien la région à l'intérieur de ces deux courbes.
On a une droite de régression y=ax+b et pour tout x on fait l'intervalle de confiance de ax+b, c'est comme ça que j'obtiens ces 2 courbes.



Désolé si ma question ne concerne pas R...

[Renaud Lancelot]

En anglais: pointwise confidence interval ou pointwise confidence region

En français, pas certain qu'il y ait une traduction "officielle" (et je serais heureux de l'apprendre). On peut tenter "intervalle de confiance point par point", ou "région de confiance point par point".

Renaud

[Logez Maxime]

Bonjour,

Borcard, régression linéaire

Maxime

[Renaud Lancelot]

Apparemment, ce doc ne répond pas précisément à la question. Voir ce post sur R-Help: http://finzi.psych.upenn.edu/R/Rhelp02a ... 67370.html

Renaud

[Logez Maxime]

Re,

il me semblait que la légende de la première figure donnait une piste.

Maxime

[Matthieu Stigler]

Hey, c'est intéressant ça... comment tu as fait pour obtenir ce polt? Merci!

[Stéphane Laurent]

Merci pour les liens.

Hey, c'est intéressant ça... comment tu as fait pour obtenir ce polt? Merci!

En bayésien, c'est simple comme bonjour. Tu obtiens des simulations de (la loi a posteriori de) a et b, ce qui te donne des simulations de ax+b, et il suffit de prendre les quantiles. Et en "mode non-informatif", ce sont bien des intervalles de confiance au sens fréquentiste.

Pour obtenir les simulations des paramètres de régression on peut utiliser le merveilleux package bayesm, "accompagné" par ce merveilleux livre http://faculty.chicagogsb.edu/peter.rossi/research/bsm.html

[Logez Maxime]

Re,

Les formules de calcul de ces intervalles de confiance sont dans tous les bouquins sur les modèles linéaires.

Code : Tout sélectionner

lm1 <- lm(y~x)
matplot(sort(x),cbind(y,predict.lm(lm1,interval="confidence"))[order(x),],type=c("p","l","l","l"),pch=1,lty=1,col="black")

Maxime