Re: [R] gamma distribution

From: Christoph Buser <buser_at_stat.math.ethz.ch>
Date: Wed 27 Jul 2005 - 17:01:50 EST

Hi

I am a little bit confused. You create two sample (from a gamma distribution) and you do a wilcoxon test with this two samples. Then you use the same monotone transformation (log) for both samples and redo the wilcoxon test.
But since the transformations keeps the order of your samples the second wilcoxon test is identical to the first one:

x<-rgamma(10, 2.5, scale = 10)
y<-rgamma(10, 2.5, scale = 10)
wilcox.test(x, y, var.equal = FALSE)
x1<-log(x)
y1<-log(y)
wilcox.test(x1, y1, var.equal = FALSE)

Maybe you can give some more details about the hypothesis you'd like to test.

Regards,

Christoph Buser



Christoph Buser <buser@stat.math.ethz.ch> Seminar fuer Statistik, LEO C13
ETH (Federal Inst. Technology)	8092 Zurich	 SWITZERLAND
phone: x-41-44-632-4673		fax: 632-1228

http://stat.ethz.ch/~buser/

pantd@unlv.nevada.edu writes:
> Hi R Users
>
>
> This is a code I wrote and just want to confirm if the first 1000 values are raw
> gamma (z) and the next 1000 values are transformed gamma (k) or not. As I get
> 2000 rows once I import into excel, the p - values beyond 1000 dont look that
> good, they are very high.
>
>
> --
> sink("a1.txt");
>
> for (i in 1:1000)
> {
> x<-rgamma(10, 2.5, scale = 10)
> y<-rgamma(10, 2.5, scale = 10)
> z<-wilcox.test(x, y, var.equal = FALSE)
> print(z)
> x1<-log(x)
> y1<-log(y)
> k<-wilcox.test(x1, y1, var.equal = FALSE)
> print(k)
> }
>
> ---
> any suggestions are welcome
>
> thanks
>
> -devarshi
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jul 27 17:07:35 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:34:02 EST