From: <pantd_at_unlv.nevada.edu>

Date: Thu 28 Jul 2005 - 05:24:43 EST

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 05:31:46 2005

Date: Thu 28 Jul 2005 - 05:24:43 EST

Hi

You are right but here I am taking into account the p values I get from the
tests on the raw and the transformed samples. And then I calculate the power of
the tests based on the # of rejections of the p values.
DO you think its a good way to determine the power of a test?

thanks

-dev

Quoting Christoph Buser <buser@stat.math.ethz.ch>:

*> 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 Thu Jul 28 05:31:46 2005

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