Nothing
##' Vector of simulated measurements
##'
##' Simulate from the unit measurement model density function for each unit
##' @param object a \code{spatPomp} object
##' @param x numeric; state at which measurements are simulated
##' @param times numeric; time at which measurements are simulated
##' @param params numeric; parameter set at which measurements are simulated
##' @param \dots additional parameters will be ignored
##' @name vec_rmeasure
##' @include spatPomp_class.R
##' @author Kidus Asfaw
##' @rdname vec_rmeasure
##' @return An array of dimension \code{length(unit_names(object))} by \code{dim(x)[2]} by \code{dim(x)[3]}
##' representing each unit's simulated measurement assessed for each replicate in \code{x} for each observation time.
NULL
setGeneric("vec_rmeasure", function(object,...)standardGeneric("vec_rmeasure"))
##' @name vec_rmeasure-spatPomp
##' @rdname vec_rmeasure
##' @aliases vec_rmeasure,spatPomp-method
##' @export
setMethod(
"vec_rmeasure",
signature=signature(object="spatPomp"),
definition=function (object, x, times, params, ...){
vec_rmeasure.internal(object=object,x=x,times=times,params=params,...)
}
)
vec_rmeasure.internal <- function (object, x, times, params, .gnsi = TRUE, ...) {
pompLoad(object)
nunits <- length(unit_names(object))
nparticles <- ncol(x)
ntimes <- length(times)
storage.mode(x) <- "double"
storage.mode(params) <- "double"
sims <- array(dim=c(nunits,nparticles,ntimes))
for(i in 1:nunits){
# for girf params is not nparams by nparticles. instead it's npars by nunits by nparticles
if(length(dim(params)) > 2){
params <- params[,i,]
}
sims[i,,] <- .Call(do_runit_measure,object,x,times,i,params,.gnsi)
}
pompUnload(object)
return(sims)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.