topic Re: Regression with restricted cubic splines in Graphics Programming

Regression with restricted cubic splines

mmazariegos — Wed, 26 Jul 2017 22:19:08 GMT

Hi all,

I want to graph the spline functions, so I create a data set to store spline basis functions, and I got the data set that is below. The issue is when I want to call the spl1-spl4 variables (In the data set appear like spl 1 spl 2 spl 3 and spl 4). The program doesn't recognize any form to call them (nor spl 1 spl 2 spl 3 spl4 neither spl_1 spl_2 spl_3 spl_4). Can anyone help me?


/* Create SplineBasis = data set that contains spline basis functions */
proc glmselect data=merged0811dersesdiet outdesign(addinputvars fullmodel)=SplineBasis;
   effect spl = spline(lactationbypar / naturalcubic basis=tpf(noint) knotmethod=percentiles(5));
   model weightchange2 = spl / selection=none;  
quit;


proc sgplot data=SplineBasis;
   series x=lactationbypar y=spl1;
   series x=lactationbypar y=spl2;
   series x=lactationbypar y=spl3;
   series x=lactationbypar y=spl4;
run;

Thanks,

Re: Regression with restricted cubic splines

Reeza — Wed, 26 Jul 2017 22:59:12 GMT

Run a proc contents and verify the variable name and label. You need the name. If the name happens to have spaces, then wrap it in qutoes and add an N at the end.

'SP 1'n

@mmazariegos wrote:

Hi all,

I want to graph the spline functions, so I create a data set to store spline basis functions, and I got the data set that is below. The issue is when I want to call the spl1-spl4 variables (In the data set appear like spl 1 spl 2 spl 3 and spl 4). The program doesn't recognize any form to call them (nor spl 1 spl 2 spl 3 spl4 neither spl_1 spl_2 spl_3 spl_4). Can anyone help me?
/* Create SplineBasis = data set that contains spline basis functions */
proc glmselect data=merged0811dersesdiet outdesign(addinputvars fullmodel)=SplineBasis;
   effect spl = spline(lactationbypar / naturalcubic basis=tpf(noint) knotmethod=percentiles(5));
   model weightchange2 = spl / selection=none;  
quit;


proc sgplot data=SplineBasis;
   series x=lactationbypar y=spl1;
   series x=lactationbypar y=spl2;
   series x=lactationbypar y=spl3;
   series x=lactationbypar y=spl4;
run;
Thanks,

Re: Regression with restricted cubic splines

mmazariegos — Thu, 27 Jul 2017 03:44:44 GMT

Thank you very much Reeza, It works!

'spl 1'n