GMRF_RW: Random walk GMRF

In andrewzm/atminv: Reproducible package for Section 4.1 of "Spatio-temporal bivariate statistical models for atmospheric trace-gas inversion"

Description

This function initialises a random walk and represents it as a Gaussian Markov Random Field with mean mu and precision matrix Q. Only random walks along the real line, first-order and second-order variants are implemented for now. Also a data frame can be specified with more details on the GMRF.

Usage

GMRF_RW(n = 10, order = 1, precinc = 1, df = data.frame())

Arguments

n	number of vertices
order	1 or 2, depending on the order of the random walk
precinc	precision constant (multiples the template precision matrix)
df	data frame of length n with more details (for example axis, covariate information)

Value

Object of class GMRF with zero mean

Examples

require(Matrix)
# Create a first-order random walk GMRF
my_RW <- GMRF_RW(n=10, order=1, precinc =2)
print(getPrecision(my_RW))
print(getMean(my_RW))
print(getDf(my_RW))

andrewzm/atminv documentation built on May 10, 2019, 11:14 a.m.