Analysis of regression of alcdep on conduct (broad), majdep, panatt ! estimating probit regression coefficients ! ! /data2/boulder/probitmult_withpairs.mx ! First group is MZF complete pairs data NI=9 NG=3 la cond majdep panatt alcdep xcond xmjdep xpnatt xlcdep wt ordinal fi=femmzf.dat definition_variables cond majdep panatt xcond xmjdep xpnatt wt / Begin matrices; O FU 1 6 ! store definition variable (i.e. vara) here K FU 6 2 fr ! regression coefficient T FU 1 1 fr ! threshold value V FU 1 1 ! variance (=unity) R FU 1 1 fr ! tetrachoric correlation (to be estimated) W FU 1 1 ! weight variable end matrices; SP O -1 -2 -3 -4 -5 -6 SP W -7 MAT V 1.00 MAT R 0.60 SP T 50 MAT T 0.1 SP K 100 0 101 0 102 0 0 100 0 101 0 102 MAT K 0.05 0 0.05 0 0.05 0 0 0.05 0 0.05 0 0.05 FREQ W; TH -(T|T)-O*K; ! Note that we now have thresholds for both twins CO V|R_ R|V; bo -5.0 5.0 t(1,1) k(1,1) K(2,1) K(3,1) bo -0.999 0.999 r(1,1) interval k(1,1) k(2,1) k(3,1) t(1,1) r(1,1) OPT func=1.E-12 OPT RS END Analysis of regression of alcdep on conduct (broad), majdep, panatt ! /data2/boulder/probitmult_withpairs.mx ! Second group is DZF complete pairs data NI=9 la cond majdep panatt alcdep xcond xmjdep xpnatt xlcdep wt ordinal fi=femdzf.dat definition_variables cond majdep panatt xcond xmjdep xpnatt wt / Begin matrices; O FU 1 6 ! store definition variable (i.e. vara) here K FU 6 2 fi ! regression coefficient T FU 1 1 fi ! threshold value V FU 1 1 ! variance (=unity) R FU 1 1 fr ! tetrachoric correlation (to be estimated) W FU 1 1 ! weight variable end matrices; SP O -1 -2 -3 -4 -5 -6 SP W -7 MAT V 1.00 MAT R 0.30 SP T 50 SP K 100 0 101 0 102 0 0 100 0 101 0 102 FREQ W; TH -(T|T)-O*K; ! Note that we now have thresholds for both twins CO V|R_ R|V; bo -0.999 0.999 r(1,1) interval r(1,1) OPT RS END Analysis of regression of alcdep on conduct (broad), majdep, panatt ! This group is for singleton women data NI=5 la cond majdep panatt alcdep wt ordinal fi=femsing.dat definition_variables cond majdep panatt wt / Begin matrices; O FU 1 3 ! store definition variable (i.e. vara) here K FU 3 1 fr ! regression coefficient T FU 1 1 fr ! threshold value V FU 1 1 ! variance (=unity) W FU 1 1 ! weight variable end matrices; SP O -1 -2 -3 SP W -4 MAT V 1.00 SP T 50 MAT T 0.1 SP K 100 101 102 MAT K 0.05 0.05 0.05 FREQ W; TH -T-O*K; CO V; OPT func=1.E-12 OPT RS END