こんな記事がありました。「席替えの数理」 # クラスの人数 ns <- 2^(1:10) # 席替え試行数 n.iter <- 10^4 # 席が変わらなかった人数の格納用行列 no.changers.mat <- matrix(0,length(ns),n.iter) for(j in 1:length(ns)){ n <- ns[j] ori <- 1:n no.chang…

こちらから Wiki ガンマ分布を見てみる は離散的に定義されているけれど、ガンマ関数を使って、その制約をとってやったものをポアソン分布の拡張とすれば「積分して１」になるための定数項を無視すれば l<-2 plot(function(x)l^x*exp(l)/gamma(x+1),from=0,t…