data923
attach(data923)

#多重共線性のチェック
zx <- data923[7:8]
r <- cor(zx)
VIF <- diag(solve(r))
tolerance = 1/VIF
data.frame(tolerance, VIF)

​

​

#モデルの比較　前回参照
step(model3)

#sem パス図

library("sem")

data4 <- cbind(Y,X6,X7)

model4 <- specifyEquations()  

Y = b1*X6 + b2*X7

​

sem4 <- sem(model=model4, S=cov(data4) , N=nrow(data4), fixed.x=c("X6","X7"))

summary(sem4, fit.indices=c("GFI", "AGFI", "CFI", "NFI", "RMSEA", "SRMR","AIC"))

standardizedCoefficients(sem4)

​

#作図

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

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

 

 

​

​