topic Re: Proc transreg in Statistical Procedures

Proc transreg

Abimal_Zippi — Thu, 24 Dec 2020 02:09:13 GMT

Can you help how to write the design matrix to a SAS data set using proc transreg to obtain a basis spline function and also how to specify an interaction between spline(x) and a variable?

Re: Proc transreg

WarrenKuhfeld — Thu, 24 Dec 2020 05:15:53 GMT

https://documentation.sas.com/?docsetId=statug&docsetTarget=statug_transreg_details06.htm&docsetVersion=15.2&locale=en

https://documentation.sas.com/?docsetId=statug&docsetTarget=statug_transreg_details28.htm&docsetVersion=15.2&locale=en

Re: Proc transreg

Rick_SAS — Thu, 24 Dec 2020 11:46:25 GMT

If you are willing to use PROC GLMSELECT instead of TRANSREG, see

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

Re: Proc transreg

Abimal_Zippi — Fri, 25 Dec 2020 02:31:24 GMT

Thanks Rick, I have seen the blog and it is very helpful. Also, have seen your other blog on "visualize a regression with spline". First, I would like to output the design matrix to sas dataset using spline function (aim is to get the weighting function - not the spline effects) and then run the interaction with spline(x). I found the effect statement not helpful to do this. Can you suggest codes how to do it using proc GLMSELECT/TRASNREG?

Re: Proc transreg

Rick_SAS — Fri, 25 Dec 2020 11:31:23 GMT

I don't know what "get the weighting function" means. If you post some example data and state what you are hoping to get, that would be helpful.

Re: Proc transreg

Abimal_Zippi — Sun, 27 Dec 2020 23:31:51 GMT

Thanks again, Rick. The aim is to model the cumulative effects of a time-varying exposure as the weighted sum of past doses (e.g. past 5 years exposures)- a weighted cumulative exposure. In this method, weights, will be modeled with cubic regression B-splines. I.e., the spline basis functions will represent the weighting function.

Re: Proc transreg

Rick_SAS — Fri, 01 Jan 2021 13:05:51 GMT

You still haven't posted data. From what you've said, it sounds like you want to run a two-step process. First, model the weights and score that model at the time points in your data. Then, run the regression and use the WEIGHT statement to use the weights from the first step.

Re: Proc transreg

Abimal_Zippi — Tue, 05 Jan 2021 00:08:20 GMT

Happy New Year Rick, find below sample data - for each ID, the exposure time points and dose in the last 5 years -

data have;
input ID exposure_time_yr dose $;
cards;
1   0.1  15 
1   0.4  7.5
1   1.9  10
1   2    100
1   3.1  100
1   3.8  10
1   3.9  140
1   4.1  140
1   4.8  50
2   0.1 75
2   0.5  150
2   0.2  75
2   0.2  50
2   0.3 75
2   0.3 150
2   0.4  25
2   0.5  225
2   1    25
2   3.5  200
2   4.5  200
3   0.4  5
3   1.5  20
3   2.1  25
3   3    200
3   4.9  100

;

Re: Proc transreg

Rick_SAS — Tue, 05 Jan 2021 00:40:29 GMT

1. Why is dose a character variable? I'll assume that is a mistake.

2. How are you using the ID variable in the model?

3. Where is the weighting variable?

I still don't know what you are trying to accomplish. Can you sketch what you want? Are we supposed to be fitting some line or curves to the following scatter plot?

proc sgplot data=have;
 scatter y=exposure_time_yr x=dose / group=ID markerattrs=(symbol=CircleFilled size=12);
run;

Re: Proc transreg

Abimal_Zippi — Tue, 05 Jan 2021 01:17:54 GMT

Yes, the dose is not a character!

As mentioned in the previous threads, I wanted the weighted sum of pas does over a 5-year window - the data shows the dose exposure and time points over 5 year for each ID. Weights, will be modeled with cubic regression B-splines. I.e., the spline basis functions will represent the weighting function.

Re: Proc transreg

Abimal_Zippi — Tue, 05 Jan 2021 01:20:41 GMT

* the weighted sum of past doses over a 5-year window

Re: Proc transreg

Rick_SAS — Tue, 05 Jan 2021 11:15:32 GMT

Sorry, but I do not understand. Do you have a reference (textbook? journal article?) that provides details?

I am giving up. Perhaps someone else will be able to help you. Good luck.

Re: Proc transreg

WarrenKuhfeld — Tue, 05 Jan 2021 11:27:31 GMT

Define what you mean by "weights." Look at examples of the WEIGHT statement in regression procs. Is that what you mean? Or do you mean regression coefficients? The TRANSREG documentation has examples of scoring as do many other procedures.

Re: Proc transreg

Abimal_Zippi — Tue, 05 Jan 2021 12:27:11 GMT

Thanks Rick, here is the link for the paper that I am referring to :

Modeling of cumulative effects of time-varying drug exposures on within-subject changes in a continuous outcome - Coraline Danieli, Therese Sheppard, Ruth Costello, William G Dixon, Michal Abrahamowicz, 2020 (sagepub.com)

Re: Proc transreg

Abimal_Zippi — Tue, 05 Jan 2021 12:40:13 GMT

In the paper mentioned above, " the weight function assigns the differential importance weights to past doses, depending on the time elapsed since the dose was taken" - weighted according to their relative proximity in time.