Re: [R] gamma distribution

From: <pantd_at_unlv.nevada.edu>
Date: Thu 28 Jul 2005 - 08:42:51 EST

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@statistik.uni-dortmund.de>:
> 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 Thu Jul 28 08:47:33 2005

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