「統計学関連なんでもあり」の過去ログ--- 040

No.02887　変数間の線形関係の強さ　　【suzuki】　2007/03/01(Thu) 13:09

「x, y間の線形関係」と「x, y^2間の線形関係」の何れが強いのかを検定したいと考えています。以下のように「同じサンプルからの相関係数の差」（http://aoki2.si.gunma-u.ac.jp/R/diff_r.html）の検定を行えばよろしいでしょうか？
x = c(0.44, 0.51, 0.34, 0.46)
y = c(1.95, 2.25, 1.79, 2.72)
diff.r(4, cor(x,y), cor(x,y^2), cor(y,y^2))

以下，エラーになった箇所をご報告します。
http://aoki2.si.gunma-u.ac.jp/R/diff_r2.html
【修正前】rtest.kano <- function(x)
【修正案】diff.r2 <- function(x)
http://aoki2.si.gunma-u.ac.jp/R/diff_r.html
【修正前】
t0 <- (abs(rxy-rvy)*sqrt((n-1)*(1+rxv)))/s # 検定統計量
qrt(2*detR*(n-1)/(n-3)+(rxy+rvy)^2*(1-rxv)^3/4)
【修正案】
t0 <- (abs(rxy-rvy)*sqrt((n-1)*(1+rxv)))/ # 検定統計量
sqrt(2*detR*(n-1)/(n-3)+(rxy+rvy)^2*(1-rxv)^3/4)

No.02892　Re: 変数間の線形関係の強さ　　【青木繁伸】　2007/03/01(Thu) 14:07

ありがとうございます。
コメントを付ける作業をしているときにチョンボをしてしまったのですね。
修正しました。