Displaying 2 results from an estimated 2 matches for "mixmodel".
2011 Jun 27
0
cld object did not plot
...omp)
library(gplots)
####### Mixed Effects Model #########
data <- read.table("AJmix.txt",header=TRUE, sep="\t")
attach(data)
y<-cbind(positive,negative)
treatment<-factor(treatment)
mouse<-factor(mouse)
data$obs<-1:nrow(data)
names(data)
detach(data)
attach(data)
Mixmodel<-lmer(y~treatment+(1|mouse:treatment)+(1|obs) ,family=binomial)
summary(Mixmodel)
model.mc<-glht(Mixmodel,linfct=mcp(treatment="Tukey"))
summary(model.mc, test=adjusted(type="Westfall"))
model.ci <- confint(model.mc, level= 0.95)
model.ci
par(mar=c(5,10,3,2)+.1)
plot(mo...
2013 Feb 27
1
metafor - interpretion of QM in mixed-effects model with factor moderator
...I've read previous posts about QM in
metafor, but I'm still a bit confused. I've dummy-coded the factors:
MData$k.animal<-ifelse(MData$Type=="Animal",1,0)MData$k.plant<-ifelse(MData$Type=="Plant",
1,0)
I've then run the model with the following script:
MixModel<-rma(z, v.z., mods=cbind(k.animal, k.plant), intercept=FALSE,
data=MData, subset=(Response=="Tolerance"))MixModel
I receive the output:
Test of Moderators (coefficient(s) 1,2):
QM(df = 2) = 8.5758, p-val = 0.0137
Model Results:
estimate se zval pval ci.lb...