サンプル数の最大値の分布

Question

私は、10個の独立した二項分布ランダム変数Xiのサンプルを引っ張って、それらのうち最大のもののインデックスを見つける。私はRに次のコードを使用HERESにhistogram

、しかし、ヒストグラムは異なるパターンを示し、それぞれが最大であることが均等にそうだろうと思った：

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) 

これは本当に奇妙に思える、任意の提案に感謝！

Answer 1

which.maxは、の数値ベクトルの最大値の位置を決定するという問題があります。ベクトルxで

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 

は3つの結ばれた最大値が存在するが、which.maxはあなたの最初の位置のみを与える： たとえば、このコードを試してみてください。もちろん、この問題は、最大値のインデックスの分布を変更します。

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