1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28
|
head(dat)
traitement origine tot_inflo L_mm En_fruits code_plante bloc
1 O BQ 13 25 0 BQ_04_01 2
2 O BQ 13 25 0 BQ_04_01 2
3 O BQ 13 25 0 BQ_04_01 2
4 O BQ 13 25 0 BQ_04_01 2
5 O BQ 13 25 0 BQ_04_01 2
6 O BQ 13 37.700000000000003 0 BQ_04_01 2
num_hampe_homogene
1 H1
2 H1
3 H1
4 H1
5 H1
6 H2
glmer_fruit_set =
glmer(cbind(En_fruits, 1-En_fruits)
~ traitement*origine*L_mm + traitement*origine*tot_inflo
+ (1|code_plante)
+ (1|num_hampe_homogene)
+ (1|bloc),
data = dat,
family = "binomial",
control=glmerControl(optimizer="bobyqa",
optCtrl=list(maxfun=100000)))
glmer_fruit_set |
Partager