Diff_GridOperator: Difference operator for a grid

In James-Thorson/movement_tools: Tools for movement operators in GMRF models

Usage

Diff_GridOperator(var = "x", grid_dataframe, h = 1)

Arguments

var
grid_dataframe
h

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (var = "x", grid_dataframe, h = 1) 
{
    DeltaX = mean(diff(sort(unique(grid_dataframe[, "x"]))))
    DeltaY = mean(diff(sort(unique(grid_dataframe[, "y"]))))
    d <- approx_equals(as.matrix(dist(grid_dataframe)), c(x = DeltaX, 
        y = DeltaY)[var]) * 0.5 * outer(grid_dataframe[[var]], 
        grid_dataframe[[var]], "-")
    diag(d) <- -colSums(d)
    as(d/h, "dgCMatrix")
  }

James-Thorson/movement_tools documentation built on May 7, 2019, 10:19 a.m.