Discussion :

[lilly74]

Bonjour,
j'ai des données (Y en fonction de age) que je trace avec une échelle logarithmique sur l'axe des ordonnées:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

plot(age, Y, log="y")

J'obtiens alors quelque chose qui ressemble à une droite.
J'ai besoin de faire une régression linéaire là-dessus et de récupérer les coefficients de la droite.
La fonction lm ne semble pas indiquée dans ce cas.
J'ai trouvé un code pour le graphique:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

plot(age, Y, log="y", panel.first=abline(lm(log10(Y) ~ age), col="red", lwd = 2), pch="+")

mais je ne sais pas récupérer ces fameux coefficients
Quelqu'un aurait une idée svp?
Merci d'avance

[lilly74]

Re-bonjour,
ce n'était pas compliqué, mais je cherchais une fonction alors que 2 minutes de réfléxion sont parfois plus efficaces
Au cas où d'autres têtes en l'air auraient ce problème :
Avec

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

lm(log10(Y) ~ age)

on obtient l'Intercept et le Beta1.
D'où :
Y = 10^(Intercept + Beta1 * age)