similar to: multcomp and lme4

Hi R users, I sent a message yesterday about NA in model estimates ( http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html). If I use aov() instead of lm() I get no NA in model estimates and I use gmodels::estimable() without problems. Ok! Now I'm performing a lot of contrasts and I need correcting for multiplicity. So, I can use multcomp::glht() for this.

Hello everyone, I've conducted a Type III repeated-measures ANOVA using Anova() from the car package, based on the suggestions at http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/(option 3) and http://languagescience.umd.edu/wiki/EEG#ERP_ANOVA_in_R. My ANOVA has two factors: Condition (3 levels) and Region (6 levels) and their interaction. Below is code to run the Anova

Hi R Forum I am a newbie to R and I have been amazed by what I can get my team to accomplish just by implementing Scripting routines of R in all my team's areas of interest.. Recently i have been trying to adopt R scripting routine for some analysis with longitudanal data.. I am presenting my R script below that I have tried to make to automate data analysis for longitudanal data by employing

Hi All, I am trying use 'multcomp' for multiple comparisons after my ANOVA analysis. I have used the following code to do ANOVA: dat <- matrix(rnorm(45), nrow=5, ncol=9) f <- gl(3,3,9, label=c("C", "Tl", "T2")) aof <- function(x) { m <- data.frame(f, x); aov(x ~ f, m) } amod <- apply(dat,1,aof) Now, how can I use

Using a factor with 'number' levels the straightforward symbolic formulation of a contrast in 'glht' of the 'multcomp' package fails. How can this problem be resolved without having to redefine the factor levels? Example: #A is a factor with 'number' levels #B similar factor with 'letter' levels dat<-data.frame(y=1:4,A=factor(c(1,1,2,2)),

Dear R-users, I'm currently using the multcomp package to produce plots of means with 95% confidence intervals i.e. mult<-glht(lm(response~treatment, data=statdata), linfct=mcp(treatment="Means")) plot(confint(mult,calpha = sig)) Unfortunately the y-axis on the plot appears to be fixed and hence if the labels on the y-axis (treatment levels) are too long, then they are not

I am running a linear model with achiev as the outcome and major as my iv (5 levels). The lm statement runs fine, but for the glht command I get the following error. I noted that someone else asked the same question a while back but received no reply. I am hoping someone might know what is happening. anovaf2<-lm(achiev ~ major, data=data_mcp) > pairwise<- glht(anovaf2,linfct =

I am trying to do Tukey HSD comparisons on a repeated measures expt. I found the following example on r-help and quoted approvingly elsewhere. It is broken. Can anyone please tell me how to get it to work? I am using R 2.4.1. > require(MASS) ## for oats data set > require(nlme) ## for lme() > require(multcomp) ## for multiple comparison stuff > Aov.mod <- aov(Y ~ N + V +

I used the multcomp package sometime back for doing multiple comparisons. I see that it has been updated and the methods like simint are no longer supported. When I run the program it prompts to me to use glht. How do I get the lower and upper conf int and the pValues using glht? Does anyone have an example? Thanks ../Murli [[alternative HTML version deleted]]

Hi, I am experimenting with using glht() from multcomp package together with coxph(), and glad to find that glht() can work on coph object, for example: > (fit<-coxph(Surv(stop, status>0)~treatment,bladder1)) coxph(formula = Surv(stop, status > 0) ~ treatment, data = bladder1) coef exp(coef) se(coef) z p treatmentpyridoxine -0.063 0.939 0.161

Dear R members, I tried to fit an lme model and to use the glht function of multcomp. However, the glht function gives me some errors when using weights=varPower(). The glht error makes sense as glht needs factor levels and the model works fine without weights=. Does anyone know a solution so I do not have to change the lme model? Thanks Sibylle --> works fine

Hi, I have a fairly simple repeated measures-type data set I've been attempting to analyze using the lme function in the nlme package. Repeated searches here and other places lead me to believe I have specified my model correctly. However, I am having trouble with post-hoc tests. From what I gather, other people are successfully using the glht function from the multcomp package to

Hi, I would like to derive p-values for pair-wise comparison (Tukey's) of effects when the response is a count. I am trying a test case where y ~ Po( lambda(x) ). x has three levels : A, B and C with lambda(x) = 10, 20 and 20 respectively. Hence, p-values for the contrast C - B should distribute uniformally. I have implemented this test case as below but do not get uniform

I am trying to use multcomp to do a Tukey posthoc on growth increments among genetic crosstypes. #Fixed effect model m1 <- lm(inc ~ 0 + Age+ Crosstype + Sex, data = Data.age) summary(m1) RESULTS of the model: summary(m1) Call: lm(formula = inc ~ 0 + Age + Crosstype + Sex, data = Data.age) Residuals: Min 1Q Median 3Q Max -0.87180 -0.34002 -0.02702 0.27710 2.17820

Hi I have some experimental data where I have counts of the number of insects collected to different trap types rotated through 5 different location (variable -location), 4 different chemical attractants [A, B, C, D] were applied to the traps (variable - semio) and all were trialled at two different CO2 release rates [1, 2] (variable CO2) I also have a selection of continuous variables

Hi R users: I have de following data frame (called "Sx") Descripcion Nitratos Cont85g 72.40 Cont85g 100.50 Cont85g 138.30 Cont80g 178.33 Cont80g 79.01 Cont80g 74.16 Cont75g 23.70 Cont75g 15.80 Cont75g 16.20 Patron80g

i ran this model > model2<-glm(rojos~ageandsex+sector+season+sector:season,quasipoisson) > glht(model2,linfct=mcp(ageandsex="Tukey")) General Linear Hypotheses Multiple Comparisons of Means: Tukey Contrasts Linear Hypotheses: Estimate M - H == 0 0.2898 SUB - H == 0 -0.2261 SUB - M == 0 -0.5159 I tried to do the same changing factor season

Hi, I am trying to test some contrasts, using glht() in multcomp package on fixed effects in a linear mixed model fitted with lme() in nlme package. The command I used is: ## a simple randomized block design, ## type is fixed effect ## batch is random effect ## model with interaction dat.lme<-lme(info.index~type, random=~1|batch/type, data=dat) glht(dat.lme, linfct = mcp(type

Hi, I have been trying glht() from multcomp package and contrast() from contrast package to test a contrast that I am interested in. With the following simulated dataset (fixed effect "type" with 3 levels (b, m, t), and random effect "batch" of 4 levels, a randomized block design with interaction), sometimes both glht() and contrast() worked and gave nearly the same p values;

Hi, I am running a simple one-way ANOVA with an independent factot variable "treat" (3 levels: a, b and c) and a response variable "y". I want to test a linear relationship of the response among the 3 levels of the variable "treat" (ordered a->b->c). I used glht() from multcomp package. Later I found out I need to exclude the situation where the response at the