R错误：“要替换的项目数不是替换长度的倍数”

Question

我对 R 很陌生。当值设置得较低时，我的脚本可以正常工作。我需要计算概率并且需要至少 1,000,000 个样本。我的第一个循环返回此错误：

Error in new.matrix[i, ] <- as.vector(table(million.rep[, i])) : number of items to replace is not a multiple of replacement length

我的脚本：

actual.prob <- c(.14, .14, .16, .13, .19, .24)
number.in.bag <- c(8,6,20,9,10,5)

million.rep <- replicate(1000000, sample(1:6, 58, replace= T, actual.prob))



new.matrix <- matrix(nrow= 1000000, ncol=6)
# Note to other readers... the above line throws the error.

for(i in 1:1000000){
new.matrix[i,] <- as.vector(table(million.rep [, i]))
}


y <- c()
for(i in 1:1000000){
y[i] <- sum(new.matrix[i,] == number.in.bag)
}
y
sum(y == 6)

Answer 1

table 计算所有出现的元素。如果缺少任何1:6，则输出的长度将小于 6，例如：

set.seed(1)

## first call; table output is length == 5
table(sample(1:6, size=10, replace=TRUE))
1 2 3 4 6 
1 2 1 3 3 

## second call; table output is length == 4
table(sample(1:6, size=10, replace=TRUE))
2 3 5 6
2 3 4 1

为了规避这种行为，您可以使用 factors 和定义的 levels (factor(sample, levels=1:6))：

## third call using factors; table output is always length == 6
## please notice the "levels=1:6"
table(factor(sample(1:6, size=10, replace=TRUE), levels=1:6))
1 2 3 4 5 6 
2 2 3 1 0 2 

对于您的示例，您可以使用：

new.matrix[i,] <- as.vector(table(factor(million.rep [, i]), levels=1:6))