10.07.2015 Views

Homework 2 (Attendance 3) for Statistics 512 Applied Regression ...

Homework 2 (Attendance 3) for Statistics 512 Applied Regression ...

Homework 2 (Attendance 3) for Statistics 512 Applied Regression ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

DATA CONCVTIME3;SET CONCVTIME2;ARRAY OBS1[15] OBS101-OBS115;ARRAY OBS2[15] OBS201-OBS215;ARRAY T1[15] T101-T115;ARRAY T2[15] T201-T215;ARRAY T3[15] T301-T315;ARRAY T4[15] T401-T415;ARRAY T5[15] T501-T515;DO I = 1 TO 15;CONCT1 = T1[I];CONCT2 = T2[I];CONCT3 = T3[I];CONCT4 = T4[I];CONCT5 = T5[I];TIME = OBS2[I];OUTPUT;END;KEEP CONCT1 CONCT2 CONCT3 CONCT4 CONCT5 TIME;RUN;PROC PRINT DATA=CONCVTIME3;VAR CONCT1 CONCT2 CONCT3 CONCT4 CONCT5 TIME;RUN;PROC REG DATA=CONCVTIME3;TITLE '3.16(B-1) BOX-COX';TITLE1 'Y TRANSFORM LAMBDA -0.2';MODEL CONCT1 = TIME;RUN;PROC REG DATA=CONCVTIME3;TITLE '3.16(B-2) BOX-COX';TITLE1 'Y TRANSFORM LAMBDA -0.1';MODEL CONCT2 = TIME;RUN;PROC REG DATA=CONCVTIME3;TITLE '3.16(B-3) BOX-COX';TITLE1 'Y TRANSFORM LAMBDA 0';MODEL CONCT3 = TIME;RUN;PROC REG DATA=CONCVTIME3;TITLE '3.16(B-4) BOX-COX';TITLE1 'Y TRANSFORM LAMBDA 0.1';MODEL CONCT4 = TIME;RUN;PROC REG DATA=CONCVTIME3;TITLE '3.16(B-5) BOX-COX';TITLE1 'Y TRANSFORM LAMBDA 0.2';MODEL CONCT5 = TIME;RUN;*3.16(C) LOG-10 TRANSFORMATION;PROC REG DATA=CONCVTIME OUTEST=EST;TITLE '3.16(C) Y TRANSFORM LOG 10';MODEL CONCLOG10 = TIME;OUTPUT OUT=OUTPLOT PREDICTED=PRED RESIDUAL=RESID;RUN;*3.16(D) LOG-10 TRANSFORM DATA AND REG PLOT;SYMBOL1 V=STAR C=BLACK;SYMBOL2 V=DOT C=BLACK I=R;PROC GPLOT DATA=OUTPLOT;TITLE '3.16(D) LOG-10 TRANSFORM';TITLE1 'SCATTER PLOT AND REGRESSION LINE';PLOT CONCLOG10*TIME PRED*TIME/ OVERLAY;RUN;*3.16(E) LOG-10 TRANSFORM;* RESIDUALS AND NORMAL PROBABILITY PLOT;SYMBOL1 V=STAR C=BLACK;PROC GPLOT DATA=OUTPLOT;TITLE '3.16(E-1) Y LOG-10 TRANSFORM';TITLE1 'RESIDUALS VERSUS PREDICTED';PLOT RESID*PRED;RUN;PROC SORT DATA=OUTPLOT;BY RESID;RUN;DATA OUTPLOT;SET OUTPLOT NOBS=NOBS;QUANTILE = PROBIT( (_N_- (3/8)) / (NOBS + (1/4)) );RUN;DATA OUTPLOT2;IF _N_ = 1 THEN SET EST;SET OUTPLOT;EXPRESIDUAL = _RMSE_*QUANTILE;RUN;PROC GPLOT DATA=OUTPLOT2;TITLE '3.16(E-2) Y LOG-10 TRANSFORM';TITLE1 'RESIDUALS VS EXPECTED RESIDUALS';PLOT RESID*EXPRESIDUAL;RUN;QUIT;

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!