Simulated capture-mark-recapture (CMR) data for Bayesian Survival Trajectory Analysis (BaSTA).
data("bastaCMRdat")A data frame with 500 observations on the following 25 variables (see details).
This dataset was created by stochastically simulating a hypothetical population with different mortality patterns between males and females and with proportional decreases in mortality as a function of a hypothetical continuous covariate (e.g. birth weight). The population was simulated for 100 years, at each one of which 100 individuals were born. The number of females per generation was randomly drawn from a binomial distribution using function rbinom with probability of 0.5 (i.e. 1:1 sex ratio). The individual continuous covariate was randomly drawn from a random normal distribution (with function rnorm) with mean parameter equal to 0 (e.g. anomaly of weights) and standard deviation equal to 1. The time of death for each individual was inversed sampled from a Gompertz CDF of ages at death. The Gompertz parameters for females were: \(b_0 = -4\) and \(b_1 = 0.15\); and for males at \(b_0 = -3\) and \(b_1 = 0.15\). The gamma parameter for the continuous covariate was \(\gamma = 0.2\).
The study was assumed to start at year 51 and to finish at year 70. Recapture probability was set to 0.6 and thus each observation per individual was randomly drawn from a Bernoulli trial with parameter \(p = 0.6\). Captures at birth and recoveries were randomly drawn from a Bernoulli trial with parameters \(p = 0.5\) and \(p = 0.2\), respectively.
Therefore, the resulting dataset includes individuals where the data are left-truncated and/or right-censored. This is typical of capture mark recovery datasets.
Here is the description of each column:
id: ID for each individual.
birth: Integer times of birth or NA when unknown.
death: Integer times of death or NA when unknown.
Y.51 - Y.70: Columns of the recapture matrix
Sex: Character string indicating the sex for the individual record (i.e., “Female” and “Male”)
Weight: Simumated body weights at birth.
## Load data:
data("bastaCMRdat", package = "BaSTA")
## Check data consistency:
checkedData <- DataCheck(bastaCMRdat, dataType = "CMR", studyStart = 51,
studyEnd = 70)
Run the code above in your browser using DataLab