[R] gamma distribution
Uwe Ligges
ligges at statistik.uni-dortmund.de
Wed Jul 27 09:18:26 CEST 2005
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
More information about the R-help
mailing list