multiflexscan

Analyze spatial count data for detecting multiple disease clusters using the information criterion and scan statistic approach developed by Takahashi and Shimadzu (2018, 2020). The package builds on rflexscan for candidate cluster generation and adds information-criterion-based selection of the number of clusters and a global Monte Carlo test.

Install

install.packages("multiflexscan")

Development version from GitHub:

install.packages("devtools")
devtools::install_github("tkhrotn/multiflexscan")

Basic usage

library(multiflexscan)

fit <- multiflexscan(
  x = x, y = y, name = area_id,
  observed = observed, expected = expected,
  nb = neighbors,
  scanmethod = "FLEXIBLE",
  stattype = "RESTRICTED",
  clustertype = "HOT",
  clustersize = 10,
  maxclusters = 10,
  ralpha = 0.2,
  simcount = 999,
  cores = 2
)

print(fit)
summary(fit)
plot(fit)

multiflexscan() returns the candidate clusters generated by rflexscan, the information-criterion values used to select the number of clusters, and a Monte Carlo p-value for the selected cluster set. If an sf object containing the same regions is available, selected clusters can be mapped with:

choropleth(fit, regions = regions_sf)

References

• Takahashi K. and Shimadzu H. (2018). Multiple-cluster detection test for purely temporal disease clustering: integration of scan statistics and generalized linear models. PLoS ONE 13(11):e0207821.
• Takahashi K. and Shimadzu H. (2020). Detecting multiple spatial disease clusters: information criterion and scan statistic approach. Int J Health Geogr 19, 33.
• Otani T. and Takahashi K. (2021). Flexible scan statistics for detecting spatial disease clusters: The rflexscan R package. J Stat Softw 99(13), 1–22.