rprocess-pomp {pomp}R Documentation

Simulate the process model of a partially-observed Markov process

Description

The method rprocess runs simulations of the the process-model portion of partially-observed Markov process.

Usage

rprocess(object, xstart, times, params, ...)
## S4 method for signature 'pomp':
rprocess(object, xstart, times, params, ...)

Arguments

object an object of class pomp.
xstart a rank-2 array containing the starting state of the system. Columns of xstart correspond to states; rows to state variables. If there is more than one column of xstart, multiple independent simulations will be performed, one corresponding to each column. Note that in this case, params must have the same number of columns as xstart.
times a numerical vector containing times. The first entry of times is the initial time (corresponding to xstart). Subsequent times are the additional times at which the state of the simulated processes are required.
params a rank-2 array of parameters with the parameters corresponding to the columns of xstart.
... at present, these are ignored.

Details

This function is essentially a wrapper around the user-supplied rprocess slot of the pomp object. For specifications on writing such a function, see pomp.

Value

rprocess returns a rank-3 array with rownames. Suppose x is the array returned. Then dim(x)=c(nvars,nreps,ntimes), where nvars is the number of state variables (=nrow(xstart)), nreps is the number of independent realizations simulated (=ncol(xstart)), and ntimes is the length of the vector times. x[,j,k] is the value of the state process in the j-th realization at time times[k]. In particular, x[,,1] must be identical to xstart. The rownames of x must correspond to those of xstart.

Author(s)

Aaron A. King (kingaa at umich dot edu)

See Also

pomp-class, pomp


[Package pomp version 0.22-6 Index]