EpiModel (version 0.95)

dissolution.coefs: Calcuate Dissolution Coefficients for Stochastic Network Models

Description

This function calculates dissolution coefficients from a duration to pass as terms to an ERGM/STERGM fit and applies a coefficient correction to account for deaths.

Usage

dissolution.coefs(dissolution, duration, d.rate = 0)

Arguments

dissolution

a right-hand sided STERGM dissolution formula (see epiNet.est); currently limited only to simple ~offset(edges) dissolution models.

duration

an edge duration in arbitrary time units.

d.rate

background death rate in the absence of disease.

Value

• A list of with the following elements:

Details

This function performs two duties prior to the estimation phase of epiNet class stochastic network models.

1. Transforming duration to coefficient form:an average duration of edges in a population must be mathematically transformed to logit coefficients. The theory and mathematics of the transformation are further developed in thenetworkDynamicpackage vignette.
2. Adjusting coefficients to account for death:In a dynamic network simulation in an open population (in which there are births and deaths), it is necessary to adjust the dissolution coefficient for the STERGM simulations to account for the death as a competing risk to edge dissolution.

Future releases of this software will allow for more flexibility in the possible dissolution models that may be calculated here, including models that have heterogenous dissolution probabilities conditional on nodal or edge attributes.

See Also

For a detailed explanation of this function, please see the HTML vignette: ../doc/epiNetUtils.html{EpiModel Network Utility Functions}.

Examples

dissolution <- ~offset(edges)
duration <- 25
dissolution.coefs(dissolution, duration)
dissolution.coefs(dissolution, duration, d.rate=0.001)