Displaying 2 results from an estimated 2 matches for "femp".
Did you mean:
fcmp
2005 Oct 27
2
F tests for random effect models
...ary(aov(Rendement ~ Pollinisateur * Lignee, data = mca2)) gives wrong tests of random
effects. But, summary(aov1 <- aov(Rendement ~ Error(Pollinisateur * Lignee), data = mca2)) gives no
test at all, and I have to do it like this :
tab1 <- matrix(unlist(summary(aov1)), nc=5, byrow=T)[,1:3]
Femp <- c(tab1[1:3, 3]/tab1[c(3,3,4), 3])
names(Femp) <- c("Pollinisateur", "Lignee", "Interaction")
1 - pf(Femp, tab1[1:3,1], tab1[c(3,3,4),1])
With "lme4" package (I did'nt succeed in writing a working formula with lme from "nlme" package)...
2005 Oct 28
2
Random effect models
...Lignee
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 4 18.0294 4.5074
Error: Pollinisateur:Lignee
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 36 5.1726 0.1437
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 50 3.7950 0.0759
# F tests :
> Femp <- c(tab1[1:3, 3]/tab1[c(3,3,4), 3])
> names(Femp) <- c("Pollinisateur", "Lignee", "Interaction")
> Femp
Pollinisateur Lignee Interaction
9.258709 31.370027 1.893061
> 1 - pf(Femp, tab1[1:3,1], tab1[c(3,3,4),1])
Pollinisateur...