Return to Answer

added 80 characters in body

edited May 15, 2024 at 14:42

85.6k
32
216
694

Do you mean you want to do all pairwise comparisons for the three factors?

mod1<-lme(Variable~Sediment*Hydrology*Depth, data=mydata, random=~1|Site/Hydrology/Depth) mydata$SHD<-interaction(mydata$Sediment,mydata$Hydrology,mydata$Depth) mod2<-lme(Variable~-1+SHD, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2,linfct=mcp(SHD="Tukey")))

mod1 <- lme(Variable ~ Sediment*Hydrology*Depth, data=mydata, random= ~ 1|Site/Hydrology/Depth) mydata$SHD <- interaction(mydata$Sediment, mydata$Hydrology, mydata$Depth) mod2 <- lme(Variable ~ -1 + SHD, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2, linfct=mcp(SHD="Tukey")))

Do you mean you want to do all pairwise comparisons for the three factors?

mod1<-lme(Variable~Sediment*Hydrology*Depth, data=mydata, random=~1|Site/Hydrology/Depth) mydata$SHD<-interaction(mydata$Sediment,mydata$Hydrology,mydata$Depth) mod2<-lme(Variable~-1+SHD, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2,linfct=mcp(SHD="Tukey")))

Do you mean you want to do all pairwise comparisons for the three factors?

mod1 <- lme(Variable ~ Sediment*Hydrology*Depth, data=mydata, random= ~ 1|Site/Hydrology/Depth) mydata$SHD <- interaction(mydata$Sediment, mydata$Hydrology, mydata$Depth) mod2 <- lme(Variable ~ -1 + SHD, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2, linfct=mcp(SHD="Tukey")))

Source Link

answered Feb 21, 2012 at 17:10

lockedoff

Do you mean you want to do all pairwise comparisons for the three factors?

mod1<-lme(Variable~Sediment*Hydrology*Depth, data=mydata, random=~1|Site/Hydrology/Depth) mydata$SHD<-interaction(mydata$Sediment,mydata$Hydrology,mydata$Depth) mod2<-lme(Variable~-1+SHD, data=mydata, random=~1|Site/Hydrology/Depth) summary(glht(mod2,linfct=mcp(SHD="Tukey")))