Afficher les p-value de lm() sur un graphe

[Baptiste drouin]

Bonjour, je réalise plusieurs corrélation linéaire avec lm() sur différents graphiques. Je souhaiterais mettre automatiquement sur ces dernier les valeurs des p-value correspondante au lm(). Je parviens à le faire pour le R^2 comme ci dessous...

Code : Tout sélectionner

ptm=read.table("p_tjuin.txt")
y<-tp[1:175,2]
x<-RIm[94:268]
lin<-lm(x~y)
plot(y,x,xlim=c(0,1),main="Ratio juin (precipitation / température) /croissance des cernes")
summary(lin)
abline(lin)
r2 <- summary(lin)$r.squared
mtext(bquote(  R^2 == .(r2)))


J'ai cherché comment la commande lm appelle les p-value pour pouvoir sortir les données de summary(lm) mais je n'est pas réussi .
Si quelqu'un à la solution ou une autre astuce je suis preneur.
Merci

[Maxime Hervé]

Bonjour,

les p-values sont dans le tableau 'coefficients', accessible de cette façon :

Code : Tout sélectionner

summary(lin)$coefficients

Ce qui est retourné est une matrice, donc vous pouvez vous ballader dedans comme dans n'importe quelle matrice. La colonne des p-values est ainsi :

Code : Tout sélectionner

summary(lm.D9)$coefficients[,"Pr(>|t|)"]

Maxime

[Baptiste drouin]

merci beaucoup ca marche