250 R -- 複数の相関係数の同等性の検定と点推定値 青木繁伸 2002/01/14 (月) 00:24
250. R -- 複数の相関係数の同等性の検定と点推定値 青木繁伸 2002/01/14 (月) 00:24 |
eq.corr <- function(n, r) { if (any(n <= 3)) { stop("n <= 3") } k <- length(n) v <- n-3 z <- 0.5*log((1+r)/(1-r)) sv <- sum(v) svz <- sum(v*z) chi <- sum(v*z*z)-svz^2/sv p <- 1-pchisq(chi,k-1) result <- c(chi, k-1, p) name <- c("chi sq.", "d.f.", "P value") if (p > 0.05) { zm = svz/sv result <- c(result, (exp(2*zm)-1)/(exp(2*zm)+1)) name <- c(name, "Estimated r") } names(result) <- name result } n <- c(10, 16, 8, 29, 36) r<- c(0.658, 0.285, 0.569, 0.427, 0.374) eq.corr(n, r) |
● 「統計学関連なんでもあり」の過去ログ--- 017 の目次へジャンプ
● 「統計学関連なんでもあり」の目次へジャンプ
● 直前のページへ戻る