search for: femp

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...