ACEt (version 1.3)

plot_acet: Plot variance curves for various ACE(t) or ACE(t)-p models

Description

Plot variance curves of the A, C and E components with respect to age modelled by B-splines or P-splines.

Usage

plot_acet(acet)

Arguments

acet

An object obtained from the B-splines or P-splines functions. For the P-splines functions, an object from the MCMC method must be used.

References

Liang He, Mikko J. Sillanp��, Karri Silventoinen, Jaakko Kaprio, Janne Pitk�niemi, Estimating modifying effect of age on genetic and environmental variance components in twin models. Genetics, 2016

Examples

# data(data_ace)
# result <- AtCtEp(data_ace$mz, data_ace$dz, knot_a = 5, knot_c = 4)
# result_mc <- acetp_mcmc(result, iter_num=10000, burnin = 500)
# plot_acet(result_mc)

# result <- AtCtEt(data_ace$mz, data_ace$dz, model=c('d','c','c'), knot_a = 5)
# plot_acet(result)