[R] quantiles and dataframe

[Gabor Grothendieck]

Sorry, try this instead.  It creates a data frame of 3 columns in which
each column equals RQ[,1] except it has NAs where the columns of RQ[-1]
has NAs.  Perform the quantile operation on that.

sapply(sign(RQ[-1]) * RQ[,1], quantile, probs = c(0, .05, .95, 1), na.rm = TRUE)


On 9/14/07, Gabor Grothendieck <ggrothendieck at gmail.com> wrote:
> Try this:
>
> sapply(RQ[-1], quantile, probs = c(0, .05, .95, 1), na.rm = TRUE)
>
>
>
> On 9/14/07, Anders Bjørgesæter <anders.bjorgesater at bio.uio.no> wrote:
> > Hi
> >
> > I have a dataframe, RQ, like this:
> >
> > A    B1    B2    B3
> > 1    NA    112    12
> > 2    NA    123   123
> > 3    NA    324    13
> > 4    3     21    535
> > 5    4     12    33
> > 6    7     1     335
> > 7    4     NA    3535
> > 8    4     NA    NA
> > 9    NA    NA    NA
> > 10    5    NA    NA
> > 12    4    NA    NA
> > 15    2    NA    NA
> > 17    3    NA    1
> > 63    1    NA    1
> > 75    NA   NA    NA
> > 100   NA   NA    NA
> > 123   NA   NA    NA
> > 155   NA   NA    NA
> > 166   NA   NA    NA
> > 177   NA   NA    NA
> >
> > I want to extract min, max, 5% and 95% from A based on the range of the Bs.
> >
> > Using this:
> >
> > s1<-A[min(which(!is.na(B1))):max(which(!is.na(B1)))]
> > q1<-quantile(s1,probs=c(0,5,95,100,NA)/100)
> >
> > I manage to get this by changing the B1 manually for each B
> >
> > B1    B2        B3
> > 4.0    1.00     1.00    (min)
> > 63.0   6.00     63.00   (max)
> > 4.5    4.5      1.65    (5%)
> > 40.0   6.00     63.00   (95%)
> >
> > I tried to use apply like this: s1<-apply(RQ,2,function(x)
> > {A[min(which(!is.na(RQ[,2:4]))):max(which(!is.na(RQ[,2:4])))] })
> >
> > to get the range of each B but that doesn't work.
> >
> > Also as you see, s1 includes the A where the B's are NA, e.g. for B1 I
> > get the 9 at row 9 (4,5,6,7,8,9,10,12,15,17,63) and not
> > (4,5,6,7,8,10,12,15,17,63), which I would prefer.
> >
> > BUT the main question is how can I extract min, max etc. from each B in
> > dataframe RQ without using a loop?
> >
> > Any help is greatly appreciated!
> >
> > Best Regards
> > Anders
> >
> > ______________________________________________
> > R-help at r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> >
>