The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

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

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.