R aggregate 求不同组平均数，如何处理缺失值

今天碰到的一个问题是，用 aggregate 求不同组平均数的时候，缺失值不好处理。下面通过用 iris 这个经典数据说一下如果处理缺失值。

iris[sample(nrow(iris),5),] # 随即抽几行，大致浏览数据格式

## Sepal.Length Sepal.Width Petal.Length Petal.Width Species

136 7.7 3.0 6.1 2.3 virginica

85 5.4 3.0 4.5 1.5 versicolor

33 5.2 4.1 1.5 0.1 setosa

50 5.0 3.3 1.4 0.2 setosa

48 4.6 3.2 1.4 0.2 setosa

<code …