昨日、検出力と非心カイ二乗分布のことを書いた 中で使用している"RandomSphere()"はこちら Ntrial<-100 Ma<-30 as<-seq(from=0,to=Ma,length.out=Ntrial) # 検定の有意水準はa^2できることにする N<-10000 # 乱数点の数指定 k<-sample(2:10,1) # 次元を指定…