setGMRF: Set GMRF

In shazhe/mvst0: Multi-variate ST modelling on globe

Description

This function takes an object of class Graph_2nodes and changes the GMRF (process) component accordingly. This can be used to change the mean in of the process component, or the precision matrix when chaining Gaussian updates.

Usage

setGMRF(Graph, obj)

## S4 method for signature 'Graph_2nodes,GMRF'
setGMRF(Graph, obj)

Arguments

Graph	object of class Graph_2nodes.
obj	object of class GMRF.

Value

Object of class Graph_2nodes.

Examples

## Not run: 
require(Matrix)
data(icesat)
data(surf_fe)

## First create observation object
icesat_obs <- Obs(df=icesat,
                 abs_lim = 5,
                 avr_method = "median",
                 box_size=100,
                 name="icesat")

## Now create GMRF defined over some FE basis
Mesh <- initFEbasis(p=surf_fe$p,
                    t=surf_fe$t,
                    M=surf_fe$M,
                    K=surf_fe$K)

mu <- matrix(0,nrow(Mesh),1)
Q <- sparseMatrix(i=1:nrow(surf_fe$p), j = 1:nrow(surf_fe$p), x = 1)

my_GMRF <- GMRF(mu = mu, Q = Q,name="SURF",t_axis = 0:6)
SURF <-GMRF_basis(G = my_GMRF, Basis = Mesh)

L1 <- link(SURF,icesat_obs)
e <- link_list(list(L1))
v <- block_list(list(O = icesat_obs, G = SURF))
G <- new("Graph",e=e,v=v)
G_reduced <- compress(G)

mu2 <- matrix(1,nrow(Mesh),1)
my_GMRF2 <- GMRF(mu = mu2, Q = Q,name="SURF",t_axis = 0:6)

G_reduced <- setGMRF(G_reduced, obj = my_GMRF2)


## End(Not run)

shazhe/mvst0 documentation built on May 29, 2019, 9:20 p.m.