options nodate nonumber ps=55 ls=90;

title 'Random Coefficients Regression';

data one;
infile 'H:\public_html\consult\rcreg1.dat';
input county $ 1-33 region $ pop1-pop7;
* x=ranuni(1234);
run;

data two;
set one;
do decade=0 to 6;
if decade=0 then do; pop=pop1; output; end;
if decade=1 then do; pop=pop2; output; end;
if decade=2 then do; pop=pop3; output; end;
if decade=3 then do; pop=pop4; output; end;
if decade=4 then do; pop=pop5; output; end;
if decade=5 then do; pop=pop6; output; end;
if decade=6 then do; pop=pop7; output; end;
end;
drop pop1-pop7;

data three; set two;
pop=log(pop);
deccls=decade;
run;




proc sort;
by county decade;
run;

* proc print;

proc mixed cl;
class region county deccls;
model pop= region decade region*decade / solution;
random intercept decade / type=un subject=county solution;
repeated deccls / type=simple subject=county r rcorr;
run;


quit;