2011-02-08から1日間の記事一覧

メモ

R

こちら 確率分布の勝負(こちら) 離散事象の確率の掛け算 # 事象1:nの確率分布が2つある n<-6 library(MCMCpack) p1<-rdirichlet(1,rep(1,n)) p2<-rdirichlet(1,rep(1,n)) # (1:n)x(1:n)の同時確率は p12<-t(p1)%*%p2 win1<-sum(p12[lower.tri(p12)]) win2<-…

メモ

wiki記事 # 2変量フーリエ n<-2 nf<-10 rs<-runif(nf) rs<-sort(rs,decreasing=TRUE) ps<-runif(nf)*2*pi #ps<-rep(0,nf) Niter<-1000 dt<-1 t<-seq(from=0,to=1,length.out=Niter)*2*pi*10 xy<-matrix(0,Niter,n) for(i in 1:nf){ #xy[,1]<-xy[,1]+rs[i]*c…

メモ

R

チーム対決ソース(ポアソン分布で近似する) # チーム2の得点を動かして # 勝敗アイコ率をとる a1<-15 N<-1000 a2s<-seq(from=0,to=30,length.out=N) win1<-win2<-tie<-rep(0,N) for(i in 1:N){ a2<-a2s[i] p1<-dpois(0:100,a1) p2<-dpois(0:100,a2) ylim<-c…