>>612
このモデルでのシミュレーション結果

https://i.imgur.com/TlJotnB.png

# 観察された台数が不明なときのシミュレーション
sim <- function(){
M=m=0 # m:観察台数 M:最大番号 (初期値0)
while(M!=60){　 # M=60でないなら
N=sample(60:100,1) # タクシー総数Nを60 ~ 100から選ぶ
m=sample(1:min(N,60),1)# 観察する台数mを1 ~ min(N,60)から選ぶ
M=max(sample(1:N,m)) # N台からm台選択して最大値をMにいれる
}
return(N) # タクシー総数を返す
}
re=replicate(1e4,sim()) # 1万回繰り返して平均値(期待値)を算出
mean(re)