[過去ログ]
臨床統計もおもしろいですよ、その1 [無断転載禁止]©2ch.net (747レス)
臨床統計もおもしろいですよ、その1 [無断転載禁止]©2ch.net http://egg.5ch.net/test/read.cgi/hosp/1493809494/
上
下
前次
1-
新
通常表示
512バイト分割
レス栞
抽出解除
必死チェッカー(本家)
(べ)
レス栞
あぼーん
このスレッドは過去ログ倉庫に格納されています。
次スレ検索
歴削→次スレ
栞削→次スレ
過去ログメニュー
202: 卵の名無しさん [sage] 2017/10/27(金) 15:54:12 ID:dzxKDqmi # n発r中の狙撃手がN発狙撃するときの命中数を返す Golgo.sim <- function(N, n, r, k=10^3,Print=TRUE){ # k:シミュレーション回数 f <-function(S,.N=N,.n=n){ # 成績サンプル:命中S個、外れ(N-S)個 y=c(rep(1,S),rep(0,.N-S)) sum(sample(y,.n)) # その成績サンプルからn個数取り出したときの命中数 } xx=r:N # r未満ではr個命中することはないので除外 SS=NULL # 容れ子 for(i in 1:k){ x=sapply(xx,f) # 命中数の配列 SS=append(SS,which(x==r)-1+r) # 命中数がr個のときの成績サンプルの命中数Sの配列をつくる } print(summary(SS)) print(quantile(SS,probs=c(0.025,0.05,0.95,0.975))) print(c(mode=names(which.max(table(SS)))),quote=FALSE) if(Print) { hist(SS,xlim=c(0,N),freq=FALSE,col=sample(colors(),1),main='',xlab='Hits') lines(density(SS))} invisible(SS) } http://egg.5ch.net/test/read.cgi/hosp/1493809494/202
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.110s*