列出置信区间的实际均值为 10

Question

我有一个table的Confident Interval如下图，题中有多少个Interval包含10的实际均值，并列出来。

我只知道数多少

Inside <- 0
for (i in 1:100)
{
if ((CIs[i,1]<10) & ( 10 < CIs[i,2]))
 {
 Inside <- Inside+1
 }
 else {}
}
Inside

你能帮我把它们列成矩阵吗？谢谢！

Answer 1

subset(CIs, lcl > 10 & ucl < 10)

这应该列出 lcl 大于 10 且 ucl 小于 10 的所有行。您需要将 "CIs" 替换为您的 data.frame 的名称。

计数可以比使用更容易：

nrow(subset(CIs, lcl > 10 & ucl < 10))

Answer 2

不清楚你想要什么。有一种可能。

# simulates a 'CIs' matrix
set.seed(666)
lwr <- rpois(20, 5)
upr <- lwr + rpois(20, 5)
CIs <- cbind(lwr = lwr, upr = upr)

# find intervals containing 10
contains <- apply(CIs, 1, function(lu) lu[1] < 10 && lu[2] > 10)

cbind(CIs, good = contains)
#       lwr upr good
#  [1,]   7  17    1
#  [2,]   3   8    0
#  [3,]  10  11    0
#  [4,]   3   6    0
#  [5,]   4  11    1
#  [6,]   6  10    0
#  [7,]  10  12    0
#  [8,]   5  10    0
#  [9,]   1   6    0
# [10,]   3  10    0
# [11,]   7  12    1
# [12,]   1   5    0
# [13,]   2   9    0
# [14,]   3   8    0
# [15,]   3   9    0
# [16,]   7  12    1
# [17,]   1   5    0
# [18,]   8  16    1
# [19,]   5  12    1
# [20,]   5  10    0

列出置信区间的实际均值为 10

List out the Confident Interval have actual means of 10

r

confidence-interval