ngme space-time non-separable model specification

Source: R/complex-models.R 
spacetime(
  mesh,
  lambda = 1,
  alpha = 2,
  cc = 1,
  kappa = 1,
  fix_gamma = FALSE,
  theta_gamma_x = 0,
  theta_gamma_y = 0,
  shared_theta_gamma = FALSE,
  B_gamma_x = matrix(1, nrow = mesh[[2]]$n, ncol = 1),
  B_gamma_y = matrix(1, nrow = mesh[[2]]$n, ncol = 1),
  B_gamma_x_list = NULL,
  B_gamma_y_list = NULL,
  stabilization = TRUE,
  ...
)

Arguments

mesh

A list of two objects:

  • mesh_t - The temporal mesh

  • mesh_s - The spatial mesh

lambda

The spatial damping parameter.

alpha

2 or 4, SPDE smoothness parameter.

cc

Parameter c in the SPDE.

kappa

Kappa parameter from Matern SPDE.

fix_gamma

TRUE if fix gamma (advection term), FALSE if estimate gamma.

theta_gamma_x

The x component of the advection term: `gamma_x = B_gamma_x

theta_gamma_yThe y component of the advection term: `gamma_y = B_gamma_y shared_theta_gammaTRUE if share the same theta_gamma for all time nodes. (theta_gamma_x and theta_gamma_y will be the same)B_gamma_xThe design matrix for the x component of the advection term.B_gamma_yThe design matrix for the y component of the advection term.B_gamma_x_listA list of design matrices for the x component of the advection term on every time node, length(B_gamma_x_list) == nt-1.B_gamma_y_listA list of design matrices for the y component of the advection term on every time node, length(B_gamma_y_list) == nt-1.stabilizationTRUE if using a stabilization term (for implicit Euler)....Additional arguments (ignored). ngme_operator object. Given a spatial and temporal model, build a non-separable space-time model.