2010-03-23から1日間の記事一覧

RでFDR(BH法)

R FDR

こちらも #まず、pの小さめなデータを作る Niter<-1000 library(Rassoc) st<-rep(0,Niter) p<-rep(0,Niter) for(i in 1:Niter){ af<-runif(1)*0.6+0.2 delta<-rnorm(1) af1<-af+af*0.05*delta af2<-af-af*0.05*delta case<-sample(c(0,1,2),1000,c(af1^2,2*a…