topic Re: Segmented line graph in SAS Programming

Segmented line graph

Levi_M — Tue, 26 Jul 2022 17:15:11 GMT

Hi! I am putting together a segmented regression with visualization. However, plotting the data, I can not seem to get the two periods to touch. Could anyone help provide assistance with this?

proc autoreg data=a1;

model rates = QTRs aft_cnt / method=ml nlag=12 backstep dwprob loglikl;;

output out=auto_reg predicted=pred;

by Patient_Age_Category; /*This line of code will treat each age_cat separately*/

quit;

proc sgplot data=glimmix_out noautolegend;

block x=QTRs block=period / transparency=.75;

scatter x=QTRs y=rates / markerattrs=(symbol=CircleFilled color=black);

series x=QTRs y=pred / group=period lineattrs=(thickness=3);

by Patient_Age_Category; /*This line of code will treat each age_cat separately*/

run;

Thank you!!

Re: Segmented line graph

PaigeMiller — Tue, 26 Jul 2022 17:23:10 GMT

I don't know what you mean by this: "I can not seem to get the two periods to touch" ?? What should the plot look like if the two periods touched? (It looks like they touch to me, separated by a vertical line)

Re: Segmented line graph

Levi_M — Tue, 26 Jul 2022 17:29:17 GMT

Thank you for your reply. What I meant to say was the two regression lines.

EDIT: for example, when I run this same process in R, the two regression lines meet at the joint. (with the same results)

Re: Segmented line graph

PaigeMiller — Tue, 26 Jul 2022 17:34:36 GMT

But your blue line and your red line do not intersect each other at the boundary between blue and red. They would intersect to the left of the boundary, in the blue region.

Re: Segmented line graph

Levi_M — Tue, 26 Jul 2022 17:36:32 GMT

I see what you are saying, thanks. I just need to drop that boundary region.

Thanks.

Re: Segmented line graph

svh — Tue, 26 Jul 2022 17:38:07 GMT

My guess is that the regression line in the left panel is only extending across the range of the data in the "before" period. The point on the vertical line would seem to belong to the "after" period.

Re: Segmented line graph

Ksharp — Wed, 27 Jul 2022 13:05:05 GMT

https://blogs.sas.com/content/iml/2017/04/05/nonsmooth-models-spline-effects.html https://blogs.sas.com/content/iml/2020/12/14/segmented-regression-sas.html