MONTECARLO:

names are y1-y6;
nobservations = 375;
nreps = 1000;
seed = 56798;

!The patmiss command lists each of the missing data patterns present in the;
!population, where 1's denote missing data.;
patmiss = 
y1(0) y2(0) y3(0) y4(0) y5(1) y6(1)|
y1(0) y2(0) y3(0) y4(1) y5(0) y6(1)|
y1(0) y2(0) y3(1) y4(0) y5(0) y6(1)|
y1(0) y2(1) y3(0) y4(0) y5(0) y6(1)|
y1(1) y2(0) y3(0) y4(0) y5(0) y6(1)|
y1(0) y2(0) y3(0) y4(1) y5(1) y6(0)|
y1(0) y2(0) y3(1) y4(0) y5(1) y6(0)|
y1(0) y2(1) y3(0) y4(0) y5(1) y6(0)|
y1(1) y2(0) y3(0) y4(0) y5(1) y6(0)|
y1(0) y2(0) y3(1) y4(1) y5(0) y6(0)|
y1(0) y2(1) y3(0) y4(1) y5(0) y6(0)|
y1(1) y2(0) y3(0) y4(1) y5(0) y6(0)|
y1(0) y2(1) y3(1) y4(0) y5(0) y6(0)|
y1(1) y2(0) y3(1) y4(0) y5(0) y6(0)|
y1(1) y2(1) y3(0) y4(0) y5(0) y6(0);

!The patprobs command assigns a percentage of the population to each missing data;
!pattern. In this case, an equal percentage of the population is assigned to each;
!of the 15 missing data patterns.;
patprobs = .066666|.066666|.066666|.066666|.066666|.066666|.066666|.066666|.066666|
.066666|.066666|.066666|.066666|.066666|.066666;

MODEL POPULATION:

icept linear | y1@0 y2@1 y3@2 y4@3 y5@4 y6@5;
[icept*20.67 linear*-0.66];
icept*191.55 linear*3.30;
icept with linear*-7.92;
y1-y6*169;

ANALYSIS:

type = missing;
coverage = .0000001;

MODEL: 

icept linear| y1@0 y2@1 y3@2 y4@3 y5@4 y6@5;
y1-y6;