TITLE:
Power Simulation based on reading gxe results;
MONTECARLO:
NAMES = mod x1 x2;
genclasses = cg(2);
classes = cg(2);
NOBSERVATIONS = 100;
NREPS = 1000;
SEED = 293;
REPSAVE = ALL;
SAVE = sim*.dat;
PATMISS =
x1(0) x2(0)|
x1(1) x2(0)|
x1(0) x2(1);
PATPROBS = .94 | .03 | .03;
ANALYSIS:
TYPE = MIXTURE RANDOM;
ALGORITHM IS INTEGRATION;
MODEL IS NOCOVARIANCES;
MODEL POPULATION:
%OVERALL%
[x1-x2@-.151];
a1 BY x1@.327;
a1 BY x2@.327;
a2 BY x1@0; !a2 doesn't exist for MZs
!a2 BY x2@0; !a2 doesn't exist for MZs
c BY x1@.678;
c BY x2@.678;
e1 BY x1@-.422;
e2 BY x2@-.422;
x1@.001; x2@.001;
a1@1; a2@1; c@1; e1@1; e2@1;
[a1@0 a2@0 c@0 e1@0 e2@0 mod@0];
mod@1; !.989
x1 on mod@.480;
x2 on mod@.480;
a1 WITH a2@0; !value for MZ twins
e1 WITH e2@0;
ax1 | a1 xwith mod;
ax2 | a2 xwith mod;
cx | c xwith mod;
ex1 | e1 xwith mod;
ex2 | e2 xwith mod;
x1 on ax1@.05; !vary this
x2 on ax1@.05; !vary this
x2 on ax2@0;
x1 on cx@-.015;
x2 on cx@-.015;
x1 on ex1@.034;
x2 on ex2@.034;
[cg#1*-.944]; !log odds of being in MZ group (ln(.28/.72))
%cg#1%
%cg#2%
x2 on ax2@.05; !vary this
x2 on a1@0;
x2 on a2@.327;
a1 WITH a2@.5; !Value for DZ twins
x2 on ax1@0;