[R] gamma distribution

[Uwe Ligges]

pantd at 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 at 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