2007-10-27 メモ カイ自乗 > r1<-runif(100) > r2<-runif(100) > c1<-qchisq(r1,1) > c2<-qchisq(r2,1) > d1<-dchisq(c1,1) > d2<-dchisq(c2,1) > dmult<-d1*d2 > c12<-c1+c2 > d12<-dchisq(c12,2) > x12<-d12/(pi*sqrt(c1*c2)) > plot(dmult,x12)