similar to: Is using glht with "Tukey" for lme post-hoc comparisons an appropriate substitute to TukeyHSD?

Hi everyone, I am fairly new to R, and I am aware that others have had this problem before, but I have failed to solve the problem from previous replies I found in the archives. As this is such a standard procedure in psychological science, there must be an elegant solution to this...I think. I would much appreciate a solution that even I could understand... ;-) Now, I want to calculate a

Hi, I want to export the result of glht in R into a LaTeX table, such as that result: Linear Hypotheses: Estimate Std. Error z value Pr(>|z|) Group1 - Group2 == 0 -0.14007 0.01589 -8.813 <0.001 "***" Group1 - Group3 == 0 -0.09396 0.01575 -5.965 <0.001 *** --- Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05

I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function

Hi, I'm trying to make sense of the options for multiple comparisons options in R. I've found the following options: pairwise.t.test, which provides standard t-tests, with options for choosing an appropriate correction for multiple comparisons TukeyHSD, which provides the usual Tukey test glht(package multcomp), which provides a variety of options >From the help list, it appears

hey, how can i extract the values from the CI's when i use following code for a tukey test? the output shows three CI's and i know it should work with 'names but i don't really know how... thanks library(multcomp) data1$soort<-as.factor(data1$soort) amod<-aov(waarde~soort,data=data1) g<-glht(amod, linfct=mcp(soort = "Tukey")) confint(g) -- View this message in

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

Hi, I am seeking help with an error when running the example from R Documentation for TukeyHSD. The error occurs with any example I run, from any text book or website. thank you... > plot(TukeyHSD(fm1, "tension")). Error in plot(confint(as.glht(x)), ylim = c(0.5, n.contrasts + 0.5), ...) : error in evaluating the argument 'x' in selecting a method for function

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]]

I have a question about TukeyHSD and the glht function because I'm getting different answers when a covariate is included in model for ANCOVA.? I'm using the cabbages dataset in the 'MASS' package for repeatability.? If I include HeadWt as a covariate, then I get different answers when performing multiple comparisons using TukeyHSD and the glht function. The difference appears

I am not sure whether there is a bug. When I tested the example given for "glht" in the help, I entered the following error: Running commands: amod <- aov(minutes ~ blanket, data = recovery) rht <- glht(amod, linfct = mcp(blanket = "Dunnett"), alternative = "less") Errors are: Error in try(coef.(model)) : could not find function

Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, "Four Quartets: Burnt Norton" -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

Could you please help me on the following ANCOVA issue? This is a part of my dataset: sampling dist h 1 wi 200 0.8687212 2 wi 200 0.8812909 3 wi 200 0.8267464 4 wi 0 0.8554508 5 wi 0 0.9506721 6 wi 0 0.8112781 7 wi 400 0.8687212 8 wi 400 0.8414646 9 wi 400 0.7601675 10 wi 900 0.6577048 11 wi 900

Dear list members, I'm running some tests looking at differences between means for various levels of a factor, using Tukey's HSD method. I would like to plot the data as boxplots or dotplots, with horizontal significance lines indicating which groups are statistically significantly different, according to Tukey HSD. Here's a nice image showing an example of such a graphical

Dear *, I would like to obtain for each factor of my anova model the "response variable vs factor" plot with means and 95% Tukey HSD intervals. I would appreciate any information on how to do that. Cheers -------------------------------------------------------------------- Max MANFRIN Tel.: +32 (0)2 650 3168 IRIDIA - CoDE, CP 194/6

Hello, I'm a totally newbie to R and I'm taking a class using S+. In the class we use the multcomp command which takes a aov object and calculates confidence intervals for all pairwise differences by the Fisher least significant differences method. How can I do this in R. Thank you for taking the time with such a basic question. I've been looking on the net for a few days and I

Hi R helpers! I have used a glht as a post-hoc test on an lmer with: -2 treatments (A & B) -1 categorical explanatory variable (song type) -1 response variable (latency to respond) I wanted to make comparisons between the categorical variables depending on treatment. At the moment the glht simply returns comparisons of each of the (3) categorical explanatory variables with each other

Hi, I'm wondering if i can use the function "TukeyHSD" to perform the all pairwise comparisons of a "aov()" model with one factor (e.g., GROUP) and one continuous covariate (e.g., AGE). I did for example: library(multcomp) data('litter', package = 'multcomp') litter.aov <- aov(weight ~ gesttime + dose, data = litter) TukeyHSD(litter.aov, which =

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

Greetings, In the following one-way ANOVA I am attempting to calculate the means of each treatment along with their 95% Tukey confidence intervals for the data shown below using a routine from the HH package. library(HH) options(digits=10) # load data treat voltage 1 130 1 74 1 155 1 180 2 150 2 159 2 188 2 126 3 138 3 168 3 110 3 160 4 34

Dear R people, I have a couple of questions about post-doc analyses for 2 by 2 within subjects ANOVA. I conducted a psycholinguistic study that combined a 2 by 2 design and a latin square design. Specifically, I had 32 items each of which generated 4 conditions. Participants saw each of the 32 items only once: 8 in Condition A, 8 in B, 8 in C, and 8 in D. The table below serves as an example.