Plot variance curves for various ACE(t) or ACE(t)-p models with 95% CIs

Description

Plot variance curves or a heritability curve (with 95% CIs) of the A, C and E components with respect to age modelled by B-splines or P-splines.

Usage

1
plot_acet(acet, boot = FALSE, heri = FALSE, xlab, ylab, main, col, legend = TRUE)

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. The default is FALSE.

xlab

The 'xlab' argument in the plot funciton. The default is 'Age'.

ylab

The 'ylab' argument in the plot funciton. The default is 'Variance'.

main

The 'main' argument in the plot funciton. The default is 'Variance curves of the A, C, and E components' for variances and 'Dynamic heritability' for heritability.

col

The 'col' argument in the plot funciton.

legend

An logical indicator of whether the default legend is plotted.

Author(s)

Liang He

References

He, L., Sillanpää, M.J., Silventoinen, K., Kaprio, J. and Pitkäniemi, J., 2016. Estimating Modifying Effect of Age on Genetic and Environmental Variance Components in Twin Models. Genetics, 202(4), pp.1313-1328.

Examples

1
2
3
4
5
6
7
# 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 = 9)
# plot_acet(result)