#############################

# 多母集団の同時分析（R: lavaanパッケージ）

# Input

#############################

 

# データの読み込み

dat <- read.csv("chapter8.csv", fileEncoding = "cp932")

 

# lavaanパッケージの読み込み

library(lavaan)

 

## モデル 1 :群間で制約を課さないモデル

model01 <-'

              情動知能の成長感 ~ 自己活用接近対処 + 自己活用回避対処 + 他者活用対処 + 現浪 + 性別 + 学年

              自己活用接近対処 ~ 挑戦 + 回避 + 現浪 + 性別 + 学年

              自己活用回避対処 ~ 挑戦 + 回避 + 現浪 + 性別 + 学年

              他者活用対処 ~ 挑戦 + 回避 + 現浪 + 性別 + 学年

              自己活用接近対処 ~~ 自己活用回避対処 + 他者活用対処

              自己活用回避対処 ~~ 他者活用対処

              挑戦 ~~ 回避

              現浪 ~~ 0*性別 + 0*学年

              性別 ~~ 0*学年

              情動知能の成長感 ~~ 情動知能の成長感

              自己活用接近対処 ~~ 自己活用接近対処

              自己活用回避対処 ~~ 自己活用回避対処

              他者活用対処 ~~ 他者活用対処

              '

              fit <- sem(model01, data= dat, group="入試形式", fixed.x = FALSE)

              summary(fit, fit.measures = TRUE)

             

 

## モデル 2:すべてのパス係数に等値制約を課すモデル

model02 <-'

              情動知能の成長感 ~ c(p1,p1)*自己活用接近対処 + c(p2,p2)*自己活用回避対処 + c(p3,p3)*他者活用対処 + c(p4,p4)*現浪 + c(p5,p5)*性別 + c(p6,p6)*学年

              自己活用接近対処 ~ c(p7,p7)*挑戦 + c(p8,p8)*回避 + c(p9,p9)*現浪 + c(p10,p10)*性別 + c(p11,p11)*学年

              自己活用回避対処 ~ c(p12,p12)*挑戦 + c(p13,p13)*回避 + c(p14,p14)*現浪 + c(p15,p15)*性別 + c(p16,p16)*学年

              他者活用対処 ~ c(p17,p17)*挑戦 + c(p18,p18)*回避 + c(p19,p19)*現浪 + c(p20,p20)*性別 + c(p21,p21)*学年

              自己活用接近対処 ~~ 自己活用回避対処 + 他者活用対処

              自己活用回避対処 ~~ 他者活用対処

              挑戦 ~~ 回避

              現浪 ~~ 0*性別 + 0*学年

              性別 ~~ 0*学年

              情動知能の成長感 ~~ 情動知能の成長感

              自己活用接近対処 ~~ 自己活用接近対処

              自己活用回避対処 ~~ 自己活用回避対処

              他者活用対処 ~~ 他者活用対処

              '

              fit <- sem(model02, data= dat, group="入試形式", fixed.x = FALSE)

              summary(fit, fit.measures = TRUE)

 

 

## モデル 3 :すべてのパス係数と誤差分散に等値制約を課すモデル

model03 <-'

              情動知能の成長感 ~ c(p1,p1)*自己活用接近対処 + c(p2,p2)*自己活用回避対処 + c(p3,p3)*他者活用対処 + c(p4,p4)*現浪 + c(p5,p5)*性別 + c(p6,p6)*学年

              自己活用接近対処 ~ c(p7,p7)*挑戦 + c(p8,p8)*回避 + c(p9,p9)*現浪 + c(p10,p10)*性別 + c(p11,p11)*学年

              自己活用回避対処 ~ c(p12,p12)*挑戦 + c(p13,p13)*回避 + c(p14,p14)*現浪 + c(p15,p15)*性別 + c(p16,p16)*学年

              他者活用対処 ~ c(p17,p17)*挑戦 + c(p18,p18)*回避 + c(p19,p19)*現浪 + c(p20,p20)*性別 + c(p21,p21)*学年

              自己活用接近対処 ~~ 自己活用回避対処 + 他者活用対処

              自己活用回避対処 ~~ 他者活用対処

              挑戦 ~~ 回避

              現浪 ~~ 0*性別 + 0*学年

              性別 ~~ 0*学年

              情動知能の成長感 ~~ c(e1,e1)*情動知能の成長感

              自己活用接近対処 ~~ c(e2,e2)*自己活用接近対処

              自己活用回避対処 ~~ c(e3,e3)*自己活用回避対処

              他者活用対処 ~~ c(e4,e4)*他者活用対処

              '

              fit <- sem(model03, data= dat, group="入試形式", fixed.x = FALSE)

              summary(fit, fit.measures = TRUE)

 

 

## モデル 4 :すべてのパス係数,誤差分散,共分散に等値制約を課すモデル

model04 <-'

              情動知能の成長感 ~ c(p1,p1)*自己活用接近対処 + c(p2,p2)*自己活用回避対処 + c(p3,p3)*他者活用対処 + c(p4,p4)*現浪 + c(p5,p5)*性別 + c(p6,p6)*学年

              自己活用接近対処 ~ c(p7,p7)*挑戦 + c(p8,p8)*回避 + c(p9,p9)*現浪 + c(p10,p10)*性別 + c(p11,p11)*学年

              自己活用回避対処 ~ c(p12,p12)*挑戦 + c(p13,p13)*回避 + c(p14,p14)*現浪 + c(p15,p15)*性別 + c(p16,p16)*学年

              他者活用対処 ~ c(p17,p17)*挑戦 + c(p18,p18)*回避 + c(p19,p19)*現浪 + c(p20,p20)*性別 + c(p21,p21)*学年

              自己活用接近対処 ~~ c(r1,r1)*自己活用回避対処 + c(r2,r2)*他者活用対処

              自己活用回避対処 ~~ c(r3,r3)*他者活用対処

              挑戦 ~~ c(r4,r4)*回避

              現浪 ~~ 0*性別 + 0*学年

              性別 ~~ 0*学年

              情動知能の成長感 ~~ c(e1,e1)*情動知能の成長感

              自己活用接近対処 ~~ c(e2,e2)*自己活用接近対処

              自己活用回避対処 ~~ c(e3,e3)*自己活用回避対処

              他者活用対処 ~~ c(e4,e4)*他者活用対処

              '

              fit <- sem(model04, data= dat, group="入試形式", fixed.x = FALSE)

              summary(fit, fit.measures = TRUE)