我正在尝试用R语言进行编程。
p_1_testing <- seq(from = 0.92, to = 1, by = 0.001)
testing_matrix <- matrix(0, nrow = 1000-52, ncol = length(p_1_testing))
for (i in 53:1000) {
cr_left_testing <- qbinom(0.025, i, 0.82)
cr_right_testing <- qbinom(0.975, i, 0.82)
for (j in 1:length(p_1_testing)) {
testing_matrix[i-52, j] <- pbinom(cr_left_testing-1, i, p_1[j]) +
(1-pbinom(cr_right_testing, i, p_1[j]))
}
}
它给出下面的输出。
> testing_matrix
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 0.9999982 NA NA NA NA NA NA NA NA NA
[2,] 1.0000000 NA NA NA NA NA NA NA NA NA
[3,] 1.0000000 NA NA NA NA NA NA NA NA NA
[4,] 1.0000000 NA NA NA NA NA NA NA NA NA
[5,] 1.0000000 NA NA NA NA NA NA NA NA NA
[6,] 1.0000000 NA NA NA NA NA NA NA NA NA
[7,] 1.0000000 NA NA NA NA NA NA NA NA NA
[8,] 0.9999999 NA NA NA NA NA NA NA NA NA
[9,] 0.9999999 NA NA NA NA NA NA NA NA NA
[10,] 1.0000000 NA NA NA NA NA NA NA NA NA
[11,] 1.0000000 NA NA NA NA NA NA NA NA NA
[12,] 1.0000000 NA NA NA NA NA NA NA NA NA
[,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] NA NA NA NA NA NA NA NA NA
[2,] NA NA NA NA NA NA NA NA NA
[3,] NA NA NA NA NA NA NA NA NA
[4,] NA NA NA NA NA NA NA NA NA
[5,] NA NA NA NA NA NA NA NA NA
[6,] NA NA NA NA NA NA NA NA NA
[7,] NA NA NA NA NA NA NA NA NA
[8,] NA NA NA NA NA NA NA NA NA
[9,] NA NA NA NA NA NA NA NA NA
[10,] NA NA NA NA NA NA NA NA NA
[11,] NA NA NA NA NA NA NA NA NA
[12,] NA NA NA NA NA NA NA NA NA
这样一直到最后
[ reached getOption("max.print") -- omitted 936 rows ]
但是我不知道为什么它对几乎所有条目都给出NA，甚至[1，]条目也不是我想要的答案...
例如，我希望我的代码完成如下操作
> qbinom(0.025, 53, 0.82)
[1] 38
> qbinom(0.975, 53, 0.82)
[1] 49
> pbinom(37, 53, 0.921) + (1-pbinom(49, 53, 0.921))
[1] 0.3887079
> pbinom(37, 53, 0.92) + (1-pbinom(49, 53, 0.92))
[1] 0.3785474
例如，我期望我的代码执行的操作是将testing_matrix [1，1]设置为0.3785474，并将testing_matrix [1,2]设置为0.3887079。
为什么我的代码不起作用？请帮帮我！谢谢！