Re: [R] multiple comparison of interaction of ANCOVA

2011-12-13 Thread Jinsong Zhao
Thank you very much for the suggestions. On 2011-12-13 3:20, Richard M. Heiberger wrote: ... On Mon, Dec 12, 2011 at 7:00 AM, Jinsong Zhao jsz...@yeah.net mailto:jsz...@yeah.net wrote: Apart from the data set here, is there a way to do multiple comparison on the interaction of one way

Re: [R] multiple comparison of interaction of ANCOVA

2011-12-12 Thread Jinsong Zhao
On 2011-12-11 22:49, Bert Gunter wrote: Inline below. -- Bert Graph the data sensibly to figure out what's going on. Statistical machinationsand anova tables with P values alone are not sufficient and can be opaque or misleading. If you do not know what sensibly is (or even if you do),

Re: [R] multiple comparison of interaction of ANCOVA

2011-12-12 Thread Jinsong Zhao
On 2011-12-12 1:48, Richard M. Heiberger wrote: Thank you for you use of HH. I think the right graph for this data is the much simpler ancova function library(HH) ancova(y ~ year * Trt, data=mydata) where we see that the three treatments have totally different slopes. Thank you very much for

Re: [R] multiple comparison of interaction of ANCOVA

2011-12-12 Thread Richard M. Heiberger
... On Mon, Dec 12, 2011 at 7:00 AM, Jinsong Zhao jsz...@yeah.net wrote: Apart from the data set here, is there a way to do multiple comparison on the interaction of one way analysis of covariance? Interaction in ancova usually means the slopes are different. ... BTW, when loading HH

[R] multiple comparison of interaction of ANCOVA

2011-12-11 Thread Jinsong Zhao
Hi there, The following data is obtained from a long-term experiments. mydata - read.table(textConnection( +y year Trt + 9.37 1993 A + 8.21 1995 A + 8.11 1999 A + 7.22 2007 A + 7.81 2010 A +10.85 1993 B +12.83 1995 B +13.21 1999 B +

Re: [R] multiple comparison of interaction of ANCOVA

2011-12-11 Thread Bert Gunter
Inline below. -- Bert On Sun, Dec 11, 2011 at 4:15 AM, Jinsong Zhao jsz...@yeah.net wrote: Hi there, The following data is obtained from a long-term experiments. mydata - read.table(textConnection( +        y year Trt +     9.37 1993   A +     8.21 1995   A +     8.11 1999   A +    

Re: [R] multiple comparison of interaction of ANCOVA

2011-12-11 Thread Richard M. Heiberger
Thank you for you use of HH. I think the right graph for this data is the much simpler ancova function library(HH) ancova(y ~ year * Trt, data=mydata) where we see that the three treatments have totally different slopes. The WoodEnergy example doesn't apply here. The WoodEnergy example