Thanks for the suggested code! It’s a very nice way of displaying most
aspects of the data and the HSD tests/CIs. The graphical display is
probably to big for inclusion in journal articles, but works well for
displaying the results when working with the data.
Regards,
Karl Ove Hufthammer
ty. de
Karl,
I will look at your function later and comment on it.
RIght now, I am replying to your comment "The levels should preferably be
ordered by the group means, "
One of the features of the tabular printout from the mmc function, and of
both
the MMC plot and the plotMatchMMC plot, is that t
Den 2013-01-14 19:58 skreiv Richard M. Heiberger:
Please look at the MMC (Mean-mean Multiple Comparisons) plot in the
HH package.
It displays both the means and the differences.
install.packages("HH") ## if you don't already have it.
library(HH)
?MMC
I have now coded a very quick-and-dirty so
Karl,
Thank you for this interesting challenge. It turns out to be reasonably
straightforward.
This code does it for the first example in ?MMC.
## This example is based on the example in ?HH:::MMC
library(HH)
## pairwise comparisons
data(catalystm)
## A. individual plots
## base graphics bo
When the group means are close together, the labels will overlap.
This overlap is actually informative, indicating that the groups are close.
For this common case, we provide the tiebreaker function
matchMMC
which is also documented and illustrated on the ?MMC page.
I have never tried to put the
må. den 14. 01. 2013 klokka 13.58 (-0500) skreiv Richard M. Heiberger:
> Please look at the MMC (Mean-mean Multiple Comparisons) plot in the HH
> package. It displays both the means and the differences.
>
> install.packages("HH") ## if you don't already have it.
> library(HH)
> ?MMC
Thanks for
Please look at the MMC (Mean-mean Multiple Comparisons) plot in the HH
package.
It displays both the means and the differences.
install.packages("HH") ## if you don't already have it.
library(HH)
?MMC
Rich
On Mon, Jan 14, 2013 at 12:36 PM, Karl Ove Hufthammer wrote:
> Dear list members,
>
> I'
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
Kingston ON Canada
> -Original Message-
> From: ulrikebraeck...@hotmail.com
> Sent: Mon, 6 Aug 2012 07:19:09 -0700 (PDT)
> To: r-help@r-project.org
> Subject: [R] Tukey HSD not fully displayed in R console
>
> Dear all,
>
> I would like to test the differences in de
Dear all,
I would like to test the differences in dependent variable X depending on 2
grouping variables of each 10 levels.
I do this with a 2-way ANOVA, followed by a Tukey HSD test (TukeyHSD(x)).
However, since a lot of combinations are possible with 2 grouping variables,
each of 10 levels, the
Also you want " for(i in 2:ncol(mydf)) { ..."
Your current setup of 1:ncol(mydf)-1 picks columns 0,1,2,...,ncol-1
(missing the last), which is not what you want...
Setting i<-2 in the first line was overridden by the call to loop.
Cheers
On Thu, Feb 9, 2012 at 2:33 PM, peter dalgaard wrote:
>
>
On Feb 9, 2012, at 20:28 , jlpugh wrote:
> 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.
>
>
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 perform
Try TukeyHSD(anova). TukeyHSD is implemented in the base package, i.e., ready
to use with the base installation of R. Also, the TukeyHSD for the
interaction term should have a colon ":" not a "*". The "which" argument in
TukeyHSD() must be identical to the name of the coefficient in the summary
tab
Hi,
I need to determinate HSD value from a matrix like that
Thesis Days A1
Cx 0 66.07
Cx 0 60.24
Cx 0 42.86
Tw 0 66.07
Tw 0 60.24
Tw 0 42.86
Aa 0 66.07
Aa 0 60.24
Aa 0 42.86
Qe 0 66
Mar 25, 2011; 12:58am Simon Bate wrote:
>> I've been happily using the TukeyHSD function to produce Tukeys HSD tests
>> but have decided to try
>> out Multcomp instead. However when I carry out the test repeatedly I have
>> found that Multcomp
>> produces slightly different values each time. (se
On 2010-12-13 07:29, PGZC wrote:
TukeyHSD(fit)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an object of class "lm"
TukeyHSD(anova)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an object of class
"function"
TukeyHSD
TukeyHSD(fit)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an object of class "lm"
> TukeyHSD(anova)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an object of class
"function"
> TukeyHSD(group)
Error in UseMethod("TukeyHS
TukeyHSD(fit)
Error in TukeyHSD(fit) : object 'fit' not found
> TukeyHSD(anova)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an object of class
"function"
> TukeyHSD(group)
Error in UseMethod("TukeyHSD") :
no applicable method for 'TukeyHSD' applied to an
On 2010-12-12 10:21, PGZC wrote:
DrugUS1 US2 Aptecha
Celebrex235.54 269.99 121.02
Detrol LA 157.99 190.99 55.3
Flomax 166.00 190.99 93.45
Lipitor 174.99 200.99 137.7
Novaldex108.6 129.99 22.48
Norvasc 186.66 203.99 161.93
Plavix 107.99 106.99
DrugUS1 US2 Aptecha
Celebrex235.54 269.99 121.02
Detrol LA 157.99 190.99 55.3
Flomax 166.00 190.99 93.45
Lipitor 174.99 200.99 137.7
Novaldex108.6 129.99 22.48
Norvasc 186.66 203.99 161.93
Plavix 107.99 106.99 64.53
Prevacid117.39 134.
install.packages('sos') # if you do not already have it
library(sos)
TH <- ???TukeyHSD
#found 27 matches; retrieving 2 pages
TH # The print method displays the results in a table in a web browser
This suggests you might also be interested in "multcompView".
Unfortunately, I failed to
Please, the language of the mailing list is English. A majority of
readers do not understand what you are talking about.
The Tukey test was designed for comparing means in equal replication
ANOVA. Some "folklore" suggests that it may work reasonably for
not-too-unequal replication cases, but I wou
Hallo,
zur Analyse von Daten zum Artenreichtum von Pflanzen, habe ich ein Glm (glm)
und anschlieÃend eine Anova (anova) durchgeführt. Nun möchte ich für die
signifikanten Einflussfaktoren einen Post Hoc Tukey Test durchführen, um zu
ermitteln in wie weit die einzelnen Fakt
Tukey's pairwise comparisons test is included in the multicomp package.
Also, see the
MMC function in the HH package, available from CSAN.
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
Hello,
I am a little out of date and am still using S-Plus instead of R. I haven't
been able to find the right place to ask this question, so I thought I would
ask it here and hope that someone can help. I am unable to located the
TukeyHSD() function in S-Plus. I have been able to use the GU
Hi Christina,
>> How can this happen? How can the p-values from the Tukey become
>> significant when the lme-model wasn't?
The link below, with an explanation by Prof. Fox is relevant to your
question:
http://www.nabble.com/Strange-results-with-anova.glm()-td13471998.html#a13475563
Another wa
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
Hi Ullrich,
>> # what does '~1 | Subj/Cond' mean?
It is equivalent to your aov() error structure [ ... +Error(Subj/Cond) ].
It gives you a set of random intercepts, one for each level of your nesting
structure.
## To get some idea of what's being done it helps to have a continuous
covariate in
Great, thanks a lot, guys!
Only thing is I now have two working versions, that yield *slightly*
different results.
ACCaov <- aov(Acc ~ Cond + Error(Subj/Cond), WMU3C)
ACClme <- lme(Acc ~ Cond, random = ~1 | Subj/Cond, WMU3C) # what does
'~1 | Subj/Cond' mean?
summary(glht(ACClme, linfct=mcp(Co
Hi Ullrich,
>> The model is
>> RT.aov <- aov(RT~Cond + Error(Subj/Cond), WMU3C)
>> I understand that TukeyHSD only works with an aov object, but that
>> RT.aov is an aovlist object.
You want to use lme() in package nlme, then glht() in the multcomp package.
This will give you multiplicity adju
Thanks, Dieter,
but as far as I understand, 'glht' does not support objects of class
'aovlist' either. I mean, I know there is a "TukeyHSD" function out
there, but that's the problem: repeated measures ANOVA yields an
aovlist object, and TukeyHSD calls for an aov object.
And I don't know if
Ullrich Ecker uwa.edu.au> writes:
> 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.
>
> Now, I want to calculate a post-hoc test following up a within-subjects ANOVA.
Prob
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.
Hello everybody!
I've just drawn an interaction plot between two factors, and I would like also
to draw the Tukey HSD intervals, but I only was able to draw sd arrows.
How I can draw that with R?
Thank you very much and many thanks for your advices.
Regards
___
35 matches
Mail list logo
