We have a situation where we want to test the association between X and Y, but the change in X from baseline is more interpretable.
There are several possibilities I see for setting up the model. But I would appreciate some help in understanding the implications of each choice, like how the interpretations change.
In these 3 choices, we are modeling our outcome at all time points:
A1. $Y_{ij} = \beta_0 + \beta_1*X_{ij}$
Here we model Y and X, without any baseline adjustment
A2. $Y_{ij} = \beta_0 + \beta_1*(X_{ij}-X_{i0})$
Here we adjust for the change score of $(X_{ij}-X_{i0})$, without any baseline adjustment. At time 0, this change score will equal 0.
A3. $Y_{ij} = \beta_0 + \beta_1*X_{i0}+\beta_2*(X_{ij}-X_{i0})$
Here we adjust for the change score, but also for the baseline level of the covariate, X0.
Looking at the results from A1 and A3 models: $\beta_1 (A1) \approx \beta_2 (A3) \approx 2*\beta_1 (A2).$
But also, $\beta_1 (A3) \approx \beta_2 (A3)$, which I think corresponds to: $\beta_0 + \beta_1*X_{i0}$ at baseline, and $\beta_0 + (\beta_1-\beta_2)*X_{i0} + \beta_2*X_{ij} = \beta_0 + \beta_2*X_{ij}$ after baseline, so essentially, we are just in the A1 case.
It seems that there might be little difference in models A1 and A3, and that it would not actually be correct to interpret $\beta_2(A3)$ as the effect of the change in X on Y, right?
But would it be useful/correct to use an interaction with $(X_{ij}-X_{i0})$ (e.g. $(X_{ij}-X_{i0})*Z$) to test if post-baseline level of X differs across levels of Z? This seems analogous to a using an interaction with a spline, where we might expect things to change more or less after baseline.