- 930 名前：卵の名無しさん mailto:sage [2018/08/29(水) 13:30:50.73 ID:FkipPE46.net]
- >>887
予想通り、
数学板では負の二項分布の特殊な場合である幾何分布の期待値の公式からの答が、即答された。
寄付金の期待値をだすシミュレーションスクリプトを書いてみた。
invest <- function(p=0.30,q=3000){
win=sample(1:0,1,prob=c(p,1-p))
try=1
while(!win){
win=sample(1:0,1,prob=c(p,1-p))
try=try+1
}
return(try*q)
}
k=1e5
re.3=replicate(k,invest(0.30,3000))
hist(re.3,freq=FALSE,col='lightblue')
summary(re.3)
re.6=replicate(k,invest(0.60,8000))
hist(re.6,freq=FALSE,col='lightblue')
summary(re.6)
十万回のシミュレーションでの寄付金の値の分布のサマリーがこれ
> summary(re.3)
Min. 1st Qu. Median Mean 3rd Qu. Max.
3000 3000 6000 10025 12000 96000
> summary(re.6)
Min. 1st Qu. Median Mean 3rd Qu. Max.
8000 8000 8000 13296 16000 56000
|
 |
