[R] problem for wtd.quantile()
Jing Yang
yangj at student.ethz.ch
Thu Mar 16 15:29:11 CET 2006
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}}
More information about the R-help
mailing list