Displaying 1 result from an estimated 1 matches for "ajmix".
Did you mean:
ajmi
2011 Jun 27
0
cld object did not plot
Dear R list.,
I am running a script to get a compact letter display.
library(lme4)
library(multcomp)
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...