options ps=54 ls=76; FILENAME website HTTP "http://www.stat.ufl.edu/~winner/data/h2o_evap.csv"; proc import datafile=website out=evap1 dbms=csv replace; getnames=yes; run; ods PDF file="proj1_part4.pdf"; ods graphics on; proc reg; model evap = tgrad satdefmax wind; output out=regout r=e1; run; proc univariate normal; var e; run; proc transreg; model BoxCox(evap / lambda=-2 to 2 by 0.02) = identity(tgrad satdefmax wind); run; data evap2; set evap1; sqrtevap = sqrt(evap); run; proc reg; model sqrtevap = tgrad satdefmax wind / dw dwprob; run; proc autoreg; model sqrtevap = tgrad satdefmax wind / nlag=1; run; ods graphics off; ods PDF close; run; quit;