Hello all,
I have a longitudinal dataset to evaluate the effect of a vaccine in a population. I want to perform an ITS analysis with April 2010 as the intervention date.
I used the command --------- itsa incidence, single trperiod (2010m4) replace prais nolog figure ------ which renegerated a really interesting graph. However, in STATA, the models available are the Prais-Winsten and Cochrane-Orcutt regression or the Newey-West models.
I want to run the analysis using the negative binomial regression model. The analysis has been performed in both SAS and R; both generated same results (change in level and change in trend).
I hope to use STATA to generate same results. Any help will be greatly appreciated. Thanks!!!!
Please find dataset below-
I have a longitudinal dataset to evaluate the effect of a vaccine in a population. I want to perform an ITS analysis with April 2010 as the intervention date.
I used the command --------- itsa incidence, single trperiod (2010m4) replace prais nolog figure ------ which renegerated a really interesting graph. However, in STATA, the models available are the Prais-Winsten and Cochrane-Orcutt regression or the Newey-West models.
I want to run the analysis using the negative binomial regression model. The analysis has been performed in both SAS and R; both generated same results (change in level and change in trend).
I hope to use STATA to generate same results. Any help will be greatly appreciated. Thanks!!!!
Please find dataset below-
Code:
* Example generated by -dataex-. For more info, type help dataex clear input double(pneumonia_age_14days person_month) float date 190 1158653.355263158 576 184 1088046.5460526315 577 166 1170664.4407894737 578 138 1141594.9342105263 579 121 1187655.2631578948 580 133 1155605.9210526315 581 115 1138360.3947368423 582 74 1143168.4868421054 583 128 1111242.335526316 584 195 1152219.1447368423 585 348 1119342.6644736843 586 413 1163404.0460526317 587 222 1174004.5723684211 588 174 1066112.7960526317 589 210 1184592.0394736843 590 150 1153103.75 591 132 1198421.6447368423 592 164 1165814.6381578948 593 119 1147966.5460526315 594 78 1151160.3947368423 595 185 1117849.0460526317 596 295 1152416.8092105263 597 344 1113918.3881578948 598 453 1156184.2763157897 599 226 1164500.8223684211 600 245 1053057.730263158 601 240 1171506.6776315789 602 181 1138432.5986842108 603 140 1178651.4473684211 604 112 1142964.6381578948 605 80 1124042.0394736843 606 64 1127701.9736842106 607 121 1095950.657894737 608 187 1135917.3684210526 609 263 1105166.1513157897 610 417 1148506.3815789474 611 219 1155725.9539473685 612 213 1047862.8289473683 613 231 1167900 614 177 1138079.7039473685 615 end format %tm date
