最大样本索引的分布

时间:2017-04-24 11:04:01

标签: r

我正在提取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)

这似乎很奇怪,谢谢你的任何建议!

1 个答案:

答案 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))

enter image description here