分布を重ねる

  • 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))
hist(Xcis)
hist(Xtrans)
plot(acis,sort(Xcis),type="b",xlim=c(0,1),ylim=c(0,max(Xall)),col=2)
par(new=TRUE)
plot(atrans,sort(Xtrans),type="b",xlim=c(0,1),ylim=c(0,max(Xall)))

ecdf.cis <- ecdf(Xcis)
ecdf.trans <- ecdf(Xtrans)
plot(ecdf.cis)
plot(ecdf.trans,add=TRUE)
par(mfcol=c(1,1))