きわめて個人的なR言語のメモ

2004-09-26Rで因子分析その3 二変量の分析

回帰直線入り散布図行列

二編量の関係を吟味してみる。

散布図行列を描く

普通に散布図行列を描くには

pairs(baseball)

で良いですが、ここではRjpWikiで公開されている

pairs用 回帰直線つきpanelを使わせて頂いて、回帰直線も表示したいと思います。

panel.lsfit <- function(x,y,...)
{
f <- lsfit(x,y)$coef
xx<- c ( min (x) , max(x) )
yy <- f["X"] * xx+ f["Intercept"]
lines(xx,yy,col="red")
points(x,y)
}
pairs(baseball, panel = panel.lsfit, pch = ".") 

相関行列を算出

ここでは青木先生の関数“my.cor”を使います。

my.cor(baseball)

で、上三角行列が「単相関係数」、下三角行列が「偏相関係数」、対角要素が「重相関係数」である相関行列が得られる。

            Var 1       Var 2       Var 3       Var 4       Var 5        Var 6
Var 1  0.49835624  0.40778275  0.04897438  0.12344932  0.10101434  0.013764419
Var 2  0.46190776  0.60631408 -0.07673645 -0.04679785 -0.30636509 -0.008015155
Var 3  0.15587627 -0.12310150  0.42688692 -0.33025436 -0.14101041 -0.209701240
Var 4  0.07122171 -0.08357778 -0.11436855  0.78567585  0.48191262  0.474239488
Var 5  0.21238691 -0.34848652 -0.02442644  0.35389118  0.61150159  0.259878898
Var 6 -0.05639442  0.04290625 -0.06561933  0.41223813  0.02322316  0.538118607
Var 7 -0.01568597  0.13425732  0.01626393  0.39812289 -0.18950813 -0.174814423
Var 8 -0.02568533 -0.13160480  0.21565671 -0.42459077 -0.05901787  0.070749878
Var 9 -0.04652111  0.19419433  0.03019263 -0.10770892  0.03879625  0.205311341
            Var 7      Var 8        Var 9
Var 1  0.09508995 -0.1386073  0.043771291
Var 2  0.19210838 -0.1167896  0.217348290
Var 3 -0.12014682  0.3644515 -0.009790575
Var 4  0.36857085 -0.5891611 -0.016409118
Var 5 -0.03273129 -0.2842107 -0.045642082
Var 6  0.03870287 -0.2441358  0.167593265
Var 7  0.48812646 -0.2525798  0.054702262
Var 8 -0.01260427  0.6341374 -0.032818085
Var 9  0.05793332 -0.0293333  0.303570602

トラックバック - http://r-user.g.hatena.ne.jp/bob3/20040926