geplot: Plots 95% credibility interval of GxE interaction

In BayesTwin: Bayesian Analysis of Item-Level Twin Data

Description

This function plots the 95% credibility region of a GxE interaction effect.

Usage

geplot(var_a, samples_beta0, samples_beta1, main, xlab, ylab, col, ...)

Arguments

var_a	A single value representing the posterior point estimate for genetic variance
samples_beta0	A vector representing draws from the posterior distribution of beta0, as produced by the main function IRT_twin of this package
samples_beta1	A vector representing draws from the posterior distribution of beta1, as produced by the main function IRT_twin of this package
main	An overall title for the plot
xlab	A title for the x axis
ylab	A title for the y axis
col	Color of the lines
...	Further arguments for the default S3 plot method

Details

The 95 % credibility region for the GxE interaction effect is displayed for the entire range of estimated genotypic values in the interval of -2 and 2 standard deviations of the posterior point estimate of genetic variance.

The function expects as input the posterior point estimate for genetic variance and two vectors representing draws from the posterior distributions of beta0 and beta1 respectively, as produced by the main function IRT_twin of this package.

Value

Plot of the 95% credibility region of the GxE interaction.

Author(s)

Inga Schwabe

Examples

#Using the output to plot the 95% credibility region 
data(results)
geplot(results$var_a, results$samples_beta0, results$samples_beta1)

Example output

Loading required package: foreign
Loading required package: coda
Loading required package: matrixStats
Loading required package: rjags
Linked to JAGS 4.3.0
Loaded modules: basemod,bugs

BayesTwin documentation built on May 2, 2019, 8:05 a.m.