############## 共分散構造分析( lavaan ) ############################ コード概要 モデル名 <- # 測定変数の指定( =~ を使う) 潜在変数名 1=~ 観測変数 1+ 観測変数 2+ ・・・・・ 潜在変数名 2=~ 観測変数 1+ 観測変数 2+ ・・・・・ # 回帰式( ~ を使う) 従属変数名 1~ 説明変数 1+ 説明変数 2+ ・・・ 従属変数名 2~ 説明変数 1+ 説明変数 2+ ・・・ # 分散、共分散( ~~ を使う) 変数名 i~~ 変数名 i # 分散 変数名 j~~ 変数名 k # 共分散 モデル部分はシングルカンマ( ' )でくくる すべての変数の分散を推定する場合は分散の指定は省力できる パラメタ値の推定 lavaan オブジェクト名 <-sem (モデル名, data= データ名) または lavaan オブジェクト名 <-lavaan (モデル名, data= データ名, model.type="sem" , auto.var=TRUE , auto.fix.first=TRUE , auto.cov.lv.x=FALSE , auto.cov.y=FALSE ) #auto.var=TRUE : すべての変数の分散または残差分散を推定する #auto.fix.first=TRUE: 非標準化解において測定変数のうち最初の 1 つのパス係数を 1 に固定する #auto.cov.lv.x=FALSE: 潜在変数間の共分散を自動的には推定しない #auto.cov.y=FALSE : 従属変数の残差間の共分散を自動的には推定しない #int.ov.free : 観測変数の切片を自由にするか( FALSE なら 0 に固定) #int.lv.free : 潜在変数の切片を自由にするか( FALSE なら 0 に固定) summary ( lavaan オブジェクト名, fit.measures=TRUE , standardized=TRUE ) # 出力オプシ...
コメント
コメントを投稿