get.Phi: Compute the precision matrix Phi of observed nodes on a...

In rwc: Random Walk Covariance Models

Compute the precision matrix Phi of observed nodes on a graph.

Description

Given a Gaussian Markov random field defined by a precision matrix Q, this returns Phi, which is the precision matrix of the nodes indexed by 'obs.idx', computed using the Schur complement.

Usage

get.Phi(Q, obs.idx)

Arguments

Q	A precision matrix
obs.idx	A vector of unique indices of observed nodes in the graph defined by Q.

Value

A square matrix of dimension equal to the length of 'obs.idx'.

Author(s)

Ephraim M. Hanks

References

Hanks and Hooten 2013. Circuit theory and model-based inference for landscape connectivity. Journal of the American Statistical Association. 108(501), 22-33.

Examples

int=raster(nrow=30,ncol=30)
values(int)=1
TL.int=get.TL(int)
Q.int=get.Q(TL.int,1)
## get precision matrix of only the first 5 nodes
get.Phi(Q.int,1:5)

rwc documentation built on April 4, 2025, 4:40 a.m.