options  nocenter  ps = 65 ls = 75;

FILENAME navn URL "http://staff.pubhealth.ku.dk/~lts/basal/data/rmr.txt";

data rmrdata;
infile navn firstobs=2;
input bw rmr;

/* de to næste linier er til senere brug */
bw70=bw-70;
if rmr>2000 then delete;

label bw='body weight'
      rmr='resting metabolic rate';
run;

title 'SPM 1';

proc sgplot data=rmrdata;
scatter Y=rmr X=bw;
run;

title 'SPM 2';

proc sgplot data=rmrdata;
reg Y=rmr X=bw;
run;

proc sgplot data=rmrdata;
loess Y=rmr X=bw;
run;

title 'SPM 2+5+6';

ods graphics on;
proc glm plots=(diagnosticspanel residuals(smooth)) data=rmrdata; 
  model rmr=bw / solution clparm;
run;

title 'SPM 3';

proc glm data=rmrdata; 
  model rmr=bw / solution clparm;
estimate 'rmr ved 70 kg' intercept 1 bw 70;
run;

title 'SPM 4';

proc glm data=rmrdata; 
  model rmr=bw / solution clparm;
estimate 'rmr ved 80 kg' intercept 1 bw 80;
run;

/**********************************************************/



FILENAME navn URL 
    "http://staff.pubhealth.ku.dk/~lts/basal/data/sundby_lille.txt";

title 'SPM 2';

data sundby;
infile navn firstobs=2;
input kon v75 v76 v17 v24af;

if kon<0 then delete;

if kon=2 then gender='female';
if kon=1 then gender='male';

vaegt=v75;
hoejde=v76/100;

log10vaegt=log10(vaegt);
log10hoejde=log10(hoejde);

bmi=vaegt/hoejde**2;
run;

proc means data=sundby;
run;

title 'SPM 3';

proc sort data=sundby; by gender;
run;

proc sgplot data=sundby;
reg Y=log10vaegt X=log10hoejde / group=gender;
run;

proc sgpanel data=sundby;
panelby gender;
reg Y=log10vaegt X=log10hoejde;
run;

proc sgpanel data=sundby;
panelby gender;
reg Y=bmi X=hoejde;
run;

title 'SPM 3';

proc glm plots=(diagnosticspanel residuals(smooth)) data=sundby; by gender;
      model log10vaegt=log10hoejde / noint solution;
run;

proc corr spearman data=sundby; by gender;
    var bmi hoejde;
run;