R: formule utilisée dans predict()

Bonjour à tous,

Je suis coincé autour d'un problème de calcule d'un intervalle de prédiction après un modèle linéaire.

Je calcule ces intervalles pour (1) une régression simple du type Y ~ a0 + a1.X et (2) un modèle linéaire avec un paramètre supplémentaire (i.e. Y ~ a0 + a1.X1 + a2.X2).

J'ai trouvé la fonction predict() dans R qui en spécifiant interval=predition me calcule bien l'intervalle de prédiction, mais pas moyen de mettre la main sur le formule utilisée. D'une maniere générale, je n'ai pas trouvé non plus de réponse claire sur Google.

D'où : quelqu'un saurait-il (1) la formule utilisée par predict() et/ou (2) la formule communément acceptée pour le calcule d'un intervalle de prédiction pour ma régression linéaire et mon modèle linéaire a 2 variables indépendantes ?

Merci

Bonjour,

En tapant "predict.lm" dans R vous pouvez accéder au code de la fonction.
Je sais que ceci ne répond pas directement à votre question, mais peut-être qu'en se penchant un peu sur le code il est possible de trouver la formule...

Bonne continuation


Cordialement,

A.D.

Il s'agit ni plus ni moins de la formule théorique associée au modèle linéaire, regarde ce pdf par exemple :
http://amath.colorado.edu/courses/7400/2010Spr/lecture6.pdf

tu trouveras aussi la formule dans google books si tu tapes prediction interval.

Merci pour vos réponses qui me sont très utiles.

Cordialement,

PP.