From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>

Date: Wed 27 Jul 2005 - 17:18:26 EST

sum(erg < 0.05) # 45

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:22:14 2005

Date: Wed 27 Jul 2005 - 17:18:26 EST

> 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 Wed Jul 27 17:22:14 2005

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