Bonjour à tous,
J'ai un tableau de données avec plusieurs variables : la biomasse, la longévité (variable catégorielle), l'impact (variables continus) et l'habitat (facteurs) :
Station Biomasse Longévité Impact habitat
1 100 4 0.005 A5.1
1 52 2 0.005 A5.1
1 15 3 0.005 A5.1
2 20 2 0.6 A5.2
2 5 4 0.6 A5.2
3 54 3 0.07 A5.1
3 60 1 0.07 A5.1
J'ai ensuite fait un glm (qui me permet d'obtenir les coefficients nécessaire pour la suite de mes calculs) :
a<- glm(log(biomasse)~log(Longévité)+ habitat + impact + impact*habitat + log(longévité)*impact + log(longévité)*habitat)
J'obtiens donc des coefficients que j'utilise dans une autre équation où je m'intéresse à chaque habitat séparement :
B0 + B1*log(10) + B2_51+B3*log(0.01) + B4_51*log(0.01) + B5_51 *log(10) + B6*log(10)*log(0.01)
donc Bo = intercept
B1 = coefficient pour la variable longévité
B2_51 = coefficient pour l'habitat A5.1
....
Comme je dois faire ce calcul par habitat, j'aimerais faire une boucle dans laquelle je récupère les coefficients et je change juste les coefficients qui change pour chaque habitat (donc B2, B3 et B5). Es ce possible?
Merci d'avance pour votre aide
Cyrielle