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
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 > >
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...{{dropped}}
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. > >------------------------------------------------------------- > ----------------- > > > >
Thanks. But that's interesting. ======= At 2006-03-16, 15:15:49 you wrote: ======>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. >> >------------------------------------------------------------- >> ----------------- >> >> >> >> > > >------------------------------------------------------------------------------ >Notice: This e-mail message, together with any attachments...{{dropped}}