多項分布

多項格子

となる整数ベクトルのとる格子点を再帰的な関数で作る 計算させてみないで作る方がよいのでこれこそHaskellが得意そうだが grid.multinom<-function(L,k){ ret<-list() if(k==1){ ret[[1]]<-c(L) }else{ cnt<-1 for(i in 0:L){ tmp<-grid.multinom(L-i,k-1) …

多項分布

こちらから ヒットを「0ヒット」「1ヒット」「2ヒット」…「ヒット」とする その確率をとする Vnumは打者人数。Nout.numは1イニング変わるためのアウトの数。 最後の打者は「0ヒット」 それ以外の打者Vothers=Vnum-1人の内訳は、「kヒット」の人数をとす…