How to plot Restricted Cubic Spline in PROC LOGISTIC (BY IMPUTATION)

[asgee]

I'm trying to run an analysis where I have a continuous variable (serum) and binary outcome "par" (yes/no). 

My analyses requires that I impute my dataset (50 iterations). I'm still quite new to visualizing plots, and am having trouble trying to visualize the spline effects of my logistic model. Here's my current code:

title "Restricted Splines";
title2 "Four Internal Knots";
		ods select ParameterEstimates SplineKnots;
		proc logistic data=imputed_50; 
		effect spl = spline(serum / details naturalcubic basis=tpf(noint)                 
                knotmethod=percentilelist(5 35 65 95)); /* RESTRICTED CUBIC SPLINE BASED ON 4 KNOTS */ 
		class par (ref='0') dev (ref='0') eth (ref='1') product (ref='0') / param=ref; 
		model par (event='1')= spl age weight eth dev product / selection=none covb;  
		by _Imputation_;     /* RUNNING 50 ITERATIONS OF THIS MODEL BASED ON AN IMPUTED DATASET */ 
		ods output ParameterEstimates=Lgsparms; 
		run;

What I'm trying to replicate is this graph based on the SAS documentation of Visualizing regressions with splines:

The example above is obviously using a continuous dependent variable (MPG). I'm having trouble trying to show the spline effect / spline points of "serum" based on a binary outcome "par", particularly since my modelling is based on an imputation and not just one iteration. 

I understand there's the "effectplots" option in PROC LOGISTIC but not sure how to implement that in my modelling that is based on an imputed dataset. Any help plotting this would be very much appreciated. 

[Rick_SAS]

Do you really want to view 50 plots, one for each imputation? Typically you will want to suppress all output (tables and graphs) by using the NOPRINT option during the estimation of parameters for each set of imputed values. You then run PROC MIANALYZE to aggregate the 50 estimates into one "best" set of parameter estimates. You can then visualize the predicted probabilities for the final model.

[asgee]

Hi @Rick_SAS ! Yes I forgot to mention that. That's the next sections in my code (I run proc mianalyze and aggregate the 50 estimates). The "ParameterEstimates=Lgsparms" is that output that I sort of clean later down the line and produce one table that shows the "best" parameter estimates using the MODELEFFECTS statement:

   ods trace on;
		proc mianalyze parms=Lgsparms; 
		modeleffects Intercept spline1 spline2 spline3 age weight eth dev product; /* 3 spline effects based on 4 knots */
		run;
   ods trace off;
   ods listing close;
   ods output ParameterEstimates=mi_results

From there I use the "mi_results" output to produce a table summarizing that one "best" set of parameter estimates.

As you pointed out, my trouble really is just trying to visualize the predicted probabilities for that aggregated "best" set. I'm not sure exactly how to isolate the predicted probabilities (either through PROC LOGISTIC or PROC MIANALYZE) for that aggregated "best" set. I tried adding a "predicted=Fit" statement beside the (ods output ParameterEstimates=Lgsparms) line in PROC LOGISTIC but it gave me an error instead.

I understand that if I have an output that just has the predicted probabilities for that final model, I can use the SGPLOT (Or EFFECTPLOTS??) function to visualize those results. Not sure if I'm missing an option or a data step to get to that point...

[Rick_SAS · Accepted solution]

It sounds like you have three problems:

1. You need to evaluate the final model, which is defined by the parameter estimates table.

The logistic model is 

eta = Intercept + b1*spline1 + b2*spline2 + ... + b8*product

prob = logistic(eta)

2. Because your model is defined in terms of splines, you should output the design matrix, which will contain the spline1-spline3 variables.

3. You need to decide how to visualize a surface in 8 variables. The typical way is to create a slice plot in which you fix a value for the age weight eth dev and product variables (often a mean or median value) and then plot the predicted probability versus the serum variable. The general idea of a slice plot is explained in this blog post: https://blogs.sas.com/content/iml/2017/12/18/visualize-multivariate-regression-models-by-slicing-con...

I think you might need to create the plot manually.

There are a lot of technical details involved in doing this so I recommend asking if you really need to do this visualization. It might be sufficient to plot the predicted probability for each subject in the data, which is much easier.

[asgee]

Hi @Rick_SAS , thanks for your reply. The steps you've outlined makes sense - I think I'll just continue exploring other options for now and follow-up with your suggestion of plotting the pred probabilities for each subject in the data. I'll have a look at the design matrix link you've sent, seems like that's a good first step for now.

Thanks!

[Saifulinfs]

hello

I am following this discussion as its helpful. Did you figure out how to get the predicted probabilities for the final model?