data dyadhlm; input dyad person future sex contrib culture; x1=1-(sex+1)/2; x2=-x1+1; cards; 1 1 75 -1 -10 1 1 2 90 1 -5 1 2 1 55 -1 0 1 2 2 75 1 10 1 3 1 45 -1 -10 1 3 2 33 1 -15 1 4 1 70 -1 5 1 4 2 75 1 15 1 5 1 50 -1 0 1 5 2 40 1 -5 1 6 1 85 -1 -10 -1 6 2 90 1 20 -1 7 1 75 -1 -5 -1 7 2 80 1 0 -1 8 1 90 -1 5 -1 8 2 68 1 0 -1 9 1 65 -1 0 -1 9 2 78 1 15 -1 10 1 88 -1 -15 -1 10 2 95 1 5 -1 ; *Model on page 90; PROC MIXED COVTEST; CLASS dyad; MODEL future = contrib culture contrib*culture /Solution DDFM=SATTERTH; RANDOM intercept / subject = dyad; run; *Model on page 90 that allows for negative nonindependence; proc mixed COVTEST; class dyad; MODEL future = contrib culture contrib*culture /Solution DDFM=SATTERTH; REPEATED / TYPE=CS SUBJECT=dyad; run; *Empty model on page 93 that allows for negative nonindependence; proc mixed COVTEST; class dyad; MODEL future = /Solution DDFM=SATTERTH; REPEATED / TYPE=CS SUBJECT=dyad; run; *Model that allows for heterogeneous variances; proc mixed covtest ; class dyad sex; MODEL future = contrib culture contrib*culture /Solution DDFM=SATTERTH; REPEATED sex/ TYPE=csh SUBJECT=dyad; run; *Model for a distinguishable variable; proc mixed covtest; class dyad; MODEL future = contrib culture contrib*culture /Solution DDFM=SATTERTH; REPEATED / TYPE=cs SUBJECT=dyad; run;quit; *A better model for a distinguishable variable (that includes sex); proc mixed covtest; class dyad; MODEL future = sex contrib culture contrib*culture sex*contrib /Solution DDFM=SATTERTH; REPEATED / TYPE=CSH SUBJECT=dyad; run; * "Two intercept" model; proc mixed covtest; class dyad sex; MODEL future = sex sex*contrib sex*culture sex*contrib*culture /Solution DDFM=SATTERTH noint; REPEATED / TYPE=CSH SUBJECT=dyad; run;quit; * Actual two intercept model with "zero" error variance; proc mixed; class dyad; MODEL future = x1 x2 contrib culture contrib*culture /Solution DDFM=SATTERTH noint; random x1 x2 / subject = dyad type=un; parms 2,1,2,0.0001/hold=4; run;quit;