Skip to contents

gllvmTMB currently fits isotropic SPDE/GMRF spatial fields. Consequently, these helpers draw a circle whose radius is the fitted practical range, sqrt(8) / kappa. Equal axes describe the model assumption (H = I); they are not evidence that anisotropy was estimated. Scope: the fitted isotropic-range display is covered by focused tests; the broader spatial-family evidence remains partial; directional anisotropy, delta/spatiotemporal fields, and new likelihoods are rejected from this plotting contract.

Usage

plot_anisotropy(object, return_data = FALSE)

plot_anisotropy2(object, model = 1L)

Arguments

object

A fitted gllvmTMB_multi model containing an intercept-only or random-slope spatial field.

return_data

If TRUE, return the range geometry instead of a plot.

model

Retained for compatibility. Only model = 1 is supported because gllvmTMB has one shared spatial range parameter.

Value

plot_anisotropy() returns a ggplot2 object, or a data frame when return_data = TRUE. plot_anisotropy2() draws the same geometry with base graphics and invisibly returns a list containing data, range, kappa, H = diag(2), and anisotropy_estimated = FALSE.

Scope

These functions do not estimate or diagnose directional anisotropy. Delta and spatiotemporal spatial models are not supported by this plotting contract and fail explicitly.