options ls=80 nodate pageno=1 NOFMTERR;  

 * radon.sas  reproduces examples in Gelman and Hill, Ch 12 ;
 * CSV data saved as: arhm.radon n = 919 observations from Minnesota
   variables: county_number  radon floor uranium;
   
* both radon and uranium (ppm) are on log scale;

 *proc glm data=	arhm.radon; * p 255;
 *  model radon=floor / solution;

 proc mixed  data=	arhm.radon;	* p 259;
   class  county_number;
   model radon=floor / solution   ddfm=bw;
   random intercept / subject=county_number v vcorr solution;
   ODS output  SolutionR =a;

 proc print data=a(obs=5);

 proc mixed  data=	arhm.radon;	* p 266;
   class  county_number;
   model radon=floor uranium/ solution   ddfm=bw;
   random intercept / subject=county_number v vcorr solution;
   ODS output SolutionR = two;

 proc print data=two(obs=5);

 run; quit;