Bonjour,
Je débute dans R, j'ai encore du mal à réaliser des boucles. Ici j'aimerai faire un glm en boucle pour ma variable de 3 échantillons, puis ressortie le glht sous format png. Le code, hors boucle, fonctionne parfaitement, avec boucle, cela ne fonctionne pas, je n'ai aucune sortie png :
for (echantillon in 1:3) {
data_bino2<-data_bino[data_bino$N_SERIE==serie]
data_bino2$Con<-as.factor(data_bino2$Con)
data_bino2$Antago<-as.numeric(data_bino2$Antago)
M1<-glm("Antago~Con",family=binomial(link=logit), data=data_bino2)
ANO<-anova(M1, test="Chi")
A<-summary(glht(model=M1, linfct = mcp(Con = "Tukey")))
B<-table_glht(A)
C<-data.frame(B)
C$Significativité<-NA
C[1,5]<-ETOIL(C[1,4]) # Fonction pour ajouter les étoiles de significativité
C[2,5]<-ETOIL(C[2,4])
C[3,5]<-ETOIL(C[3,4])
C[4,5]<-ETOIL(C[4,4])
C[5,5]<-ETOIL(C[5,4])
C[6,5]<-ETOIL(C[6,4])
tt<- ttheme_minimal(
core=list(bg_params = list(col="grey"),
fg_params=list(fontface=10), fontface=4),
colhead=list(fg_params=list(col="black", fontface=4)),
rowhead=list(fg_params=list(col="black", fontface=10L)))
png(paste(filename = "STAT_5_T3_CON_GLHT_", HE,".png"), width = 600, height = 200)
grid.table(C, theme=tt)
dev.off()
}
D’avance merci !