rprocess runs simulations of the the process-model portion of partially-observed Markov process.
This function is part of the low-level interface to pomp objects.
This help page does not give instructions on the implementation of models: see pomp for instructions.rprocess(object, xstart, times, params, ...)
## S3 method for class 'pomp':
rprocess(object, xstart, times, params, offset, \dots)pomp.xstart correspond to states; rows to state variables.
If there is more than one column of xstart, multiple independent simulations will be performed, times is the initial time (corresponding to xstart).
Subsequent times are the additional times at which the state of the simulated processes are required.xstart.offset times.rprocess returns a rank-3 array with rownames.
Suppose x is the array returned.
Then dim(x)=c(nvars,nreps,ntimes-offset), 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+offset].
The rownames of x must correspond to those of xstart.rprocess slot of the pomp object.
For specifications on writing such a function, see pomp.pomp-class, pomp