DySeq (version 0.13)

EstTime: EstTime

Description

• Simulates k state-transition tables (see: StateTrans) based one state-transtion table containing expected population frequencies using the EstFreq function.
• Repeats simulations for a number of defined time intervals or until a termination criterion is reached.
• The estimated number of cells with low or zero frequencies are computed.
• The proportion of low/zero frequencies is plotted against the number of time points/intervalls.

Usage

EstTime(x, t = NA, crit = 0.05, zero = T, min.cell = 5, k = 20000, pos = "bottomleft", smoothed = T, show.it = 10, max.it = 10000)

Arguments

x

a matrix containing the assumed probabilities (see EstFreq for further detail!)

t

optional: a vector of time intervalls for which frequencies should be simulated. Can be used to decrease simulation time!

crit

optional: but must be sepcified if t is not; Simulations will end if the relative frequency of zero cells is less than crit

zero

optional: if FALSE simulations will end if number of low frequencies (instead of zero frequencies) is less than crit

min.cell

a single integer defines what counts as a low frequency (lower than 5 by convention)

k

Number of simulations (at least 20.000 is recommended)

pos

position of the output's legend. Options are: "bottomleft", "bottomright", "upperleft", and "upperright".

smoothed

logical value. If true, output lines will be smoothed!

show.it

single integer that defines which steps of iteration protocol should be shown. Only active if t is not defined else iteration protocoll is replaced with a progression bar

max.it

single integer that defines the maximum number of iterations if t is not specified.

Value

EstTime object; a list of three vectors if printed, provides a plot of expected number of low and zero cell frequencies

Details

First vector represents time points, second vector provides rel. frequency of cases with zero's, third vector rel. frequency of cases with low cell frequencies.

Examples

## Not run: 
# my.trans.table<-matrix(c(0.57, 0.13,0.05,0.05,0.05, 0.05,0.05,0.05),4,2)
# my.cellproblems<-EstTime(my.trans.table, k=500) # low k only for examplification purposes!
# my.cellproblems
# ## End(Not run)