NEWS
fastLISA 1.0.1
- Vignette: corrected the "Reproducibility" note, which incorrectly stated that
multi-core permutation results "may differ between runs". Permutation inference
is thread-invariant: because the RNG is re-seeded per observation from
iseed
and the observation index alone, the pseudo-p-values and z-scores are
bit-for-bit identical for any n.cores and under any OpenMP schedule.
fastLISA 1.0.0 (2026-07-06)
- First release.
- Local Indicators of Spatial Association computed by a plain-C backend with
optional OpenMP multi-threading and a
xoshiro256++ permutation RNG:
local_moran() -- univariate local Moran's I
local_moran_bv() -- bivariate local Moran's I
local_moran_eb() -- Empirical-Bayes-rate local Moran's I
local_geary() -- univariate local Geary's C
local_multigeary() -- multivariate local Geary's C
local_g() / local_gstar() -- Getis-Ord G and G*
- Accepts any
spdep listw weights object, including custom and
non-contiguity (distance-decay) weights.
- Returns compact statistic-specific matrices containing the observed
statistic, permutation Z-score, and pseudo-p-value, with optional
permutation-moment columns.
- Folded pseudo-p-values are returned for Moran, G, and G* statistics;
Geary statistics return tail-adaptive pseudo-p-values.
p.method selects the permutation pseudo-p-value rule: the standard "count"
(default) counts permutations at least as extreme as the observed value, while
"rank" uses spdep's ties-averaged rank. The two differ only under exact
ties; both return a folded (smaller-tail) value.
- Cluster codes follow
rgeoda conventions, including an Isolated category
for observations with no neighbours.
n.cores defaults to 1L; raise it to use multiple OpenMP threads.
- The permutation RNG is re-seeded per observation, so for a fixed
iseed the
pseudo-p-values are identical for any n.cores (verified in
tests/multicore-reproducible.R).