topic Re: Visualizing result from regression with restricted cubic splines for longitudinal continuous dat in Statistical Procedures

Visualizing result from regression with restricted cubic splines for longitudinal continuous data

Ming_Yen_Lin — Tue, 14 Mar 2023 02:53:21 GMT

I try to use proc glimmix to model the relationship of one dependent variable and one independent variable with restricted cubic splines (RANGEFRACTIONS(0.05 0.275 0.5 0.725 0.95)). The dependent variable is one repeated continuous outcome. From previous posts, I understand it is not intuitive when we would display the effect of an independent variable with restricted cubic splines through parameters estimated by regression on one dependent variable. Are there any online experts who can offer advice on presenting such a relationship?

Re: Visualizing result from regression with restricted cubic splines for longitudinal continuous dat

acordes — Tue, 14 Mar 2023 11:24:50 GMT

Look at these blogs from Rick Wiclin:

https://blogs.sas.com/content/iml/2020/06/08/interactions-spline-regression.html

Visualize a regression with splines - The DO Loop (sas.com)

Re: Visualizing result from regression with restricted cubic splines for longitudinal continuous dat

StatDave — Wed, 22 Mar 2023 23:37:16 GMT

I don't think this is a problem. You could proceed similarly to as discussed and illustrated for a logistic model in this note - see the latter section titled "Odds Ratio for Increase in Predictor in a Spline Model". While GLIMMIX does not have an EFFECTPLOT statement, you can simply save the model with a STORE statement and then use the saved model in PROC PLM which has an EFFECTPLOT statement. Using example longitudinal data for the Getting Started: Generalized Estimating Equations section of the GENMOD documentation, the following displays the fitted model using a random intercepts model from GLIMMIX.

proc glimmix data=six;
effect sage=spline(age/naturalcubic);
model wheeze =  sage / dist=bin s;
random intercept / subject=case ;
store glim;
run;
proc plm restore=glim;
effectplot / link;
run;

Re: Visualizing result from regression with restricted cubic splines for longitudinal continuous dat

Ming_Yen_Lin — Fri, 24 Mar 2023 09:11:21 GMT

Thank you for the valuable suggestions. It helps me a lot!
My further question is
As you provided SAS codes:
I would reflect the modeling results through your codes. Should I just put "age" or "sage" in your codes to reflect the natural cubic spline relationship between age and wheeze?

proc plm restore=glim; effectplot fit(x=sage) / link; run;

Re: Visualizing result from regression with restricted cubic splines for longitudinal continuous dat

StatDave — Fri, 24 Mar 2023 14:00:05 GMT

Yes, you can add the FIT option using the original variable, AGE, as you showed to plot the fitted spline model. But specifying FIT isn't necessary if there is only a single predictor in the model.