F分布

メモ

上の記事でF分布についてメモした 少しANOVAを。こちらも参考 n.iter <- 1000 ret <- rep(0,n.iter) n.g <- 2 n.s <- 200 for(i in 1:n.iter){ X <- rnorm(n.s) P <- sample(1:n.g,n.s,replace=TRUE) s1 <- which(P==1) s2 <- which(P==2) ret[i] <- t.test(…

F分布

F分布はカイ二乗分布の比と関係がある # 2つのカイ二乗分布乱数を発生させる d1 <- 3 d2 <- 5 q1 <- rchisq(1000,d1) q2 <- rchisq(1000,d2) # カイ二乗分布乱数を自由度で割ったものの比がF統計量 f <- (q1/d1)/(q2/d2) # F統計量は構成元の二つの自由度の…