2011-09-15から1日間の記事一覧

仮置き場

R

n1<- 20 n2<- 20 t_mean<- 6 t_sd<- 0.8 t_death<- 6 Niter<- 3 p_values<- rep(0,Niter) G1<- array(0,c(n1,2,Niter)) G2<- array(0,c(n2,2,Niter)) p1<- c(0.5,0.5) p2<- c(0.5,0.5) for(i in 1:Niter){ s1<- sample(0:1,size=n1,replace=TRUE) s2<- sampl…