[fpr 2413] R のpackage SEM使用経験

Masahiko AIDA

相田@ミシガン大学です。

Rのsem packageだとさらに外部のソフトと連携して、図まで自動でかけます。
SEMのソフトウェアで自動で図がかけるものってあまり無いのではないでしょう
か。(Lisrelは書けると聞きましたが)

私が個人的に思うわざわざRを使うメリットはBootstrapなんかと組み合わせた
推定も柔軟にプログラムできる点でしょうか。AMOSでもその点はできますけど、
AMOS Basic使うぐらいならMplusを使った方がいい気がしますし、GUIは手が疲れ
ますしね。ただまだRのsem packageはできることが限られているのが難点ですね。
ま、ソースも公開されているから関心があれば付け足せばいいのでしょうけど。

夏にミシガン大学のサマースクールでJohn FoxがRの授業をしています。$2100
払えばKenn Bollenの授業も受けられますが。。。(5日間)

On Wed, 21 May 2003 17:00:50 +0900
Mitsuo Igarashi <mitsu5 (at) ruby.famille.ne.jp> wrote:

> Mx でお世話になった五十嵐です。
> 
> R の package SEM を使ってみましたのでご報告させてください。
> 私みたいな初心者でもSEMに入門するためのソフトとしては結構良い感じがします。
> 
> 以下が使用例です。
> 豊田先生の 「共分散構造分析 入門編」第4章 確認的因子分析モデル趣味の印象
> を分析してみました。 入力データはcovarianceにしてから入れてあります。
> 
> ## プログラム ##
> ## これを R のコンソールにコピーすれば結果が出力されます。##
> 
> library(sem)
> 
> syumi4.mh<-matrix(c(
> 0.862,0,0,0,0,0,
> 0.489,1.089,0,0,0,0,
> -0.056,-0.189,0.606,0,0,0,
> -0.122,-0.222,0.406,0.606,0,0,
> 0.016,-0.044,-0.089,-0.022,0.262,0,
> 0.067,0.133,-0.200,-0.167,0.200,0.533
> ),ncol=6,byrow=T)
> 
> model.mh<-matrix(c(
>          'F1 -> V1', 'a1',NA,
>          'F1 -> V2', 'a1',NA,
>          'F2 -> V3', 'a2',NA,  
>          'F2 -> V4', 'a2',NA,
>          'F3 -> V5', 'a3',NA,
>          'F3 -> V6', 'a3',NA,
>          'V1 <-> V1','e1',NA,
>          'V2 <-> V2','e1',NA,
>          'V3 <-> V3','e2',NA, 
>          'V4 <-> V4','e2',NA,
>          'V5 <-> V5','e3',NA,
>          'V6 <-> V6','e3',NA,
>          'F1 <-> F1',NA,1, 
>          'F2 <-> F2',NA,1, 
>          'F3 <-> F3',NA,1, 
>          'F1 <-> F2','f12',NA, 
>          'F2 <-> F3','f23',NA, 
>          'F1 <-> F3','f13',NA  
>  ),ncol=3,byrow=T)
> 
> obs.vars.mh <- c('V1','V2','V3','V4','V5','V6')
> rownames(syumi4.mh) <- colnames(syumi4.mh) <- obs.vars.mh
> 
> sem.mh <- sem(model.mh, syumi4.mh, 30)
> 
> summary(sem.mh)
> 
> ## 以下が出力結果です。##
> > summary(sem.mh)
> 
>  Model Chisquare =  10.601   Df =  12 Pr(>Chisq) = 0.56335
>  Goodness-of-fit index =  0.90117
>  Adjusted goodness-of-fit index =  0.82705
>  RMSEA index =  0   90 % CI: (0, 0.17089)
>  BIC =  -51.714 
> 
>  Normalized Residuals
>      Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
> -1.32e+00 -4.67e-01 -1.08e-07  1.34e-07  3.31e-01  1.32e+00 
> 
>  Parameter Estimates
>     Estimate Std Error  z value   Pr(>|z|)           
> a1   0.69929  0.144910  4.82565 1.3954e-06 V1 <--- F1
> a2   0.63718  0.106321  5.99301 2.0600e-09 V3 <--- F2
> a3   0.44721  0.092411  4.83938 1.3024e-06 V5 <--- F3
> e1   0.48650  0.127811  3.80639 1.4101e-04 V1 <--> V1
> e2   0.20000  0.052575  3.80412 1.4231e-04 V3 <--> V3
> e3   0.19750  0.051917  3.80415 1.4229e-04 V5 <--> V5
> f12 -0.33047  0.235806 -1.40147 1.6107e-01 F2 <--> F1
> f23 -0.41936  0.224736 -1.86602 6.2039e-02 F3 <--> F2
> f13  0.13750  0.274794  0.50037 6.1681e-01 F3 <--> F1
> 
>  Iterations =  48 
> ## 以上出力結果、終了。##
> 
> ###   ###
> CONSTRAINTS は下のようになるように 、model.mh<-matrix(c( のところで設定してあります。
> V1<- F1 = V2<- F1,  V3<- F2 = V4<- F2,  V5<- F3 = V6<- F3
> V1-error=V2-error, V3-e=V4-e, V5-e=V6-e
> ###   ###
> 
> Free software が魅力であるとともに、結構使いやすい気がします。
> ご検討のほどよろしくお願いいたします。
> 
> まだ、sem の未熟者ですので、誤りがあったらご勘弁ください。
> --------========----------
> Mitsuo Igarashi
> mitsu5 (at) ruby.famille.ne.jp

--------
Masahiko Aida <aida (at) highway.ne.jp>
MSN: masa_aida (at) hotmail.com
http://www102.sakura.ne.jp/~spicetea/masahiko/


スレッド表示 著者別表示 日付順表示 トップページ

ここは心理学研究の基礎メーリングリストに投稿された過去の記事を掲載しているページです。