Hi
 I am a new user the R and I am very grateful for all your help but.......

I have a problem and I can't resolve yet. I am trying to get the contrasts for a binomial 3-way GLM (T= 4 temperature, t= 2 time and c= 2 substrate levels, plus treatment control) in total they are 17 treatments.

  

I have tried with the glht but this function only work for 1-way GLM, 

 

acacia<-cbind(g,N-g)

> acacia1<-glm(acacia~tratamiento,binomial)

 

> acacia2<-glht(acacia1, linfct = mcp(tratamiento = "Tukey"))

> summary(acacia2)

 

         Simultaneous Tests for General Linear Hypotheses

 

Multiple Comparisons of Means: Tukey Contrasts

 

 

Fit: glm(formula = acacia ~ tratamiento, family = binomial)

 

Linear Hypotheses:

            Estimate Std. Error z value p value    

B - A == 0    2.9073     1.0395   2.797  0.2417    

C - A == 0    4.0860     1.0227   3.995   <0.01 ** 

D - A == 0    4.8458     1.0207   4.748   <0.01 ***

E - A == 0    3.0445     1.0362   2.938  0.1727    

F - A == 0    3.3399     1.0306   3.241  0.0744 .........etc. 

 

If the GLM is with interactions don't work

 

Ų       acacia3<-glm(acacia~T*t*c,binomial)

Ų       > glht(acacia3, linfct = mcp(T = "Tukey"))

Ų       Error in modelparm.default(model, ...) : 

Ų         dimensions of coefficients and covariance matrix don't match

Ų       Warning message:

Ų       In mcp2matrix(model, linfct = linfct) :

Ų         covariate interactions found -- default contrast might be inappropriate

 

 

Other function could be simtes but this was deprecated from package multcomp

 

There are others options to get comparisons contrasts for model GLM with interactions? or this test don't exist ?

 

I hope help me 

 


	[[alternative HTML version deleted]]

