共分散構造分析[R編]ってないの？[第３章6]その2

「共分散構造分析［Amos編］－構造方程式モデリング－」の63頁図3.18をRで計算します。
セミナー評価という潜在因子は誤差を考慮せずに観測変数に重みをかけて足し上げた合成指標ということです。図3.17のセミナー評価は重回帰分析のように誤差を計算させましたが、図3.18では誤差(因子独自の分散)を0に制約するということになるはずです。
モデルは次のように書き、 semi.3.18.txtとして保存しました。

テキスト -> セミナー評価, p1, NA
プレゼン -> セミナー評価, p2, NA
ペース   -> セミナー評価, p3, NA
講師対処 -> セミナー評価, p4, NA
セミナー評価 -> 満足度,   p5, NA
セミナー評価 -> 理解度,   p6, NA
セミナー評価 -> 目的一致, NA, 1　　　　　 #ここは1に制約
テキスト <-> テキスト, v1, NA
プレゼン <-> プレゼン, v2, NA
ペース <-> ペース,     v3, NA
講師対処 <-> 講師対処, v4, NA
満足度 <-> 満足度,     e1, NA
理解度 <-> 理解度,     e2, NA
目的一致 <-> 目的一致, e3, NA
セミナー評価 <-> セミナー評価, NA, 0 　　    #ここは0に制約
テキスト <-> プレゼン, c1, NA
テキスト <-> ペース,   c2, NA
テキスト <-> 講師対処, c3, NA
プレゼン <-> ペース,   c4, NA
プレゼン <-> 講師対処, c5, NA
ペース <-> 講師対処,   c6, NA

ファイルを読み込んで、実行します。

semi.mdl <- specifyModel("semi.3.18.txt")
semi.sem <- sem(semi.mdl, data=semi.data)
summary(semi.sem)
stdCoef(semi.sem)
round(1-diag(semi.sem$P[-8,-8] / semi.sem$S),2)
round(cov2cor(semi.sem$P),2)