How do I preform a interrupted time series analysis using PROC GENMOD with a binary outcome?

Fable — Mon, 03 Jul 2023 03:24:56 GMT

Hi, I am trying to conduct an interrupted time series analysis looking at lab test rates (outcome) and how they changed due to the COVID shelter-in-place order in March 2020.

I am making use of this paper to guide me: Time After Time: Difference-in-Differences and Interrupted Time Series Models in SAS

I have a individual level dataset of patients with multiple records per patient, one for each month they are in the cohort. From what I understand, since I am looking at testing rates, I could either model this with the individual data as a binary outcome (ie have an flag set to 1 if the patient in a particular month had a test and 0 if they did not) OR I could aggregate the patient level data so that I have the overall testing rate as a percentage per month, meaning I could then model it as a continuous outcome.

Does anyone know what the difference would be for each approach?

From what the linked paper makes it seem like, I should use the GENMOD Procedure if I want to do it as a binary outcome, and use the MIXED Procedure if I want to do it as an aggregated continuous outcome.

The paper provides an example of how to use PROC MIXED for a continuous outcome, but how would I use PROC GENMOD for a binary one? The paper says:

These proportions can be estimated using PROC GENMOD as odds or probabilities on the log scale (relative risks) or as proportions (absolute risks) for each group at each time point.

How would one specify if they want relative risks or absolute risks with PROC GENMOD? How does one interpret each?

Thank you.

Re: How do I preform a interrupted time series analysis using PROC GENMOD with a binary outcome?

SteveDenham — Fri, 07 Jul 2023 18:11:12 GMT

This is a case where the MIXED approach can be generalized (see what I did there) by using GLIMMIX with a binomial distribution.

SteveDenham

topic Re: How do I preform a interrupted time series analysis using PROC GENMOD with a binary outcome? in Statistical Procedures

How do I preform a interrupted time series analysis using PROC GENMOD with a binary outcome?

Re: How do I preform a interrupted time series analysis using PROC GENMOD with a binary outcome?