[R] gamma distribution
pantd@unlv.nevada.edu
pantd at unlv.nevada.edu
Thu Jul 28 00:42:51 CEST 2005
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
>
More information about the R-help
mailing list