[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