[Rd] meaning of "trim" in mean()
Duncan Murdoch
murdoch at stats.uwo.ca
Sun Oct 28 15:01:55 CET 2007
On 28/10/2007 4:43 AM, Patrick Burns wrote:
> If the sentence in question were amended to:
>
> Values of trim outside that range ...
>
> then I think it would rule out the misinterpretation of
> the sentence.
Good suggestion. I've changed it.
Duncan Murdoch
>
> Pat
>
>
> Prof Brian Ripley wrote:
>
>> There is only one _range_ mentioned, (0, 0.5). I don't see how you can
>> construe 'that range' to be a reference to anything other than (0, 0.5).
>>
>> And why do you suppose the description for argument 'trim' is referring to
>> 'values' of a different argument?
>>
>> It is telling you what happens for values of trim < 0 or > 0.5: that is
>> not information that it is appropriate to excise.
>>
>>
>> On Thu, 25 Oct 2007, Peter Dalgaard wrote:
>>
>>
>>
>>> Liaw, Andy wrote:
>>>
>>>
>>>> (I see this in both R-patched r43124 and R-devel r43233.)
>>>> In the Argument section of ?mean:
>>>>
>>>> trim the fraction (0 to 0.5) of observations to be trimmed from each
>>>> end of x before the mean is computed. Values outside that range are
>>>> taken as the nearest endpoint.
>>>>
>>>> Then in the Value section:
>>>>
>>>> If trim is non-zero, a symmetrically trimmed mean is computed with a
>>>> fraction of trim observations deleted from each end before the mean is
>>>> computed.
>>>>
>>>> The description in "trim" to me sounds like Windsorizing, rather than
>>>> trimming. Should that be edited?
>>>>
>>>>
>>>>
>>>>
>>> I think so:
>>>
>>>
>>>
>>>> x <- sort(rnorm(10))
>>>> mean(x,trim=.1)
>>>>
>>>>
>>> [1] -0.6387413
>>>
>>>
>>>> mean(x[2:9])
>>>>
>>>>
>>> [1] -0.6387413
>>>
>>>
>>>> mean(x[c(2,2:9,9)]) # Winsorizing
>>>>
>>>>
>>> [1] -0.6204222
>>>
>>> So yes, it is trimming, not Winsorizing, and the last sentence in the
>>> description of "trim" is misleading and should be, well..., trimmed.
>>>
>>>
>>>
>>>
>>
>>
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
More information about the R-devel
mailing list