2013-09-04から1日間の記事一覧

分布を重ねる

N数が違う分布が似ているかどうか、目で見る Ncis <- 10000 Ntrans <- 100 Xall <- rchisq(Ncis+Ntrans,3,5) Xcis <- Xall[1:Ncis] Xtrans <- Xall[(Ncis+1):(Ncis+Ntrans)] acis <- ppoints(Ncis,0.5) atrans <- ppoints(Ntrans,0.5) par(mfcol=c(2,2)) his…

分布を重ねる