Dans l'exemple ci-dessous je voudrais emboîter le facteur "dil" dans l'interaction "Opérateur*Equipement" :
Code : Tout sélectionner
> lmer(log(response) ~ (1|Operator) + (1|Equipement) + (1|Operator:Equipement) + (1|dil/Operator:Equipement), data=dat1)
Linear mixed model fit by REML
Formula: log(response) ~ (1 | Operator) + (1 | Equipement) + (1 | Operator:Equipement) + (1 | dil/Operator:Equipement)
Data: dat1
AIC BIC logLik deviance REMLdev
-61.86 -50.78 37.93 -79.46 -75.86
Random effects:
Groups Name Variance Std.Dev.
Operator:Equipement:dil (Intercept) 0.00000000 0.000000
Operator:Equipement (Intercept) 0.00000000 0.000000
dil (Intercept) 0.00000000 0.000000
Equipement (Intercept) 0.00039392 0.019847
Operator (Intercept) 0.00818002 0.090443
Residual 0.00536712 0.073261
Number of obs: 36, groups: Operator:Equipement:dil, 12; Operator:Equipement, 4; dil, 3; Equipement, 2; Operator, 2
Fixed effects:
Estimate Std. Error t value
(Intercept) 3.5039 0.0666 52.61
En sortie j'obtiens le facteur "dil" dans les Random efects, pourquoi ? Avec SAS je ne l'obtiens pas :
Code : Tout sélectionner
Operator 0.008180
Equipement 0.000394
Operator*Equipement 0
dil(Operato*Equipem) 0
Residual 0.005367