A Bayesian Hierarchical Spatial Model to Correct for Misreporting in Count Data: Application to State-Level COVID-19 Data in the United States

Bayesian hierarchical spatiotemporal models for prediction of (under)reporting rates and cases: COVID-19 infection among the older people in the United States during the 2020-2022 pandemic

The gap between the reported and actual COVID-19 infection cases has been an issue of concern. Here, we present Bayesian hierarchical spatiotemporal disease mapping models for state-level predictions of COVID-19 infection risks and (under)reporting rates among people aged 65 and above during the first two years of the pandemic in the United States. With prior elicitation based on recent prevalence studies, the study suggests that the median state-level reporting rate of COVID-19 infection was 90% (interquartile range: [78%, 96%]). Our study uncovers spatiotemporal variations and dynamics in state-level infection risks and (under)reporting rates, suggesting time-varying associations between higher population density, higher percentage of minorities, and higher percentage of vaccination and increased risks of COVID-19 infection, as well as an association between more easily accessible tests and higher reporting rates. With sensitivity analyses, we highlight the impact and importance of incorporating covariates information and objective prior references for evaluating the issue of underreporting.

Visible minority status and occupation were associated with increased COVID-19 infection in Greater Vancouver British Columbia between June and November 2020: an ecological study

Background

The COVID-19 pandemic has highlighted health disparities, especially among specific population groups. This study examines the spatial relationship between the proportion of visible minorities (VM), occupation types and COVID-19 infection in the Greater Vancouver region of British Columbia, Canada.

Methods

Provincial COVID-19 case data between June 24, 2020, and November 7, 2020, were aggregated by census dissemination area and linked with sociodemographic data from the Canadian 2016 census. Bayesian spatial Poisson regression models were used to examine the association between proportion of visible minorities, occupation types and COVID-19 infection. Models were adjusted for COVID-19 testing rates and other sociodemographic factors. Relative risk (RR) and 95% Credible Intervals (95% CrI) were calculated.

Results

We found an inverse relationship between the proportion of the Chinese population and risk of COVID-19 infection (RR = 0.98 95% CrI = 0.96, 0.99), whereas an increased risk was observed for the proportions of the South Asian group (RR = 1.10, 95% CrI = 1.08, 1.12), and Other Visible Minority group (RR = 1.06, 95% CrI = 1.04, 1.08). Similarly, a higher proportion of frontline workers (RR = 1.05, 95% CrI = 1.04, 1.07) was associated with higher infection risk compared to non-frontline.

Conclusion

Despite adjustments for testing, housing, occupation, and other social economic status variables, there is still a substantial association between the proportion of visible minorities, occupation types, and the risk of acquiring COVID-19 infection in British Columbia. This ecological analysis highlights the existing disparities in the burden of diseases among different visible minority populations and occupation types.

Evaluating active leprosy case identification methods in six districts of Nepal

BACKGROUND

Nepal has achieved and sustained the elimination of leprosy as a public health problem since 2009, but 17 districts and 3 provinces with 41% (10,907,128) of Nepal's population have yet to eliminate the disease. Pediatric cases and grade-2 disabilities (G2D) indicate recent transmission and late diagnosis, respectively, which necessitate active and early case detection. This operational research was performed to identify approaches best suited for early case detection, determine community-based leprosy epidemiology, and identify hidden leprosy cases early and respond with prompt treatment.

METHODS

Active case detection was undertaken in two Nepali provinces with the greatest burden of leprosy, Madhesh Province (40% national cases) and Lumbini Province (18%) and at-risk prison populations in Madhesh, Lumbini and Bagmati provinces. Case detection was performed by (1) house-to-house visits among vulnerable populations (n = 26,469); (2) contact examination and tracing (n = 7608); in Madhesh and Lumbini Provinces and, (3) screening prison populations (n = 4428) in Madhesh, Lumbini and Bagmati Provinces of Nepal. Per case direct medical and non-medical costs for each approach were calculated.

RESULTS

New case detection rates were highest for contact tracing (250), followed by house-to-house visits (102) and prison screening (45) per 100,000 population screened. However, the cost per case identified was cheapest for house-to-house visits [Nepalese rupee (NPR) 76,500/case], followed by contact tracing (NPR 90,286/case) and prison screening (NPR 298,300/case). House-to-house and contact tracing case paucibacillary/multibacillary (PB:MB) ratios were 59:41 and 68:32; female/male ratios 63:37 and 57:43; pediatric cases 11% in both approaches; and grade-2 disabilities (G2D) 11% and 5%, respectively. Developing leprosy was not significantly different among household and neighbor contacts [odds ratios (OR) = 1.4, 95% confidence interval (CI): 0.24-5.85] and for contacts of MB versus PB cases (OR = 0.7, 95% CI 0.26-2.0). Attack rates were not significantly different among household contacts of MB cases (0.32%, 95% CI 0.07-0.94%) and PB cases (0.13%, 95% CI 0.03-0.73) (χ2 = 0.07, df = 1, P = 0.9) and neighbor contacts of MB cases (0.23%, 0.1-0.46) and PB cases (0.48%, 0.19-0.98) (χ2 = 0.8, df = 1, P = 0.7). BCG vaccination with scar presence had a significant protective effect against leprosy (OR = 0.42, 0.22-0.81).

CONCLUSIONS

The most effective case identification approach here is contact tracing, followed by house-to-house visits in vulnerable populations and screening in prisons, although house-to-house visits are cheaper. The findings suggest that hidden cases, recent transmission, and late diagnosis in the community exist and highlight the importance of early case detection.

Estimating state-specific population-based hospitalization rates from in-state hospital discharge data

OBJECTIVE

To develop weights to estimate state population-based hospitalization rates for all residents of a state using only data from in-state hospitals which exclude residents treated in other states.

DATA SOURCES AND STUDY SETTING

Agency for Healthcare Research and Quality, Healthcare Cost and Utilization Project (HCUP), State Inpatient Databases (SID), 2018-2019, 47 states+DC.

STUDY DESIGN

We identified characteristics for patients hospitalized in each state differentiating movers (discharges for patients hospitalized outside state of residence) from stayers (discharges for patients hospitalized in state of residence) and created weights based on 2018 data informed by these characteristics. We calculated standard errors using a sampling framework and compared weight-based estimates against complete observed values for 2019.

DATA COLLECTION/EXTRACTION METHODS

SID are based on administrative billing records collected by hospitals, shared with statewide data organizations, and provided to HCUP.

PRINCIPAL FINDINGS

Of 34,186,766 discharged patients in 2018, 4.2% were movers. A higher share of movers (vs. stayers) lived in state border and rural counties; a lower share had discharges billed to Medicaid or were hospitalized for maternal/neonatal services. The difference between 2019 observed and estimated total discharges for all included states and DC was 9402 (mean absolute percentage error = 0.2%). We overestimated discharges with an expected payer of Medicaid, from the lowest income communities, and for maternal/neonatal care. We underestimated discharges with an expected payer of private insurance, from the highest income communities, and with injury diagnoses and surgical services. Estimates for most subsets were not within a 95% confidence interval, likely due to factors impossible to account for (e.g., hospital closures/openings, shifting consumer preferences).

CONCLUSIONS

The weights offer a practical solution for researchers with access to only a single state's data to account for movers when calculating population-based hospitalization rates.