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

No.14663　Rで共分散構造分析　　【看護研究にて】　2011/05/21(Sat) 19:11

Rで，相関行列から共分散構造分析をしようとしています。
次の通りの相関行列とプログラムで実施しようとしたのですが，警告が出ます。モデルがおかしいのかと思いましたが，自由度は2で，問題ないように見えるのですが，なにかまずいところがあるでしょうか。

nc <- 47
dat <- gendat(nc, matrix(c(
1,-0.376716550308732,0.228151596610790,-0.158092112183398,-0.381075749056570,
-0.376716550308732,1,-0.657586848928479,0.520742827734458,0.429423509822718,
0.228151596610790,-0.657586848928479,1,-0.637686215418466,-0.313031396683772,
-0.158092112183398,0.520742827734458,-0.637686215418466,1,0.623882553078004,
-0.381075749056570,0.429423509822718,-0.313031396683772,0.623882553078004,1
),ncol=5))
colnames(dat) <- c("COM","RSE","STT","SPQ","PPQ")

model <- specify.model()
SPQ < QOL,path1,1
PPQ < QOL,path2,NA
QOL < RSE,path3,1
QOL < STT,path4,NA
QOL < COM,path5,NA
RSE < COM,path6,NA
STT < COM,path7,NA
QOL <> QOL,s1,NA
RSE <> RSE,s2,NA
STT <> STT,s3,NA
PPQ <> PPQ,s4,NA
SPQ <> SPQ,s5,NA
COM <> COM,NA,1
RSE <> STT,s6,NA

ans <- sem(model,cor(dat),nc)
summary(ans)

プログラムは以上です。SPQとPPQは，QOL（潜在変数）を測定する尺度得点です。
COMからQOLへの直接効果と，RSE，STTをそれぞれ介した間接効果の検討をしています。

警告は，
In sem.default(ram = ram, S = S, N = N, param.names = pars, var.names = vars, :
Could not compute QR decomposition of Hessian.
Optimization probably did not converge.
です。

何卒宜しくお願い申し上げます。

No.14664　Re: Rで共分散構造分析　　【青木繁伸】　2011/05/21(Sat) 20:51

回答ではないのですが，実際の相関係数行列があるのに gendat を使う意図がわからない。単に元データから相関係数行列を計算してそれを使えばよいし，それが不可能だとしても今回のプログラムで使用された相関係数行列を read.moments で読めば良いだけだと思いますが。

No.14665　Re: Rで共分散構造分析　　【青木繁伸】　2011/05/21(Sat) 22:04

解の一部は求まっているようですね。それが正しいかどうかはよく分かりませんが。
print(ans) や path.diagram(ans, ...) で結果の一部分は見えますが，それが，正しい結果かどうかは私には分かりません。
> print(ans)

 Model Chisquare =  -157821410   Df =  2 

      path1       path2       path3       path4       path5 
 0.57904151  0.18915127  0.55857612 -0.74364507 -0.25308604 
      path6       path7          s1          s2          s3 
-0.37671655  0.22815160  0.36148966  0.85808464  0.94794685 
         s4          s5          s6 
-0.01820046 -0.41884166 -0.57163837 

 Iterations =  26  

No.14667　Re: Rで共分散構造分析　　【看護研究にて】　2011/05/22(Sun) 11:25

そうでした。gendatについて，特に意図はありません。read.momentsですべきでした。
仰る通り，係数は計算されているのですが，適合度指標などがなく，正しいかどうか分からないのです。