Укажите Уверенный Интервал, у которого есть фактические средства 10

Question

У меня есть таблица доверительного интервала, как показано на рисунке ниже, вопросы о том, сколько интервалов содержит фактическое среднее значение 10, и перечислите их.

Я знаю только, как посчитать, сколько

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

Не могли бы вы помочь мне перечислить их в виде матрицы? Спасибо!

Answer 1

Не ясно, чего ты хочешь. Вот возможность.

# 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

Answer 2

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

в этом списке должны быть перечислены все строки с lcl больше 10 и ucl меньше 10. Вам нужно переставить "CI" с именем вашего data.frame.

Подсчет может быть намного проще, чем:

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