plot_effect: Genome-wide effects plot

In SingerLab/gac: Genetic Analysis of Cells

Genome-wide effects plot

Description

Plot estimated effects genome wide.

Usage

plot_effect(
  cnr,
  effect.column = grep("Estimate", names(cnr$Y), value = TRUE)[1],
  type = "h",
  ylab = expression(hat(italic(Beta))),
  xlab = "HSA Genome",
  ...
)

Arguments

cnr	a cnr bundle
effect.column	character, name of the effect column to plot
type	character, plot type, defaults to "h". See plot
ylab	character, y-axis label
xlab	character, x-axis lable
...	additional arguments passed to plot

Value

base graphics histogram like plot containing the estimated genome-side genotypic effects.

Examples

data(cnr)

cnr <- histo_logit(cnr, trait = "binary1",
   pheno0 = 0, pheno1 = 1)

plot_effect(cnr, effect.column = "0.vs.1.lr.Estimate")  ## effect from logistic reg

SingerLab/gac documentation built on March 23, 2024, 5:15 a.m.