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

No.20905　Ｒのパス解析でのエラーについて　　【ゆう】　2014/02/23(Sun) 00:58

Rを勉強中のものです。

現在，Rでパス解析を行っているのですが，よく分からないエラーが出て困っております。
下記の通りエディタに入力し，実行したところ，相関行列，共分散行列，「Read 4 records」まではちゃんと出力されるのですが，その後に見慣れないエラーが出てしまいます。
いろいろなサイト等拝見したのですが，なぜこのようなエラーが出るのか分かりませんでした。別のモデルで試してみても，相関行列を直接入力しても同じだったので，何か根本的なミスをしているように思います。

恐らく初歩的な質問かとは思いますが，ご教授いただけますと幸いです。

-----------
dat<-read.csv("Book1.csv")
r<-cor(dat)
r
cov<-cov(dat)
cov

library(sem)

model<-specify.model()
F1 -> F2, p1, NA
F2 -> F3, p2, NA
F2 <-> F2, e2, NA
F3 <-> F3, e3, NA

result<-sem(model, r, N=10)
summary(result)
-----------

実行すると，「Read 4 records」の後，以下のエラー

-----------
> result<-sem(model, r, N=10)
以下にエラー csem(model = model.description, start, opt.flag = 1, typsize = typsize, :
The matrix is non-invertable.
> summary(result)
以下にエラー summary(result) : オブジェクト 'result' がありません
-----------

No.20906　Re: Rのパス解析でのエラーについて　　【青木繁伸】　2014/02/23(Sun) 09:35

model<-specify.model()
F1 -> F2, p1
F2 -> F3, p2
F2 <-> F2, e2, NA
F3 <-> F3, e3, NA
F1 <-> F1, e1, NA

F1 <-> F1, e1, NA が必要なんじゃないですか？？
F1 <-> F1, NA, 1 の方がいいのかもしれない

No.20908　Re: Rのパス解析でのエラーについて　　【ゆう】　2014/02/23(Sun) 23:48

青木先生

ご回答をいただき，ありがとうございます。

ご指摘いただいたとおり入力すると，正しく出力されるようになりました。
モデルの概念図を見て，F1 <-> F1はいらないものと思っていました。
勉強不足を痛感しました。

日曜にもかかわらずご返信をいただき，誠にありがとうございました。