这是我目前正在运行的代码:
n <- 7
N <- 52
r <- 13
reps <- 1000000
deck <- rep(c('h','d','c','s'), each = r)
diamonds <- rep(NA, length.out = reps)
pos <- sample(x = 1:52, size = 7, replace = FALSE)
for(i in 1:reps) {
hand <- sample(x = deck, replace = FALSE)[pos]
diamonds[i] <- sum(ifelse(hand == 'd', 1, 0))
}
barplot(table(diamonds), col = 'red', xlab = '# of diamonds',
ylab = paste('frequency out of',reps,'trials'),
main = paste('Positions:',pos[1],pos[2],pos[3],pos[4],
pos[5],pos[6],pos[7]))
我真正喜欢的是能够通过以下内容为barplot提供标题
barplot(..., main = paste('Positions:',pos))
并标题为“Positions:p1 p2 p3 p4 p5 p6 p7”,其中p1,p2,...,p7是pos的元素。
对于任何有兴趣的人,此代码随机从52中选择7个位置,然后在每次洗牌后的7个位置内计算钻石数量('d')1000000 shuffle。然后绘制了这7张卡片中钻石数量的经验分布。
答案 0 :(得分:3)
使用collapse
中的paste
折叠包含基本测试和pos
的矢量中的多个元素,
paste(c('Positions:', pos), collapse=" ")
否则,当您将“位置:”粘贴到pos
时,您会将前者回收到pos
的长度。