Re: [R] gamma distribution

From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>
Date: Wed 27 Jul 2005 - 17:18:26 EST

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

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