#表7-4の分析

library(lavaan)

#変数名

varnames<-c("t1","t2")

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

mz <- '

0.675,

0.513,0.714'

#1卵性の共分散行列の作成

mz.cov <- getCov(mz,names=varnames)

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

dz <- '

0.621,

0.434,0.623'

#2卵性の共分散行列の作成

dz.cov <- getCov(dz,names=varnames)

#各群の標本サイズの指定

N<-list(mzG=102, dzG=97)

#各群の共分散行列をまとめる

twin.cov<-list(mzG=mz.cov, dzG=dz.cov)

#フルモデル（ACEモデル）の指定（表7-4）

model.table7.4ace <- '

fA1 =~ c(a1,a1)*t1

fA2 =~ c(a1,a1)*t2

fC =~ c(c1,c1)*t1+c(c1,c1)*t2

fE1 =~ c(e1,e1)*t1

fE2 =~ c(e1,e1)*t2

fA1 ~~ 1*fA1

fA2 ~~ 1*fA2

fA1 ~~ c(1,0.5)*fA2

fC ~~ 1*fC

fE1 ~~ 1*fE1

fE2 ~~ 1*fE2

t1 ~~ 0*t1

t2 ~~ 0*t2

'

#lavaanを実行する

fit.table7.4ace <- lavaan(model.table7.4ace, sample.cov=twin.cov, sample.nobs=N)

summary(fit.table7.4ace, standardized=T, rsquare=T, fit.measures=TRUE)

#遺伝・非共有環境モデル（AEモデル）の指定（表7-4）

model.table7.4ae <- '

fA1 =~ c(a1,a1)*t1

fA2 =~ c(a1,a1)*t2

fE1 =~ c(e1,e1)*t1

fE2 =~ c(e1,e1)*t2

fA1 ~~ 1*fA1

fA2 ~~ 1*fA2

fA1 ~~ c(1,0.5)*fA2

fE1 ~~ 1*fE1

fE2 ~~ 1*fE2

t1 ~~ 0*t1

t2 ~~ 0*t2

'

#lavaanを実行する

fit.table7.4ae <- lavaan(model.table7.4ae, sample.cov=twin.cov, sample.nobs=N)

summary(fit.table7.4ae, standardized=T, rsquare=T, fit.measures=TRUE)

#共有環境・非共有環境モデル（CEモデル）の指定（表7-4）

model.table7.4ce <- '

fC =~ c(c1,c1)*t1+c(c1,c1)*t2

fE1 =~ c(e1,e1)*t1

fE2 =~ c(e1,e1)*t2

fC ~~ 1*fC

fE1 ~~ 1*fE1

fE2 ~~ 1*fE2

t1 ~~ 0*t1

t2 ~~ 0*t2

'

#lavaanを実行する

fit.table7.4ce <- lavaan(model.table7.4ce, sample.cov=twin.cov, sample.nobs=N)

summary(fit.table7.4ce, standardized=T, rsquare=T, fit.measures=TRUE)

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

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