忽略置信区间的 NA (CI)

Question

我想获得 R 中向量的置信区间，但我有一个 NA 值并且想忽略 NA。我有 na.rm = T 但没用。

这是我的向量和我尝试过的：

myvector <- c(60.558,60.680,60.757,60.750,60.708, 60.738, NA, 61.103,60.702,60.436)
CIvector <- CI(myvector,na.rm=T)

我收到一条错误消息：

'Error in 1 - ci : non-numeric argument to binary operator'

据我所知是因为 CI 函数不支持 NA.RM。如果我不使用 NA,RM 我只对所有区间值使用 NA

有人知道解决方法吗？

Answer 1

只需将向量包裹在 na.omit():

> CIvector <- CI(na.omit(myvector))
> CIvector
   upper     mean    lower 
60.85270 60.71467 60.57663

Answer 2

使用complete.cases。 Rmisc::CI

中没有 na.rm 参数

Rmisc::CI(myvector[complete.cases(myvector)])
#   upper     mean    lower 
#60.85270 60.71467 60.57663

Ignoring NA for confidence interval (CI)