如何找到r中数据子集的方差？

【问题标题】：How to find the variance of a subset of the data in r?如何找到r中数据子集的方差？
【发布时间】：2020-07-09 09:14:24
【问题描述】：

我正在尝试查找整个数据 (dat) “pollutionData.csv”的子集的方差。

我想要下雨 = 0 时 PM2.5 水平的方差。

var(PM2.5, data=subset(dat, RAIN == 0))

上面的代码不起作用。

aggregate(dat[, 6], list(dat$RAIN==0), var, na.rm=TRUE)

上面的代码在rain = 0和rain > 0时输出方差，但我想对方差进行假设检验，所以这没有帮助。

任何帮助将不胜感激！

【问题讨论】：

【解决方案1】：

我们可以将“RAIN”为 0 的“PM2.5”作为子集，然后取var

with(dat, var(PM2.5[RAIN == 0], na.rm = TRUE))

另一种选择是replace 'PM2.5 中的值'RAIN' 为 0 和NA，然后应用var

with(dat, var(replace(PM2.5,  RAIN == 0, NA), na.rm = TRUE))

aggregate 在我们想要分组操作时是必需的。在这里，我们只得到了“RAIN”为 0 的“PM2.5”的var

【讨论】：