Displaying 1 result from an estimated 1 matches for "variedade2".
Did you mean:
variedade
2012 Nov 05
0
Diference in results from doBy::popMeans, multcomp::glht and contrast::contrast for a lme model
...nt 35 28 27 30 20 22 28 16 18 36 ...
# fixed effect model
m0 <- lm(y~bloco+variedade, da)
summary(m0)
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 35.7556 1.0004 35.742 < 2e-16 ***
## bloco10 -0.3778 1.2729 -0.297 0.770447
## ...
## variedade2 -9.2000 0.9262 -9.933 3.01e-08 ***
## variedade3 -8.0667 0.9262 -8.710 1.81e-07 ***
## ...
require(contrast)
contr <- contrast(m0, type="average", list(bloco=levels(da$bloco),
variedade="1"))
contr
## Contrast S.E. Lower Upper t df Pr(>|t|)...