Displaying 3 results from an estimated 3 matches for "lmes".
Did you mean:
lme
2004 May 13
3
GLMMs & LMEs: dispersion parameters, fixed variances, design matrices
Three related questions on LMEs and GLMMs in R:
(1) Is there a way to fix the dispersion parameter (at 1) in either glmmPQL (MASS) or GLMM (lme4)?
Note: lme does not let you fix any variances in advance (presumably because it wants to "profile out" an overall sigma^2 parameter) and glmmPQL repeatedly calls lme, so I...
2010 Aug 30
1
Help With Post-hoc Testing
...uot;,idvar="iobs",sep="",
dir="long")
> solution <- factor(solution,levels(1:17),)
> stack[1:5,]
iobs solution soi
32.1 32 1 7.00
70.1 70 1 5.95
78.1 78 1 3.00
104.1 104 1 3.75
117.1 117 1 5.00
> Lmes.mod <- lme(soi ~ solution + iobs, random = ~1 | iobs/solution, data =
stack)
> anova(Lmes.mod)
numDF denDF F-value p-value
(Intercept) 1 2383 2894.8342 <.0001
solution 1 2383 0.0003 0.9870
iobs 1 147 0.0126 0.9109
> summary(glht(Lme...
2010 Aug 05
3
eval-parse and lme in a loop
...=rnorm(9, mean=0.5, sd=0.7),
c3=rnorm(9, mean=0.9, sd=1)
);
## this works
means.lmeNoCovariate=lme(Mean_1 ~ Group, data=df, random = ~ 1 | Subject)
print(summary(means.lmeNoCovariate))
print(anova(means.lmeNoCovariate))
for (meanCol in paste("Mean_", 1:3, sep="")) {
## LMEs go here
cat("####################################################################################################\n")
cat("### AOV\n")
cat("####################################################################################################\n")
## this works...