/*specifying library*/
%let direc=G:\SOE Charles\programs\;
libname here "&direc" ;
libname there "&direc.data8" ;
option spool;

data this;
set here.finaldata;
*if gender2="FEMALE";
length COL $ 9;
if closureX =0 or closureX=0.5 or closureX=1 then do;
COL="loner";
end;
if closureX =1.5 or closureX=2 or closureX=2.5 then do;
COL="normal";
end;
if closureX =3 or closureX=3.5 or closureX=4 then do;
COL="collegial";
end;

if closureX=2.5 or closureX =3 or closureX=3.5 or closureX=4 then do;
COL2="2collegial";
end;

if closureX =0 or closureX=0.5 or closureX=1 or closureX =1.5 or closureX=2 then do;
COL2="1loner";
end;

length z $ 25;
length comparison $ 25;
z=compress(TS_TALKA||COL2);
if z="3StalkzeroT1loner" then comparison="(1)Loner_Zero";
*if z="3StalkzeroT1loner" then comparison="(1)Zero";
if z="1StalkoneT1loner" then comparison="(2)Loner_One";
if z="2StalktwoTs1loner" then comparison="(3)Loner_Two";

if z="3StalkzeroT2collegial" then comparison="(4)Col_Zero";
*if z="3StalkzeroT2collegial" then comparison="(1)Zero";
if z="1StalkoneT2collegial" then comparison="(5)Col_One";
if z="2StalktwoTs2collegial" then comparison="(6)Col_Two";

*if comparison ne "(1)Zero";
if gender2="MALE";
run;



/*Interaction*/
%INCLUDE "&direc.glimmix.txt" /nosource;
%GLIMMIX(DATA=this,procopt=covtest NOCLPRINT, weight=F1weight,
STMTS=%STR(
class gender2 COL2
race2 urban gender 
F1getalong
TS_TALKA F1schooltype F1sch_ID F1famact F1palvalue
comparison;
MODEL F1college=
/*comparison*/
TS_TALKA
closure
TS_talkA*closure
/*TS_TALKA*COL2*/
BYcollege
/*gender2*/
race2 
f1ses
F1compo_score 
/*BF1devRlog*/
F1devR
urban
F1G10enrol  

F1schooltype  
F1famact 
F1palvalue F1clubtime parentschool
F1getalong F1Tutor
/solution ddfm=bw;
parms (.3) (1) /eqcons=2;
RANDOM INTERCEPT  /SUB=F1sch_ID;
),
ERROR=binomial,LINK=logit);
RUN;


%let soln=_soln;
%include "C:\temp\varcut.txt";