similar to: writing complex outputs to table

Dear R-developeRs, Attached follows a patch against svn 34959 that adds the printing of p-values to the TukeyHSD.aov function in stats package. I also updated the corresponding documentation file and added a 'see also' reference to the simint function of the multcomp package. As it was already brought up in a previous thread [1] in R-help, one can obtain the adjusted

Dear R-Helpers, I was calling the TukeyHSD function and not getting confidence intervals or p-values. It turns out this was caused by missing data and the fact that I had previously turned on R Commander (Rcmdr). John Fox knew that Rcmdr sets na.action to na.exclude, which causes the problem. If you have this problem, you can either exit Rcmdr before calling TukeyHSD or you can set na.action to

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

Hello, When plotting the results of a TukeyHSD multiple comparisons procedure with an ANOVA (lm) object, an extra line appears in the confidence intervals that contain 0. For example (this is straight from the TukeyHSD helpfile): > summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks)) > TukeyHSD(fm1, "tension", ordered = TRUE) > plot(TukeyHSD(fm1,

Quick question, as I attempt to learn R. For post-hoc tests 1) Is there an easy function that will take, say the results of tukeyHSD and create a grouping table. e.g., if I have treatments 1, 2, and 3, with 1 and 2 being statistically the same and 3 being different from both Group Treatment A 1 A 2 B 3 2) I've been stumbling over the proper syntax for simple effects for a tukeyHSD

Hi All, I am wondering if their is a convenient way to export the results of the TukeyHSD function to Word or Excel. I have used capture.output(tukey.contrast, file="tukey.contrast.xls") and this works, but the data are not in a table form, and so it is sort of a pain to manipulate the output. I have be unable to find a more convenient way and any assistance would be greatly

Hey, I have this TukeyHSD output from which I would like to extract only the P-values (p adj, last number). The problem is that the test output is a character list. How can I "break" this sentence to separate the Pv? Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = Fe1$Fe ~ Fe1$genotype) $`Fe1$genotype` diff lwr upr

hi list, i have to ask you again, having tried and searched for several days... i want to do a TukeyHSD after an Anova, and want to get the adjusted p-values after the Tukey Correction. i found the p.adjust function, but it can only correct for "holm", "hochberg", bonferroni", but not "Tukey". Is it not possbile to get adjusted p-values after

Dear useRs, I'm working on multiple comparison design on two factor (2 3 levels) ANOVA. Each of the tests I have tried (Tukey, multcomp package) seem to do only with one factor at a time. fm1 <- aov(breaks ~ wool * tension, data = warpbreaks) tHSD <- TukeyHSD(fm1, "tension", ordered = FALSE) $tension diff lwr upr p adj M-L -10.000000 -19.35342

Hello, I have some data in a data.frame where for each of a number of subjects, I have scores for all of a number of symptoms. Subjects are subdivided in a number of groups, which have unequal sizes. I'd like to plot between-group differences in the scores on the various symptoms. Ideally, that would be in a form as would be produced by > bwplot( Score~Symptom ) but I'm not sure

Hello, I am really new to R and it's still a challenge to me. Currently I'm working on my Master's Thesis. My supervisor works with SAS and is not familiar with R at all. I want to run an Anova, a tukey-test and as a result I want to have the tukey-grouping ( something like A - AB - B) I came across the HSD.test in the agricolae-package, but... unfortunately I do not get an output

[Moderator: This was erronously sent to R-announce, and filtered fortunately ] Hi, This is my first time to use R. I'm wondering which function I can use to do multiple comparison. I have an lm object and want to do multiple comparison based on this object. Splus has a function multicomp() but I could not find a similar one in R. Thanks in advance. Really appreciate it. Julie

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 all!

I'm new using R im trying to do a tukey test, but when i see the results the p value results in NA im guessing its because i have missing values but im not sure how to fix it AnovaModel.2 <- aov(area ~ trat, data=apilados) > summary(AnovaModel.2) Df Sum Sq Mean Sq F value Pr(>F) trat 4 11847 2961.76 9.9905 1.500e-06 *** Residuals 76 22531 296.46

Dear R-listers, I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I

Hey Folks: New to R and learning as I go, but really, I know just enough to get myself into trouble. I've waded through everything up till now, and don't see anything in the search that is directly helpful for the POS that I've created. The GOAL: All I want in the world is a program that performs 1-way ANOVA's on every column in a data set (taking the first column as the

I suppose that multcomp in R and multicomp in S-Plus are related and it appears that it is possible to use multicomp with lme in S-Plus given the following correspondence on s-news sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary and multicomp.default() Is it possible to use multicomp with lme in R and if so what is the syntax from a simple readily available

Hello, Is there an option I can easily set to put labels next to the points in a simple plot? The closest I got to it was to use the option pch, but this replaces the point symbol with the character specified, and it is not what I want. Thanks -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello everyone, So I ran an anova with aov and then I want to run post-hoc comparisons but keep receiving this message : > no applicable method for "TukeyHSD" Here is my code: > d<-read.table("d.txt") > d > Obs subj Hand Gaze RT > 1 1 s1 1 1 401.4 > 2 2 s2 1 1 363.3...... > summary(ano <-