vunit_measure: vunit_measure

In spatPomp: Inference for Spatiotemporal Partially Observed Markov Processes

vunit_measure

Description

vunit_measure evaluates the variance of a unit's observation given the entire state

Usage

## S4 method for signature 'spatPomp'
vunit_measure(object, x, unit, time, params, Np = 1)

Arguments

object	An object of class spatPomp
x	A state vector for all units
unit	The unit for which to evaluate the variance
time	The time for which to evaluate the variance
params	parameters at which to evaluate the unit variance
Np	numeric; defaults to 1 and the user need not change this

Value

A matrix with the unit measurement variance implied by the state, x, and the parameter set params for unit unit.

Author(s)

Kidus Asfaw

Examples

# Complete examples are provided in the package tests
## Not run: 
b <- bm(U=3)
s <- states(b)[,1,drop=FALSE]
rownames(s) -> rn
dim(s) <- c(3,1,1)
dimnames(s) <- list(variable=rn, rep=NULL)
p <- coef(b); names(p) -> rnp
dim(p) <- c(length(p),1); dimnames(p) <- list(param=rnp)
o <- obs(b)[,1,drop=FALSE]
vunit_measure(b, x=s, unit=2, time=1, params=p)

## End(Not run)

spatPomp documentation built on Aug. 10, 2023, 1:10 a.m.