* S3A3 Example done in class on October 26, 2017;

* Diagnostic plots;

libname S3A3 "D:\STAT 3A03\Fall 2017\Data";

PROC IMPORT Out=S3A3.Data1
    Datafile="D:\STAT 3A03\Fall 2017\Data\Mydata1.txt"
    DBMS=DLM REPLACE;
	Getnames=YES;
	Datarow=2;
	run;

* Diagnostic Plots for the first response variable;
Proc REG Data=S3A3.Data1 plots=none;
   Model Y1=X1;
   Plot Y1*X1;
   Plot Student.*Pred.;
   Plot Student.*nqq.;
run;

* Diagnostic Plots for the second response variable;
Proc REG Data=S3A3.Data1 plots=none;
   Model Y2=X1;
   Plot Y2*X1;
   Plot Student.*Pred.;
   Plot Student.*nqq.;
run;

* Diagnostic Plots for the third response variable;
Proc REG Data=S3A3.Data1 plots=none;
   Model Y3=X1;
   Plot Y3*X1;
   Plot Student.*Pred.;
   Plot Student.*nqq.;
run;

* The Real Estate Dataset;
PROC IMPORT Out=S3A3.House
    Datafile="C:\Users\Angelo\Documents\McMaster Teaching\STAT 3A03\Fall 2017\Data\Realestate.txt"
    DBMS=DLM REPLACE;
	Getnames=YES;
	Datarow=2;
	run;

* Scatterplot Matrix;
PROC SGScatter Data=S3A3.House;
   Matrix Price Tax Bath Bed Lot Space Age;
run;

* Full Model and diagnostic plots;
PROC REG Data=S3A3.House plot=none;
  Model Price=Tax Bath Bed Lot Space Age;
  Plot Price*Pred.;
  Plot Student.*Price;
  Plot Student.*nqq.;
  Plot Student.*Obs.;
run;

PROC IMPORT Out=S3A3.Data2
    Datafile="C:\Users\Angelo\Documents\McMaster Teaching\STAT 3A03\Fall 2017\Data\Mydata2.txt"
    DBMS=DLM REPLACE;
	Getnames=YES;
	Datarow=2;
run;

* Scatterplot matrix;	
PROC SGSCATTER Data=S3A3.Data2;
   Matrix y x1 x2 x3;
run;

PROC REG Data=S3A3.Data2 plots=none;
  Model Y=X1 X2 X3/partial;
run;
quit;