plot_effect: Genome-wide effects plot

View source: R/logistic_reg_histology.R

plot_effectR Documentation

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 Oct. 14, 2022, 4:09 a.m.