第１３回　構造方程式モデリング

構造方程式モデリング
SEM : Structural Equation Modeling
共分散構造分析（covariance structural analysis）

​

因子分析と重回帰分析を同時に行う
観測変数、潜在変数

確認的因子分析
　表 931
　図 932

​

観測変数

潜在変数

​

​課題　以下のモデルについて確認的因子分析をRで実行し、各因子から形容詞尺度（観測変数）への標準化係数（パス係数）を報告しなさい。

# 確認的因子分析モデル
model01 <- cfa()
 f1: AkaruiKurai,JimiHade,NibuiSurudoi,SizukaSawagasii,YoukiInnki,DoutekiSeiteki # 測定方程式（活動性）    

f2: KinntyoYurui,KaruiOmoi # 測定方程式（力量）    

f3: KitanaiUtsukusi,SukiKirai,WaruiYoi # 測定方程式（価値）

​

​

> model01
   Path                                Parameter               StartValue
1  f1 -> AkaruiKurai                   <fixed>                 1         
2  f1 -> JimiHade                      lam[JimiHade:f1]                  
3  f1 -> NibuiSurudoi                  lam[NibuiSurudoi:f1]              
4  f1 -> SizukaSawagasii               lam[SizukaSawagasii:f1]           
5  f1 -> YoukiInnki                    lam[YoukiInnki:f1]                
6  f1 -> DoutekiSeiteki                lam[DoutekiSeiteki:f1]            
7  f2 -> KinntyoYurui                  <fixed>                 1         
8  f2 -> KaruiOmoi                     lam[KaruiOmoi:f2]                 
9  f3 -> KitanaiUtsukusi               <fixed>                 1         
10 f3 -> SukiKirai                     lam[SukiKirai:f3]                 
11 f3 -> WaruiYoi                      lam[WaruiYoi:f3] 

12 f1 <-> f1                           V[f1]                             
13 f2 <-> f2                           V[f2]                             
14 f3 <-> f3                           V[f3]                             
15 f1 <-> f2                           C[f1,f2]                          
16 f1 <-> f3                           C[f1,f3]                          
17 f2 <-> f3                           C[f2,f3] 

18 AkaruiKurai <-> AkaruiKurai         V[AkaruiKurai]                    
19 JimiHade <-> JimiHade               V[JimiHade]                       
20 NibuiSurudoi <-> NibuiSurudoi       V[NibuiSurudoi]                   
21 SizukaSawagasii <-> SizukaSawagasii V[SizukaSawagasii]                
22 YoukiInnki <-> YoukiInnki           V[YoukiInnki]                     
23 DoutekiSeiteki <-> DoutekiSeiteki   V[DoutekiSeiteki]                 
24 KinntyoYurui <-> KinntyoYurui       V[KinntyoYurui]                   
25 KaruiOmoi <-> KaruiOmoi             V[KaruiOmoi]                      
26 KitanaiUtsukusi <-> KitanaiUtsukusi V[KitanaiUtsukusi]                
27 SukiKirai <-> SukiKirai             V[SukiKirai]                      
28 WaruiYoi <-> WaruiYoi               V[WaruiYoi] 

​

# 確認的因子分析モデル　モデルの当てはめ

fit01 <- sem(model=model01,S=cov(data1),N=nrow(data1))
summary(fit01,fit.indices=c("GFI","AGFI","CFI","NFI","SRMR","RMSEA","AIC"))
standardizedCoefficients(fit01)

​

#作図 # dot言語によるパス図の出力。

pathDiagram(fit01, ignore.double=FALSE, edge.labels="values", digits=2,standardize=TRUE)