Impute Blocks for Extra Stops in "epidata"
Objects
epidata_intersperse.Rd
This function modifies an object inheriting from class "epidata"
such
that it features the specified stop time points. For this purpose, the time
interval in the event history into which the new stop falls will be split
up into two parts, one block for the time period until the new stop – where
no infection or removal occurs – and the other block for the time period
from the new stop to the end of the original interval.
Main application is to enable the use of knots
in twinSIR
, which
are not existing stop time points in the "epidata"
object.
Arguments
- epidata
an object inheriting from class
"epidata"
.- stoptimes
a numeric vector of time points inside the observation period of the
epidata
.- verbose
logical indicating if a
txtProgressBar
should be shown while inserting blocks for extrastoptimes
.
See also
as.epidata.epidataCS
where this function is used.
Examples
data("hagelloch")
subset(hagelloch, start < 25 & stop > 25 & id %in% 9:13, select = 1:7)
# there is no "stop" time at 25, but we can add this extra stop
nrow(hagelloch)
moreStopsEpi <- intersperse(hagelloch, stoptimes = 25)
nrow(moreStopsEpi)
subset(moreStopsEpi, (stop == 25 | start == 25) & id %in% 9:13, select = 1:7)