# Mean vs. Trimmed mean in the normal distribution

In a simple experiment with the normal distribution in R I ran 500 iterations of a simulated normal distribution with N=100 each. For each iteration from the 500 iterations, I calculated both the mean and the trimmed mean with 20% trim (from each side), resulting in 500 values for each. Then, I compared the values of both with a boxplot:

It seems that the mean values are more "precise". I have managed to reproduce these results in almost all tries, and in the tries that I couldn't, the boxplot resulted in a similar plot for each.

This feels a bit counter-intuitive. I expected for it to be the other way around, since the 20% trim will remove results with high deviation. The only explanation I was able to think of for this observation is that the trim removes data that would otherwise "balance" the mean, though, it's not a formal explanation.

Would love some insights on this observation, thanks!

• It might help to contemplate what your procedure would do with much smaller sample sizes, such as $n=5$ or even $n=3$.
– whuber
Nov 4, 2017 at 19:54

For an exponential family like the Normal distribution, the sample average $\bar{x}$ is know to achieve the Cramér-Rao lower bound, that is the minimal possible variance among all unbiased estimators of the mean. It is thus no surprise that another estimator such as the trimmed mean is found to be more variable than $\bar{x}$.
If you look at say a $t_4$ distribution you can see some gain from trimming. If you look at a uniform on $(-k,k)$, you can see a cost from trimming (indeed, you'd be better still to average the trimmed-off values at some very small level of trimming than to use the mean).