am2 <- assess(formula= los ~ ., data=hosprog, intervention = "program",
topcode =NULL, int.time="month", regression="none", treatment= 5,
interrupt=c(5,9), did="many", its="two", newdata=TRUE, propensity=NULL)
plot(am2, "DID", add.legend="bottomleft", ylim=c(2, 8)) #DID model
plot(am2, "ITS", add.legend="top", ylim=c(2, 8)) #ITS model
Run the code above in your browser using DataLab