dsdive.obstx.matrix: Compute probability transition matrix for partially observed...

Description Usage Arguments Examples

View source: R/dsdive.obstx.matrix.R

Description

Given model parameters, this function will compute the probability transition matrix for a Continuous time Markov chain (CTMC) that is observed once, and then again at a time of tstep units of time later. If include.raw==TRUE, then raw components of the probability transition matrix will be returned so that the transition matrix can be computed for arbitrary timesteps.

Usage

1
2
3
4
5
6
7
8
9
dsdive.obstx.matrix(
  depth.bins,
  beta,
  lambda,
  s0,
  tstep,
  include.raw = FALSE,
  delta
)

Arguments

depth.bins

n x 2 Matrix that defines the depth bins. The first column defines the depth at the center of each depth bin, and the second column defines the half-width of each bin.

beta

the current depth bin transition model parameters

lambda

the current depth bin transition rate model parameters

s0

the stage for which to compute the transition matrix

tstep

Time between observations of the CTMC

include.raw

TRUE to include raw components of the probability transition matrix so that the matrix can be computed for arbitrary timesteps

delta

If delta>0, then the observation matrix and raw components computed will be for a transition matrix whose generator is perturbed to allow much faster computation.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
data('dive.sim')
attach(dive.sim)
attach(dive.sim$params)

o = dsdive.obstx.matrix(depth.bins = depth.bins, beta = beta, lambda = lambda, 
                        s0 = 1, tstep = 300, include.raw = TRUE, delta = 1e-10)


detach(dive.sim$params)
detach(dive.sim)

jmhewitt/dsdive documentation built on May 29, 2020, 5:18 p.m.