Re: [R] TukeyHSD for multiple response

Michael, Thank you very much for your answer. I finally tried lsmeans to compare what I wanted. I'll follow your advice and explore the CDA. It's probably a better solution to assess what I want. Best, Sérgio. - Mensagem original - > De: "Michael Friendly" > Para:

Re: [R] TukeyHSD for multiple response

Hi Sergio Doing those tests 30 times is going to give you a huge Type I error rate, even if there was a function that did that. There is a reason why TukeyHSD doesn't make it easy. In general, if there are useful comparisons among the species, you are better off setting up and testing

[R] TukeyHSD for multiple response

Dear all, I'm testing the effect of species and sex in my sample by using the principal component scores of a PCA analysis. I have 30 PCs and I tried to see if there is any significant difference from males to females, given that there is a significant effect of phylogeny (factor with

[R] TukeyHSD and glht differ for models with a covariate

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

Re: [R] TukeyHSD troubles

Bart, I want to thank you for your code. I was having similar problems as Amy, even after setting my numeric variable as a factor using as.factor(). I used is.factor() to confirm and received the answer as TRUE from R; however after running the TukeyHSD() my set factor in my aov() was not read

[R] TukeyHSD output

Hi all, Is there a R-function that orders Tukey results with conveniant letters, similar to the SPSS output (A, AB, ABC, C, etc.) . [[alternative HTML version deleted]]

Re: [R] TukeyHSD output

On 2012-08-24 8:58, David Douterlungne wrote: Hi all, Is there a R-function that orders Tukey results with conveniant letters, similar to the SPSS output (A, AB, ABC, C, etc.) . library(multcompView)

Re: [R] TukeyHSD output

-project.org Subject: Re: [R] TukeyHSD output On 2012-08-24 8:58, David Douterlungne wrote: Hi all, Is there a R-function that orders Tukey results with conveniant letters, similar to the SPSS output (A, AB, ABC, C, etc.) . library(multcompView) ?multcompLetters

[R] TukeyHSD not working

Dear r-help members. I would like to compare species numbers of moths between eight different forests (each sampled for six nights). I would like to do a nested anova to compare species numbers between forests and nights. For more site specific details I wanted to do a Tukey test (TukeyHSD).

Re: [R] TukeyHSD not working

Michael, Use dput() to output your data (or perhaps a small subset). Then paste the result of that call and your R code (just those lines of code that are needed to reproduce the problem) right in your message to R-help. That makes it easier for the R-help list readers to help you

Re: [R] TukeyHSD not working

Michael Eisenring eimic...@ethz.ch wrote on 07/20/2012 09:35:03 AM: Dear Jean, thanks for this email. I'm grateful for this instruction Just to make sure that I understood you properly: something like this?: Michael, Yes, this is perfect. Very helpful. My data:

Re: [R] TukeyHSD not working

Hi, Check this link. http://stackoverflow.com/questions/5904513/tukeyhsd-after-within-factors-anova A.K. - Original Message - From: Michael Eisenring eimic...@ethz.ch To: r-help@r-project.org Cc: Sent: Friday, July 20, 2012 6:15 AM Subject: [R] TukeyHSD not working Dear r-help

Re: [R] TukeyHSD not working

HI, I don't have much experience with laercio package.  If it works for you, it is good. A.K. - Original Message - From: Michael Eisenring eimic...@ethz.ch To: arun smartpink...@yahoo.com Cc: Sent: Friday, July 20, 2012 8:38 AM Subject: Re: [R] TukeyHSD not working Hi Thanks

Re: [R] TukeyHSD troubles

Thanks! I was having the same issue. my treatments were 'days' so they were in intervals of 15 days. My anova was fine but I couldn't get my tukey to produce results. Used your code and it worked. -- View this message in context:

[R] TukeyHSD plot error

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

Re: [R] TukeyHSD plot error

Hmmm, I can't reproduce, but I'm not really sure why that would happen... is there any way you can test this in a --vanilla R session? (That's the UNIX-y way to start a totally clean session; not sure exactly how to achieve that on Windows) Does this happen if you just run example(TukeHSD)

Re: [R] TukeyHSD with Type III SS

Ali S wrote: Is there a way to run TukeyHSD with Type III instead of Type I SS? I am sure Bill Venables will respond, but in the meantime read: http://r-project.markmail.org/search/?q=type+III+venables And search google for Venables exegesis. It's one of the unique terms. Dieter --

[R] TukeyHSD with Type III SS

Is there a way to run TukeyHSD with Type III instead of Type I SS? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] TukeyHSD error

All - I think I'm being dense, but for the life of me, I can't figure out why I get error message with the code below. I have data that looks like param level perc.surv asin.tran DO 3 0.6864407 0.9764544 DO 3 0.125 0.3613671 DO 3 0.8738739 1.2077299 DO 4 0.4615385 0.7468986

Re: [R] TukeyHSD error

Rescinded. Problem solved. I discovered that TukeyHSD doesn't like numeric factors; changing the level to a factor solved my problem. SR Steven H. Ranney steven.ran...@montana.edu __ R-help@r-project.org mailing list

[R] TukeyHSD responses in R

Hi All, Does anyone know how to get contrast statements or an output similar to TukeyHSD when doing survival analysis with the coxph function? Thanks, Adriana -- View this message in context: http://r.789695.n4.nabble.com/TukeyHSD-responses-in-R-tp2530120p2530120.html Sent from the R help

[R] TukeyHSD for binomial distribution

Hi, I am looking for a way to get contrast statements (similar to the output of the TukeyHSD function) for glm fit data with a binomial distribution: glm.fit=glm(y~(treat)*(geno),family=binomial) tukey.fit=TukeyHSD(aov(glm.fit)) My y is a dataframe of number successful and number failed

[R] TukeyHSD model thing

Hi, I am trying to reproduce a tukey test in R == x=c(145,40,40,120,180, 140,155,90,160,95, 195,150,205,110,160, 45,40,195,65,145, 195,230,115,235,225, 120,55,50,80,45 ) y2=c( rep(as.character(1),5),

Re: [R] TukeyHSD model thing

On 3/6/2010 4:38 PM, casperyc wrote: Hi, I am trying to reproduce a tukey test in R == x=c(145,40,40,120,180, 140,155,90,160,95, 195,150,205,110,160, 45,40,195,65,145, 195,230,115,235,225, 120,55,50,80,45 ) y2=c(

[R] TukeyHSD troubles

I've tried to run a Tukey post-hoc but keep getting this weird error, whether the aov was significant or not. treat_code is a dummy variable, but that shouldn't matter. Any suggestions? Thanks Amy summary(aov(EtoH~treat_code, mydata)) Df Sum Sq Mean Sq F value Pr(F) treat_code

Re: [R] TukeyHSD troubles

treat_code is a dummy variable, but that shouldn't matter.  Any suggestions? It does matter to TukeyHSD. If treat_code is a numeric variable with discrete values 0 and 1, then it does not have class factor. It is true that aov will give the same ANOVA table for a two-level factor as for a

Re: [R] TukeyHSD troubles

treat_code isn't a factor, but a numeric variable. You should use: summary(aov(EtoH~as.factor(treat_code), mydata)) TukeyHSD(aov(EtoH~as.factor(treat_code), mydata)) Bart -- View this message in context: http://n4.nabble.com/TukeyHSD-troubles-tp1570205p1570228.html Sent from the R help

[R] TukeyHSD no longer working with aov output?

I can prove I've done this before, but I recently installed Rexcel (and it was easiest to reinstall R and some other bits to make it work) and now its no longer working. Before I would do an ANOVA and a tukey post-hoc like this: data1.aov=aov(result~factor1*factor2, data=data1) then...

Re: [R] TukeyHSD no longer working with aov output?

Hi Clayton, I don't think you need summary(). TukeyHSD(data1.aov) should work. -Ista On Tue, Oct 20, 2009 at 4:17 PM, Clayton Coffman clayton.coff...@gmail.com wrote: I can prove I've done this before, but I recently installed Rexcel (and it was easiest to reinstall R and some other bits to

Re: [R] TukeyHSD no longer working with aov output?

I've tried that as well, and I get an error like: TukeyHSD(phenolic.aov) Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 'rep' In addition: Warning messages: 1: In replications(paste(~, xx), data = mf) : non-factors ignored: time 2: In replications(paste(~, xx), data = mf) :

Re: [R] TukeyHSD no longer working with aov output?

Problem solved: Lesson learned (I think): TukeyHSD doesn't like it when you use numbers as names for the factors. ie factor time cannot be 24 and 48 but twenty-four and forty-eight work fine. On Tue, Oct 20, 2009 at 5:36 PM, Clayton Coffman clayton.coff...@gmail.comwrote: I've tried that as

Re: [R] TukeyHSD no longer working with aov output?

I think its a problem with my data, something about how Rexcel imported it We don't have enough information to be sure. My guess is that your data in Excel is integers which are intended to be levels of a factor. Excel doesn't distinguish between integers and integers that might be factor

Re: [R] TukeyHSD no longer working with aov output?

interesting to know, that might explain it. When I imported I just saved as CSV and imported with read.csv, I've never done anything to specify the integers as factors (this is the first time I've used numbers as names). Here are the precise commands I use(d): phen=read.csv(file=phenolics.csv)

Re: [R] TukeyHSD no longer working with aov output?

Now I am worried that you have a wrong analysis. the aov function is perfectly happy using either factors or numeric variables. Are there really only two levels of time, which is what one degree of freedom for time suggests? Or are there more than two level, but since aov() sees that as a

Re: [R] TukeyHSD no longer working with aov output?

There are two factors of time, but they are evenly replicated across all the other factors/levels. The experiment is perfectly balanced except for one lost sample, which is deleted automatically in the aov. I am very certain the analysis is correct. I think its merely a discrepancy between how

Re: [R] TukeyHSD no longer working with aov output?

I think that Richard may be correct. In my experience it's always a bad idea to use numerical labels for factors, especially when importing data from Excel. But why not do a str() on your data to see whether R thinks that time is a factor or not? And if not, why not convert time to factor before

[R] TukeyHSD with Cohens d ?

hi, is there a function or an easy way to get the effect size (cohens d) while performing multiple comparisons via the TukeyHSD function? thanks! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] TukeyHSD - t-value

Hi, how do I get the t-value when I perform multiple comparisons with the TukeyHSD() - function? thanks for any help! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] TukeyHSD and 0.0000000

On Wed, Nov 5, 2008 at 5:42 AM, Fredrik Karlsson [EMAIL PROTECTED] wrote: Dear list, Sorry to ask you this, but I just ran a TukeyHSD on an model with a two thee level factors as independent variables and a numeric score dependent variable. The aov gives a significant interaction effect, and

[R] TukeyHSD and 0.0000000

Dear list, Sorry to ask you this, but I just ran a TukeyHSD on an model with a two thee level factors as independent variables and a numeric score dependent variable. The aov gives a significant interaction effect, and using the TukeyHSD, I get almost every row to be exactly 0.000. Should I

Re: [R] TukeyHSD and 0.0000000

Hi, Sorry about that. Sure, here is some further information: summary(InekeReduc) segorg intorg l2group mani score Belfast:3782 Belfast:3782 1: 2604 1 :2480 Min. : 1.000 L2 :7502 L2 :7440 2:12462 2 :2542 1st Qu.: 1.000 German

Re: [R] tukeyHSD

, -Cody Cody Hamilton Edwards Lifesciences -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of hodgess Sent: Thursday, January 03, 2008 8:18 AM To: José Alberto Monteiro; [EMAIL PROTECTED] Subject: Re: [R] tukeyHSD Suppose you have an aov model x.aov

[R] tukeyHSD

Does anybody know how to deffine in the test tukeyHSD the pairs of comparison that you want to get? It is throw the lmat commnad, but I don't know the correct usage... Thanks a lot! José -- MSc José Alberto F. Monteiro Botanisches Institut Universität Basel السلام عليكم

Re: [R] tukeyHSD

: [R] tukeyHSD Does anybody know how to deffine in the test tukeyHSD the pairs of comparison that you want to get? It is throw the lmat commnad, but I don't know the correct usage... Thanks a lot! José -- MSc José Alberto F. Monteiro Botanisches Institut Universität Basel السلام عليكÙ

[R] TukeyHSD doubts

Hello, I have some doubts on TukeyHSD application. I want to investigate the effects of depth, latitude and month variation on the length of a fish. These are orthogonal and observational data. For this, I have made an aov model (L~month+lat+prof+month*lat), after applying drop1 and step

Re: [R] TukeyHSD doubts

This is a repost of https://stat.ethz.ch/pipermail/r-help/2007-September/141727.html Please do study the posting guide to see why you did not get an answer and what to do when you do not. It is nothing to do with TukeyHSD, as the differences are there in the means. It is clear that your

[R] TukeyHSD doubts

Hello, I have some doubts on TukeyHSD application. I want to investigate the effects of depth, latitude and month variation on the length of a fish. These are orthogonal and observational data. For this, I have made an aov model (L~month+lat+prof+month*lat), after applying drop1 and step