I have a regression model using discrete time hazard regression. And I wanted to study the effect of modifier A (binary) and factor B(linear spline with 5 knots, represents by “B?" in below), and the interaction between A and B (represents by "B?int?" below) to the outcome y. There was 3 follow-up timepoints, so fu1, fu2, fu3 were used in the discrete time model.

Therefore, the regression writes: cloglog y i.A B? B?int? fu1 fu2 fu3, nocon eform

Then I want to calculate the estimates for 100 points, and use these estimates to draw a graph. And my question is highlighted in red below: I wonder to represent the coefficients for fu1, fu2, fu3 in mata structure, should I use m0? To simply my question, how to represent estimates in mata for the “covariates”.

Here is the code I wrote:

set obs 100

mata:xw=st_data(.,”B?")
mata:xn=J(100,5,0)
mata:m0=J(100,1,0)
mata:m1=J(100,1,1)

mata:x0=(xw,m0,m0,xn,m0,m0,m0) // reference group
mata:x1=(xw,m0,m0,xw,m0,m0,m0) // group 1 in A, without main effect
mata:x2=(xw,m0,m1,xw,m0,m0,m0) // group 1 in A, with main effect
mata:x3=(xn,m0,m0,xw,m0,m0,m0) // group 1 interaction

Thank you very much!