From: Christoph Buser <buser_at_stat.math.ethz.ch>

Date: Thu 28 Jul 2005 - 18:27:01 EST

Christoph Buser <buser@stat.math.ethz.ch> Seminar fuer Statistik, LEO C13

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

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 Thu Jul 28 18:30:56 2005

Date: Thu 28 Jul 2005 - 18:27:01 EST

Hi

Again to come back on the question why you don't get identical p.values for the untransformed and the transformed data.

I ran your script below and I get always 2 identical test per loop. In your text you are talking about the first 1000 values for the untransformed and the next 1000 values for the transformed.

But did you consider that in each loop there is a test for the
untransformed and the transformed, so the tests are printed
alternating.

This might be a reason why you did not get equal results.

Hope this helps,

Christoph

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 Thu Jul 28 18:30:56 2005

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