ACEt (version 1.5.3)

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

Description

Plot variance curves or a heritability curve (with 95

Usage

plot_acet(acet, boot = FALSE, heri = FALSE)

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.

boot

An logical indicator of whether the confidence bands estimated from the bootstrap method are plotted. The default is FALSE. Only available for the AtCtEt model.

heri

A logical indicator of whether to plot the dynamic heritability curve (currently only available for the B-splines). The default is FALSE.

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 <- AtCtEtp(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, mod=c('d','c','c'), knot_a = 5)
# plot_acet(result)