カイ二乗分布

d01 <- rnorm(n=100000,mean=0,sd=1)
d02 <- rnorm(n=100000,mean=0,sd=1)
d03 <- rnorm(n=100000,mean=0,sd=1)
d04 <- rnorm(n=100000,mean=0,sd=1)

dx <- d01^2 + d02^2 + d03^2 + d04^2
hist(dx,breaks=20)
table(cut(dx,breaks=c(seq(0,20,0.1))))