p-spline using sas/* applied method of p-spline */
/* 1st step: load & open data */
%web_drop_table(WORK.arwa);
FILENAME REFFILE 'D:/Rdata/arwanew.xlsx';
PROC IMPORT DATAFILE=REFFILE dbms=xlsx
out=work.arwa; getnames=yes;
run;
PROC CONTENTS DATA=WORK.arwa;
run;
%web_open_table(WORK.arwa);
/* penalized spline first method */
proc transreg data=WORK.arwa;
model identity (yi) = pbspline(xi); /* AICC criterion */
run;
proc loess data=WORK.arwa;
model yi = xi;
run;
/* penalized spline second method */
proc tpspline data=WORK.arwa;
model yi = (xi);
run;
/* plot results */
proc sgplot data=WORK.arwa;
pbspline x=xi y=yi;
loess x=xi y=yi / nomarkers; /* uses PRESEARCH option */
run;
Exploratory Data Analysis
