myscratchpad

# 狩野裕・三浦麻子 (2002). AMOS, EQS, CALISによるグラフィカル多変量解析――目で見る共分散構造分析　第2版　現代数学者　
 # p2, 表1のデータより
kk <- c(89,99,128,98,52,47,40,39,38,48,27,23) # 価格
sk <- c(4.3,1.9,5.2,5.1,4.0,4.8,8.7,8.2,3.3,3.9,8.2,7.2) # 走行距離
jn <- c(5,4,2,3,6,8,7,7,10,6,8,8) # 乗車年数
skn <- c(24,18,13,4,15,24,3,6,14,0,24,24) # 車検
dat <- data.frame(kk, sk, jn, skn)
dat2 <- scale(dat)

summary(lm(kk~jn))

library(lavaan) # 残念ながら、相関行列をデータとしてうけとることはできないらしい
# p15、パス解析モデル1
model1 <- '
 sk ~ jn
 kk ~ sk + skn
'

model1sc <- '
 sk ~ jn
 kk ~ sk + skn
'

# p16、パス解析モデル2
model2 <- '
 sk ~ jn
 kk ~ sk + skn + jn
'

model2sc <- '
 sk ~ jn
 kk ~ sk + skn + jn
'

fit1 <- sem(model1, data=dat)
fit1sc <- sem(model1sc, data=dat2)
fit2 <- sem(model2, data=dat)
fit2sc <- sem(model2sc, data=dat2)
summary(fit1)
summary(fit1sc)
summary(fit2)
summary(fit2sc)

参考
Rosseel, Y. lavaan: an R package for structural equation modeling and more