替換概率向量中的NA

Question

在R中概率向量中排除NA值的一種好方法是什麼？我正在做一些抽樣，有時在概率向量中有值，因爲總體中沒有特定值的出現。下面的代碼說明了問題：

population <- c("A", "B", "B") 
probs <- prop.table(table(population)) 
action <- c("A", "C", "D") 
test <- sample(action,1,prob=c(probs[1], probs[2], probs[3])) 

Answer 1

這就像你要指定使用​​相應元素的概率action每個元素的概率在我看來。您可以通過文字索引probs載體，具體做到這一點，索引probs與​​：

set.seed(1L); 
population <- c('A','B','B'); 
probs <- prop.table(table(population)); 
action <- c('A','C','D'); 
test <- sample(action,1L,prob=probs[population]); 
test; 
## [1] "C" 

您的代碼失敗，因爲probs只有兩個元素; probs[3]返回NA。還請注意probs[1:3]是c(probs[1],probs[2],probs[3])的更簡潔的因子分解。