# 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 |

### 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)
``` |