similar to: plot() with TukeyHSD

R Community: In case you haven't seen it, I have written a manuscript called "The R Guide" that is available in the Contributed Documents section on CRAN. It is about 50 pages in length and it is written for the beginning R user. I use the guide in my undergraduate probability and math stat sequence, but it is pretty accessible for anyone with a basic understanding of statistics.

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

Hello, Version 2.2 of "The R Guide" is available for download in the Contributed Documents section on CRAN. "The R Guide" is written for the beginning R user. I use the guide in my undergraduate probability and math stat sequence, but anyone with a basic understanding of statistics (who wants to learn R) should find it useful. This updated version includes sections on

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

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 community, I'm trying to understand this behavior of TukeyHSD. My goal is to obtain defensible, labelled multiple comparisons of an interaction term. Firstly, if I plot the TukeyHSD from the model that calculates its own interactions, then the y-axis labels appear to be reflected on their median when compared to the text output of the TukeyHSD statement. The labels are integers.

R Version 2.0.0 Under development (unstable) (2004-05-07) I have been wanting to use TukeyHSD for two and three way aov's, as they are especially simple for students to use correctly. I realize model simplification is usually a preferred methodology, but my disciplinary enertia and simplicity of TukeyHSD is also compelling. However, the lack of names on the comparisons for the higher order

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

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

Hello, A newly-revised version of my introductory manual "The R Guide" is now available in the Contributed Documents section in CRAN. Among other changes, the description of the function t.test() has been expanded as well as an introduction to methods for the creation of prediction and confidence intervals with regression. The manual is still under 60 pages in length. The current

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

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

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:

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

I'm interested in restricting the pairwise comparisons of interaction effects in a multi-way factorial ANOVA, because I find comparisons of interactions between all different variables different to interpret. For example (supposing a p<0.10 cutoff just to be able to use this example): > summary(fm1 <- aov(breaks ~ wool*tension, data = warpbreaks)) Df Sum Sq Mean Sq F

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 am wondering how to customize a pairwise comparisons plot of a factorial ANOVA, without doing a lot of manual manipulation of a TukeyHSD object. The customizations I'd like are: 1. The aov used log-transformed response data, but I'd like to plot the intervals on their original, untransformed scales 2. Plot all the main and interaction effects together, rather than in a separate

So I've come across a few cases where complex outputs from functions will not write to tables. The most recent case involves the TukeyHSD function in the stats package. If I save the TukeyHSD call and print it, that obviously goes fine, but when I try writing to a table, I get an error message that says "cannot coerce class \multicomp\" into dataframe. What does this mean, and

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