!
! MZ group: Smoking and depression from Kendler et al 1993
!
Gp 1
da ng=2 ni=4 no=415
la S1 D1 S2 D2
km
1.0
0.24 1.00
0.81 0.24 1.00
0.24 0.49 0.24 1.00
matrices
i id 2 2
b fu 2 2
g fu 4 16
p sy 16 16 
cov (i@(i-b)~)*(g*p*g')*((i@(i-b)~)')/
sp b
0 8
0 0
sp g 
2 3 0 4 0 0 0 0 0 0 0 0 0 0 0 0  
0 0 0 0 5 6 0 7 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 2 3 0 4 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 0 0 5 6 0 7
la c g
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
la c p
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
la r p
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
ma p
1
0 1
0 0 1
0 0 0 1
0 0 0 0 1
0 0 0 0 0 1
0 0 0 0 0 0 1
0 0 0 0 0 0 0 1
1 0 0 0 0 0 0 0 1
0 1 0 0 0 0 0 0 0 1
0 0 1 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 1
0 0 0 0 1 0 0 0 0 0 0 0 1
0 0 0 0 0 1 0 0 0 0 0 0 0 1
0 0 0 0 0 0 1 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
start 0.6 all
ou rs
!
! DZ group: Smoking and depression
!
Gp 2
da ni=4 no=310
la S1 D1 S2 D2
km
1.0
0.45 1.00
0.54 0.30 1.00
0.30 0.23 0.45 1.00
matrices
i id 2 2
b fu 2 2  =b(1)
g fu 4 16 =g(1)
p sy 16 16 
cov (i@(i-b)~)*(g*p*g')*((i@(i-b)~)')/
la c g
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
la c p
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
la r p
a11 c11 d11 e11 a21 c21 d21 e21 a12 c12 d12 e12 a22 c22 d22 e22
ma p
1
0 1
0 0 1
0 0 0 1
0 0 0 0 1
0 0 0 0 0 1
0 0 0 0 0 0 1
0 0 0 0 0 0 0 1
0.5 0 0 0 0 0 0 0 1
0 1 0 0 0 0 0 0 0 1
0 0 0.25 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 1
0 0 0 0 0.5 0 0 0 0 0 0 0 1
0 0 0 0 0 1 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0.25 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
ou rs df=-6 nd=8