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

No.20373　相関関係の視覚化　　【コロン】　2013/10/25(Fri) 09:23

いつもお世話になっております。

また基本的なことかも知れませんが，ご指導いただきたいことがございます。

3変数の「相関関係」をpathDiagramを使って図示したいと思っております。
しかしその前のsem関数で解を求めることに苦しんでおります。
片矢印のパスを使ったパス解析については理解しておりますが，すべて両矢印の相関関係を表す際にどのように式を設けるべきなのかわかりません。ご指導頂けますでしょうか。

例は，豊田秀樹『共分散構造分析』（入門編）のp.43に掲載されているものです。

以下途中まで行ったものです。分散共分散などをどのように設定すればよろしいのでしょうか？

N=100です。

よろしくお願いいたします。

**************************************************
r<-readMoments(names=c("乳製品", "総熱量", "直腸ガン"))
1.00
0.65 1.00
0.20 0.58 1.00

model<-specifyModel()
乳製品<->総熱量, r1, NA
総熱量<->直腸ガン, r2, NA
直腸ガン<->乳製品, r3, NA

No.20375　Re: 相関関係の視覚化　　【青木繁伸】　2013/10/25(Fri) 10:17

推定すべきパラメータがないので，普通の相関図になるだけではないでしょうか？

> r<-readMoments(names=c("乳製品", "総熱量", "直腸ガン"))
1: 1.00
2: 0.65 1.00
4: 0.20 0.58 1.00
7: 
Read 6 items
> model<-specifyModel()
1: 乳製品<->総熱量, r1, NA
2: 総熱量<->直腸ガン, r2, NA
3: 直腸ガン<->乳製品, r3, NA
4: 乳製品<->乳製品,NA,1
5: 総熱量<->総熱量,NA,1
6: 直腸ガン <->直腸ガン,NA,1
7: 
Read 6 records
> ans <- sem(model, r, N=100)
> summary(ans)

 Model Chisquare =  0   Df =  3 Pr(>Chisq) = 1
 AIC =  6
 BIC =  -13.81551

 Normalized Residuals
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 

 Parameter Estimates
   Estimate Std Error  z value   Pr(>|z|)                         
r1 0.65     0.04791380 13.566028 6.368060e-42 総熱量 <--> 乳製品  
r2 0.58     0.05660074 10.247215 1.218024e-24 直腸ガン <--> 総熱量
r3 0.20     0.08632217  2.316902 2.050907e-02 乳製品 <--> 直腸ガン

 Iterations =  0 
> 
> pathDiagram(ans, "comp-model", ignore.double=FALSE, edge.labels="both")
Running  dot -Tpdf -o comp-model.pdf  comp-model.dot

No.20376　Re: 相関関係の視覚化　　【コロン】　2013/10/25(Fri) 11:34

青木先生

ありがとうございました。私のまちがいがわかりました。
実は先生と類似したことは行っておりました。

1: 乳製品<->総熱量, r1, NA
2: 総熱量<->直腸ガン, r2, NA
3: 直腸ガン<->乳製品, r3, NA
4: 乳製品<->乳製品,NA,1
5: 総熱量<->総熱量,NA,1
6: 直腸ガン <->直腸ガン,NA,1

最初の3行は同じでしたが，次の3行を,

乳製品<->乳製品, e1，NA

のようにしていました。

ありがとうございました。