[EMAIL PROTECTED] 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