/* LAB 4B:   */

*The Data;
filename amesh url 'http://www.openintro.org/stat/data/ames_sas.csv';

proc import datafile=amesh out=work.ames dbms=csv replace;
            getnames=yes;
            guessingrows=max;
run;

proc surveyselect data=work.ames out=work.amessample sampsize=60
     method=srs ranuni;
run;

data work.amessample;
     set work.amessample;
     keep Gr_Liv_Area;
run;

*Confidence Intervals;
proc means data=work.amessample mean;
     var Gr_Liv_Area;
run;

proc means data=work.amessample mean clm alpha=0.05;
     var Gr_Liv_Area;
run;

*Confidence Levels;
proc means data=work.ames mean;
     var Gr_Liv_Area;
run;

proc sql;
     select AVG(Gr_Liv_Area) into :popmean FROM work.ames;
run;

proc surveyselect data=work.ames out=work.amessampler sampsize=60
     method=srs reps=50 ranuni;
run;

proc means data=work.amessampler mean clm alpha=0.05 noprint;
     by replicate;
     var Gr_Liv_Area;
     output out=work.reprun mean=s_mean lclm=s_lower uclm=s_upper;
run;

proc print data=work.reprun;
  var s_mean s_lower s_upper;
run;

*On Your Own;
data work.reprun;
	set work.reprun;
	captured = (s_lower le &popmean le s_upper);
run;

proc means data=work.reprun mean;
     var captured;
run;