/*Question 1*/
TITLE "HOMEWORK 5";
OPTIONS PAGENO=1;

libname labdata "c:\temp\labdata";
DATA labdata.AFIFI;
    INFILE "c:\temp\labdata\AFIFI.DAT";
    INPUT
    #1 IDNUM 	1-4 
       AGE 		5-8 
       HEIGHT 	9-12 
       SEX 		13-15 
       SURVIVE 	16 
       SHOKTYPE 17-20 
       SBP1 	21-24
       MAP1 	25-28
       HEART1 	29-32 
       DBP1 	33-36 
       CVP1 	37-40 1
       BSA1 	41-44 2 
       CARDIAC1 45-48 2 
	   APPTIME1 49-52 1
	   CIRC1 	53-56 1
       URINE1 	57-60 
	   PLASMA1 	61-64 1
	   REDCELL1 65-68 1
       HGB1 	69-72 1
	   HCT1		73-76 1
    #2 SBP2 	21-24
       MAP2 	25-28
       HEART2 	29-32 
       DBP2 	33-36 
       CVP2 	37-40 1
       BSA2 	41-44 2 
       CARDIAC2 45-48 2 
	   APPTIME2 49-52 1
	   CIRC2 	53-56 1
       URINE2 	57-60 
	   PLASMA2 	61-64 1
	   REDCELL2 65-68 1
       HGB2 	69-72 1
	   HCT2		73-76 1;

    IF SHOKTYPE=2 THEN INSHOCK=0;
    IF SHOKTYPE IN (3,4,5,6,7) THEN INSHOCK=1;
	
	shockdum2 = (shoktype=2);
	shockdum3 = (shoktype=3);
	shockdum4 = (shoktype=4);
	shockdum5 = (shoktype=5);
	shockdum6 = (shoktype=6);
	shockdum7 = (shoktype=7);

RUN;

title2 "Question 1";
proc means data=labdata.afifi;
run;
proc freq data=labdata.afifi;
  tables sex shoktype survive;
run;
proc contents data=labdata.afifi;
run;

/*Question 2*/
proc corr data=labdata.afifi;
 var map1 map2;
 title "Question 2";
run;

goptions reset=all;
goptions device=win target=winprtm;
symbol1 color=black value=dot interpol=rl line=1;
symbol2 color=black value=triangle interpol=rl line=2;
proc gplot data=labdata.afifi;
 plot map2 * map1;
 run;

 /*Question 3*/
proc reg data=labdata.afifi;
  model map2 = map1;
  plot residual.*predicted.;
  output out=regdat1 p=predict1 r=resid1;
  title "Question 3";
run; quit;

proc univariate plot normal data=regdat1;
  var resid1;
  histogram;
  qqplot / normal (mu=est sigma=est);
run;

/*Question 4: SAS/INSIGHT*/

/*Question 5: Correlations and Scatter Plot Matrix*/
title "Question 5";
proc corr;
  var map2 map1 sbp1 dbp1 cvp1 hct1 hgb1;
run;

/*Question 6: Multiple Regression*/
title "Question 6";
proc reg data=labdata.afifi;
  model map2 = map1 sbp1 dbp1 cvp1 hct1 hgb1 / tol vif collin ;
run;quit;

/*Question 7: New Regression Model*/
title "Question 7";
proc reg data=labdata.afifi;
  model map2 = map1 cvp1 hgb1 / tol vif collin;
  output out=regdat2 p=predict r=resid;
  plot residual.*predicted.;
run;quit;
proc univariate data=regdat2 plot normal;
  var resid;
  histogram;
  qqplot / normal(mu=est sigma=est);
run;