ssim_polygon(), using a separate
permutation mean for each polygon and statistic and the existing draws
that independently permute both maps._p, _q, and _sig
columns for all four local statistics. Apply BH correction across
polygons separately for each statistic. Preserve p_value,
q_value, and sig as aliases for the SSIM
results; with fdr = FALSE, _q values equal the
unadjusted p-values.q <= alpha for local polygon significance,
matching rasters; without FDR adjustment use p <= alpha.
The observed statistics, existing SSIM p/q values, global tests, and
random-number sequence are unchanged for the same inputs and seed.transform to "none" in
ssim_raster(), ssim_polygon(), and
ssim_bandwidth(). Existing analyses using an explicit
transformation retain that transformation.Require R >= 3.6.0 so the bundled fixed-seed reproduction uses
the documented Rejection sampling RNG consistently. The
exported functions retain their numerical calculations and
defaults.
Include portable, offline manuscript reproduction scripts and
compact reference values accessible with
system.file("reproduction", package = "SSIMmap"). Fixed
seeds reproduce the tables and local significance maps; the commuting
comparison has 124 of 279 areas with q < 0.01.
Document and independently verify the November 11 FWI against the official CWFIS historical coverage, including source attribution and Open Government Licence - Canada information.
Add fwi_1111_bc, the actual 11 November 2023 FWI
raster used in the revised manuscript. The original
fwi_1101_bc remains the unchanged 1 November 2023 dataset
for compatibility.
Update the FWI examples to use the November 11 data with explicit
transform = "none" and 3x3, 5x5, or 7x7 windows. The
default transformation and all three numerical algorithms are
unchanged.
Correct the Toronto polygon documentation to 279 Aggregate Dissemination Areas (2021) and identify the CIMD source dimension as residential instability, while preserving the numeric example data. Document the preserved commuting measure using Census characteristic 2601, total Census population as denominator, two-decimal DA rounding, unweighted intersecting-DA averaging, and the suppressed-source-value limitation.
Use the manuscript’s maximum bandwidth of 200 in example selection, yielding midpoint bandwidths 87 and 91.
Make vignette examples portable, use regular panel tags, and correct raster significance-layer names.
Use American spelling in authored documentation and display labels; preserve official source titles, license names, and API argument names. Numerical algorithms are unchanged.
ssim_bandwidth(): redesigned with
transform, option arguments and bias-variance
trade-off plots.ssim_polygon() and ssim_raster(): added
permutation testing (do_test, R,
fdr, alpha, seed) and
FDR-corrected local p-values.ssim_constant(). Constants k1 and
k2 are now arguments of ssim_polygon().Toronto_SSIM, fwi_0816_bc,
fwi_0818_bc, fwi_1101_bc.Toronto_SSIM (sf polygon object).fwi_0816_bc, fwi_0818_bc,
fwi_1101_bc (packed terra rasters of BC Fire Weather Index
for Aug 16, Aug 18, and Nov 1, 2023).