From: <pantd_at_unlv.nevada.edu>

Date: Thu 28 Jul 2005 - 08:42:51 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 08:47:33 2005

Date: Thu 28 Jul 2005 - 08:42:51 EST

Hi

I ran your code. I think it should give me the number of p values below 0.05
significance level (thats what i could understand from your code), but after
running your code there is neither any error that shows up nor any value that
the console displays.

thanks in advance

-dev.

Quoting Uwe Ligges <ligges@statistik.uni-dortmund.de>:

> pantd@unlv.nevada.edu wrote:

*>
**> > 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.
**>
**> He?
**> - log() transforming the data does not change the Wilcoxon statistics
**> (based on ranks!)!
**> - Why is this related to Excel?
**> - What are you going to show?
**>
**> I get
**>
**> erg <- replicate(1000, {
**> x<-rgamma(10, 2.5, scale = 10)
**> y<-rgamma(10, 2.5, scale = 10)
**> wilcox.test(x, y, var.equal = FALSE)$p.value
**> })
**> sum(erg < 0.05) # 45
**>
**> which seems plausible to me.
**>
**>
**> Uwe Ligges
**>
**>
**>
**> >
**> > --
**> > 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 08:47:33 2005

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