GET (version 0.1-3)

qdir_envelope: Global scaled maximum absolute difference (MAD) envelope tests

Description

Performs the global scaled MAD envelope tests, either directional quantile or studentised, or the unscaled MAD envelope test. These tests correspond to calling the function global_envelope_test with type="qdir", type = "st" and type="unscaled", respectively. The functions qdir_envelope, st_envelope and unscaled_envelope have been kept for historical reasons; preferably use global_envelope_test with the suitable type argument.

Usage

qdir_envelope(curve_set, ...)
st_envelope(curve_set, ...)
unscaled_envelope(curve_set, ...)

Arguments

curve_set

A curve_set (see create_curve_set) or an envelope object. If an envelope object is given, it must contain the summary functions from the simulated patterns which can be achieved by setting savefuns = TRUE when calling envelope.

...

Additional parameters to be passed to global_envelope_test.

Value

An object of class "global_envelope" and "fv" (see fv.object), which can be printed and plotted directly. See global_envelope_test for more details.

Details

The directional quantile envelope test (Myllym<U+00E4>ki et al., 2015, 2017) takes into account the unequal variances of the test function T(r) for different distances r and is also protected against asymmetry of T(r).

The studentised envelope test (Myllym<U+00E4>ki et al., 2015, 2017) takes into account the unequal variances of the test function T(r) for different distances r.

The unscaled envelope test (Ripley, 1981) corresponds to the classical maximum deviation test without scaling, and leads to envelopes with constant width over the distances r. Thus, it suffers from unequal variance of T(r) over the distances r and from the asymmetry of distribution of T(r). We recommend to use the other global envelope tests available, see global_envelope_test for full list of alternatives.

References

Myllym<U+00E4>ki, M., Grabarnik, P., Seijo, H. and Stoyan. D. (2015). Deviation test construction and power comparison for marked spatial point patterns. Spatial Statistics 11: 19-34. doi: 10.1016/j.spasta.2014.11.004

Myllym<U+00E4>ki, M., Mrkvi<U+010D>ka, T., Grabarnik, P., Seijo, H. and Hahn, U. (2017). Global envelope tests for spatial point patterns. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 79: 381<U+2013>404. doi: 10.1111/rssb.12172

Ripley, B.D. (1981). Spatial statistics. Wiley, New Jersey.

See Also

global_envelope_test

Examples

# NOT RUN {
# See more examples in ?global_envelope_test
## Testing complete spatial randomness (CSR)
#-------------------------------------------
if(require("spatstat", quietly=TRUE)) {
  pp <- spruces
  
# }
# NOT RUN {
nsim <- 999 # Number of simulations
# }
# NOT RUN {
  
# }
# NOT RUN {
  ## Test for complete spatial randomness (CSR)
  # Generate nsim simulations under CSR, calculate L-function for the data and simulations
  env <- envelope(pp, fun="Lest", nsim=nsim, savefuns=TRUE, correction="translate",
                  simulate=expression(runifpoint(ex=pp)))
  res_qdir <- qdir_envelope(env) # The directional quantile envelope test
  plot(res_qdir)
  # or (requires R library ggplot2)
  plot(res_qdir, plot_style="ggplot2")

  ## Advanced use:
  # Create a curve set, choosing the interval of distances [r_min, r_max]
  curve_set <- crop_curves(env, r_min=1, r_max=8)
  # For better visualisation, take the L(r)-r function
  curve_set <- residual(curve_set, use_theo=TRUE)
  # The directional quantile envelope test
  res_qdir <- qdir_envelope(curve_set); plot(res_qdir, plot_style="ggplot2")
  # The studentised envelope test
  res_st <- st_envelope(curve_set); plot(res_st, plot_style="ggplot2")
  # The unscaled envelope test
  res_unscaled <- unscaled_envelope(curve_set); plot(res_unscaled, plot_style="ggplot2")
}
# }