[R] problem for wtd.quantile()

Liaw, Andy andy_liaw at merck.com
Thu Mar 16 15:15:49 CET 2006


I don't know why using the fractional weights and normwt=TRUE gives that
answer, but:

> wtd.quantile(x, 5:1)
  0%  25%  50%  75% 100% 
   1    1    2    3    5 

Andy

From: Jing Yang 
> 
> Actually, I am looking for this result:
>   0%  25%  50%  75% 100% 
>    1    1    2    3    5
> -----------------------
> y <- c(1,1,1,1,1,2,2,2,2,3,3,3,4,4,5)
> quantile(y)
> -----------------------
> 
> Since x with weights w is equivalant to y.
> 
> ======= At 2006-03-16, 14:34:14 you wrote: =======
> 
> >Perhaps you're looking for this?
> >
> >> ?wtd.quantile
> >> wtd.quantile(x,weights=w, normwt=TRUE)
> >  0%  25%  50%  75% 100% 
> >   1    2    2    3    5 
> >
> >Andy
> >
> >From: Jing Yang
> >> 
> >> Dear R-users,
> >> 
> >> I don't know if there is a problem in wtd.quantile (from
> >> library "Hmisc"):
> >> --------------------------------
> >> x <- c(1,2,3,4,5)
> >> w <- c(0.5,0.4,0.3,0.2,0.1)
> >> wtd.quantile(x,weights=w)
> >> -------------------------------
> >> The output is:
> >>   0%  25%  50%  75% 100% 
> >> 3.00 3.25 3.50 3.75 4.00
> >> 
> >> The version of R I am using is: 2.1.0
> >> 
> >> Best,Jing
> >> 
> >> 
> >
> >
> >-------------------------------------------------------------
> ----------
> >-------
> >Notice:  This e-mail message, together with any attachments, 
> contains information of Merck & Co., Inc. (One Merck Drive, 
> Whitehouse Station, New Jersey, USA 08889), and/or its 
> affiliates (which may be known outside the United States as 
> Merck Frosst, Merck Sharp & Dohme or MSD and in Japan, as 
> Banyu) that may be confidential, proprietary copyrighted 
> and/or legally privileged. It is intended solely for the use 
> of the individual or entity named on this message.  If you 
> are not the intended recipient, and have received this 
> message in error, please notify us immediately by reply 
> e-mail and then delete it from your system.
> >-------------------------------------------------------------
> -----------------
> 
> 
> 
>




More information about the R-help mailing list