#図5-1の分析

library(lavaan)

#共分散行列の下三角要素の読み込み

lower <- '

14.066,

3.696, 15.535,

6.911, 2.006, 13.468'

#共分散行列の作成

dep.cov <- getCov(lower, names=c("dep","harm","friend"))

#単回帰モデルの指定（図5-1）

model.fig5.1 <- '

dep ~ harm

dep ~~ dep

'

#lavaanを実行する

fit.fig5.1 <- lavaan(model.fig5.1, sample.cov=dep.cov, sample.nobs=518)

summary(fit.fig5.1, standardized=T, rsquare=T, fit.measures=TRUE)

#====================================================

#図5-2の分析

library(lavaan)

#共分散行列の下三角要素の読み込み

lower <- '

14.066,

3.696, 15.535,

6.911, 2.006, 13.468'

#共分散行列の作成

dep.cov <- getCov(lower, names=c("dep","harm","friend"))

#重回帰モデルの指定（図5-2）

model.fig5.2 <- '

dep ~ harm+friend

dep ~~ dep

'

#lavaanを実行する

fit.fig5.2 <- lavaan(model.fig5.2, sample.cov=dep.cov, sample.nobs=518)

summary(fit.fig5.2, standardized=T, rsquare=T, fit.measures=TRUE)

#====================================================

#図5-7の分析

library(lavaan)

#共分散行列の下三角要素の読み込み

lower <- '

14.066,

3.696, 15.535,

6.911, 2.006, 13.468'

#共分散行列の作成

dep.cov <- getCov(lower, names=c("dep","harm","friend"))

#パス解析モデルの指定（図5-7）

model.fig5.7 <- '

dep ~ harm+friend

friend ~ harm

dep ~~ dep

friend ~~ friend

'

#lavaanを実行する

fit.fig5.7 <- lavaan(model.fig5.7, sample.cov=dep.cov, sample.nobs=518)

summary(fit.fig5.7, standardized=T, rsquare=T, fit.measures=TRUE)

#====================================================

#図5-8の分析

library(lavaan)

#共分散行列の下三角要素の読み込み

lower <- '

14.066,

3.696, 15.535,

6.911, 2.006, 13.468'

#共分散行列の作成

dep.cov <- getCov(lower, names=c("dep","harm","friend"))

#パス解析モデルの指定（図5-8）

model.fig5.8 <- '

dep ~ friend

friend ~ harm

dep ~~ dep

friend ~~ friend

'

#lavaanを実行する

fit.fig5.8 <- lavaan(model.fig5.8, sample.cov=dep.cov, sample.nobs=518)

summary(fit.fig5.8, standardized=T, rsquare=T, fit.measures=TRUE)

#====================================================

#2014.09.19 尾崎幸謙・荘島宏二郎