Hi.
Just looking for some advice/help.
We've run an interrupted time series analysis and wantvto report the difference between the point estimate and counterfactual at 1 year and 2 years post-intervention.
We want to report such with a p-value and CIs. But such needs to be derived from post-hoc estimation. This is currently done using the emmeans package in R. Said package offers a pairwise comparison with Tuckey adjusted p-value. I wondered whether directly reporting the above p-value was the right approach or can you use that based on the z-score of the estimate and work out the p manually? The former [Tuckey adjusted] yields non-significant results whereas a manual value derived from the z-score of the output is statistically significant.
Assume we're not the first to attempt to report such and wondered if anyone could share how they derived p-values and CIs previously for a point estimate to counterfactual comparison at a specific time point.