EstTime: EstTime

In DySeq: Functions for Dyadic Sequence Analyses

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.

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.

Value

EstTime object; a list of three vectors if printed, provides a plot of expected number of low and zero 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)

DySeq documentation built on May 1, 2019, 10:14 p.m.