#行の中で、「#」以降は、コメントアウトで処理しません。
#分析ごとにRエディタにコピーペーストして実行してください。
#ただし、一気に処理してもエラーは出ません。
#====================================================
#図1-8の分析プログラム
#パッケージlavaanを使用(要パッケージインストール)
library(lavaan)
#パッケージpsychを使用(要パッケージインストール)
library(psych)
#NAを含むデータを除外する
bfi2<-na.omit(bfi)
#モデル指定
model.fig1.8
<- '
f
=~ N1+N2+N3 #N1, N2, N3で因子fを測定
N1
~~ N1
#N1の誤差分散を推定
N2
~~ N2
#N2の誤差分散を推定
N3
~~ N3
#N3の誤差分散を推定
f
~~ 1*f
#因子fの分散を1に固定
'
#lavaanを実行する
fit.fig1.8
<- lavaan(model.fig1.8, data=bfi2[,16:18])
#結果を出力する
summary(fit.fig1.8,
standardized=T, rsquare=T, fit.measure=T)
#====================================================
#2014.09.15 尾崎幸謙・荘島宏二郎