我正在提取10个独立的二项分布随机变量Xi的样本,然后找到其中最大的索引。我认为每个都可能是最大值,但直方图显示了一个独特的模式,继承人histogram
我在R中使用了以下代码:
n=10^5
Which_maximum = function(x){
sample1 = rbinom(10,5,0.5)
return(which.max(sample1))
}
repeat1 = sapply(1:n,Which_maximum)
hist(repeat1)
这似乎很奇怪,谢谢你的任何建议!
答案 0 :(得分:0)
问题在于which.max
确定数字向量的第一个最大值的位置。
例如,尝试以下代码:
set.seed(23456)
x <- rbinom(10,5,0.5)
print(x)
[1] 2 3 3 4 3 2 3 4 4 3
which.max(x)
[1] 4
在向量x
中,有3个绑定的最大值,但which.max
仅为您提供第一个的位置。当然,这个问题会修改最大值索引的分布。
或者,您可以使用which.is.max
包的nnet
命令,其中“在向量中找到最大位置,随意断开关系”:
library(nnet)
n <- 10^5
which_maximum <- function(x) {
sample1 <- rbinom(10,5,0.5)
return(which.is.max(sample1))
}
repeat1 <- sapply(1:n,which_maximum)
barplot(table(repeat1))