[R] gamma distribution

[pantd@unlv.nevada.edu]

Hi
I ran your code. I think it should give me the number of p values below 0.05
significance level  (thats what i could understand from your code), but after
running your code there is neither any error that shows up nor any value that
the console displays.


thanks in advance

-dev.

Quoting Uwe Ligges <ligges at statistik.uni-dortmund.de>:
> 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
>