similar to: Tukey in R, extracting values

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

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

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

Hello, I am trying to determine the most appropriate way to run post-hoc comparisons on my lme model. I had originally planned to use Tukey HSD method as I am interested in all possible comparisons between my treatment levels. TukeyHSD, however, does not work with lme. The only other code that I was able to find, and which also seems to be widely used, is glht specified with Tukey:

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

Dear all, I was trying to understand how "multcomp" package works by running the examples given in the documentation. However I still don't understand when it comes to multiple comparison set by user (please refer to "Ksub" in the code). Therefore I run 2 other cases along with the original example (case 1), with the expectation I'll get the point from the output. The

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

Dear all, glht (from the multcomp package) needs a term and a model component in it's fitted model. In fitted models from e.g. repeated measurements ANOVAs I do not find neither model nor term. Is it possible to build together a model and term component myself, so that glht will work for repeated measurements ANOVAs? If so, how would I do that? Best regards, Michael Zehetleitner

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 All, Another nls related problem (for background, I'm migrating a complicated modelling package from S-plus to R). Below I've reduced this to the minimum necessary to demonstrate my problem (I think); the real situation is more complicated. Two similar selfStart functions, ssA and ssB. The 'initial' function for ssB modifies its arguments a little and then calls getInital

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

Hi, I've got an HP ZBook with C6.5 and KDE. Yesterday I was able to play music, today I can't get any sound out of my pc. lspci -v 00:1b.0 Audio device: Intel Corporation 8 Series/C220 Series Chipset High Definition Audio Controller (rev 04) Subsystem: Hewlett-Packard Company Device 1909 Flags: bus master, fast devsel, latency 0, IRQ 34 Memory at cca30000 (64-bit, non-prefetchable)

Hi All, I have performed a 3-way ANOVA analysis for my experimental data using aov function. My simple R funtion for this is: 3aof <- function(x){ m <- data.frame(R,S,T, x); anova(aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) ) } Now, I am getting P values for all the main and interactions effects. If I want to perform postdoc test on one of my main effects, say T, what method I should use (i have

Hi, I am trying to setup a Driver Domain for a NIC. I use ubuntu Jaunty Dom0 and the same kernel for the Driver Domain as well. I have successfully pass-through-ed the PCI device to the DomU ( Driver Domain). But i am not sure as to what i will have to do next to get the IDD functional ? -- -- With Regards , Amod Jaltade Good, Fast, Cheap: Pick any two (you can''t have all

Hello, I am performing a BACI analysis with ANOVA using the following glm: fit1<-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1) The summary(aov(fit1)) shows significance in the interaction; however, now I would like to determine what combinations of BA and Otter are significantly different (each factor has two levels). ID and PRIMARY substrates are categorical and included in

Dear Colleagues, I fit this model: mod1 <- lmer(x~category*comp+(1|id),data=impchiefsrm) where category has 4 levels and comp has 8 levels. These work: glht(mod1, linfct=mcp(category="Tukey") glht(mod1, linfct=mcp(comp="Tukey") What I'd like is (conceptually): glht(mod1, linfct=mcp(category:comp="Tukey") but it gives a syntax error. Any help is

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 +

Hi everyone I'm using Tukey HSD as post-hoc test following a lme analysis. I'm measuring hemicelluloses in different species treated with three different CO2 concentrations (l=low, m=medium, h=high). The whole experiment is a split-plot design and the Tukey-function from the package multcomp is suitable for lme-analysis with random factors. The analysis works fine but I get a non

Thanks in advance for your help. I am running a repeated measures ANOVA in r. The same group undergoes to four different treatment conditions. So, all individuals are treated with treatments A, B, C and D in four different occasions. Once I get a significant ANOVA, I first run a paired samples t-test using the code: t.test(X1,X2,paired=TRUE) #being x1 the punctuation after treatment 1 and x2 the