Hi



I have a question concerning the denominator of an Fvalue.



Here is a partial print of the results :



Linear mixed model fit by REML

Formula: Y ~ mod + (1 | Bloc) + (1 | Bloc:mod)

   Data: df

  AIC  BIC logLik deviance REMLdev

 3704 3750  -1840     3750    3680

Random effects:

 Groups   Name        Variance Std.Dev.

 Bloc:mod (Intercept)  593.72  24.366

 Bloc     (Intercept)    0.00   0.000

 Residual             3207.74  56.637

Number of obs: 342, groups: Bloc:mod, 18; Bloc, 2



anova(res_lmer)



Analysis of Variance Table

    Df Sum Sq Mean Sq F value

mod  8 215186   26898  8.3854





Fvalue = MS "mod" / var_residual = 26898 / 3207



As the interaction bloc:mod is random, I would prefer to test MS "mod" against MS "bloc:mod".

Am I right to prefer that ?

If yes how can I get the MS "bloc:mod" ?



Thanks



florent




	[[alternative HTML version deleted]]


