Detecting multiple disease clusters using the
information criterion and scan statistic approach developed by
Takahashi and Shimadzu (2020)
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.packages("multiflexscan")
Development version from GitHub:
install.packages("devtools")
devtools::install_github("tkhrotn/multiflexscan")
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)