The data set consists of 12666 observations at level 1. We asked participants to report their behaviour during social interactions (multiple/day) for a period of 2 weeks. Specifically, we looked at the effect of drinking alcohol on social behaviour.
In this example here, the effect of CALC (event-level drinking) on quarrelsome behaviour and how this relationship is influenced by Cheavy (between person-level drinking variable).
proc mixed data=in.finalALC METHOD=ML NOCLPRINT COVTEST ;
class id study;
model IMQURBEH = calc cheavy
calc*cheavy
/SOLUTION
DDFM=BETWITHIN htypE=1,3 ;
RANDOM intercept cdrinkyn/type=un(1) sub=id ;
repeated /type=ar(1) sub=id group=study;
run;
The full code I used for the proc template is the following:
proc template;
define table Stat.Mixed.FTests;
notes "F-tests table";
column Effect Label Var NumVar Margins ByLevel AtVar DF NumDF DenDF ChiSq
FValue ProbChiSq OrdChiSq ProbF OrdF ANOVANumDF ANOVADenDF
ANOVAFValue ANOVAProbF ANOVAProbFinf;
header anovafheader;
define Effect;
header = "Effect";
style = RowHeader;
id;
end;
define Label;
header = "Label";
id;
end;
define Var;
id;
generic;
end;
define NumVar;
translate _val_=._ into "";
just = l;
id;
generic;
end;
define Margins;
header = "Margins";
end;
define ByLevel;
header = "By Level";
end;
define AtVar;
format = 6.2;
generic;
end;
define NumDF;
header = "Num DF";
format = 10.0;
end;
define DF;
header = "DF";
format = 10.0;
end;
define DenDF;
header = "Den DF";
format = best8.;
end;
define ChiSq;
header = ";Chi-Square";
translate _val_=.I into " Infty";
format = 7.2;
end;
define FValue;
parent = Common.ANOVA.FValue;
end;
define ProbChiSq;
width = 12;
parent = Common.ProbChiSq;
end;
define OrdChiSq;
header = "Ord ChiSq";
format = 6.4;
end;
define ProbF;
parent = Common.ANOVA.ProbF;
end;
define OrdF;
header = "Ord F";
format = 6.4;
end;
define anovafheader;
text "ANOVA F";
space = 0;
expand = "-";
end = ANOVAProbFInf;
start = ANOVANumDF;
spill_margin;
end;
define ANOVANumDF;
header = "Num DF";
format = best5.;
glue = 2;
end;
define ANOVADenDF;
header = "Den DF";
format = best5.;
glue = 2;
end;
define ANOVAFValue;
header = "Value";
glue = 2;
parent = Common.ANOVA.FValue;
end;
define ANOVAProbF;
header = "Pr > F(DDF)";
glue = 2;
parent = Common.ANOVA.ProbF;
end;
define ANOVAProbFInf;
header = "Pr > F(infty)";
glue = 2;
parent = Common.ANOVA.ProbF;
end;
use_format_defaults;
end;
run;
Hope that makes things clearer / was what you were asking for.
