Simulate twin data from a linear normal ACE/ADE/AE model.
twinsim(
nMZ = 100,
nDZ = nMZ,
b1 = c(),
b2 = c(),
mu = 0,
acde = c(1, 1, 0, 1),
randomslope = NULL,
threshold = 0,
cens = FALSE,
wide = FALSE,
...
)
Number of monozygotic twin pairs
Number of dizygotic twin pairs
Effect of covariates (labelled x1,x2,...) of type 1. One distinct covariate value for each twin/individual.
Effect of covariates (labelled g1,g2,...) of type 2. One covariate value for each twin pair.
Intercept parameter.
Variance of random effects (in the order A,C,D,E)
Logical indicating wether to include random slopes of the variance components w.r.t. x1,x2,...
Treshold used to define binary outcome y0
Logical variable indicating whether to censor outcome
Logical indicating if wide data format should be returned
Additional arguments parsed on to lower-level functions
Klaus K. Holst
twinlm