In kbroman/qtlcharts: Interactive Graphics for QTL Experiments

This document serves to test the R/qtlcharts function ipleiotropy.

knitr::opts_chunk$set(fig.width=9, fig.height=6)
set.seed(69891250)
library(qtl)
library(qtlcharts)
data(grav)

grav <- calc.genoprob(grav, step=1)
out.hk <- scanone(grav, pheno.col=c(156, 241), method="hk")

ipleiotropy(grav, out.hk, lodcolumn=1:2, pheno.col=c(156, 241), chr=1)

ipleiotropy(grav, pheno.col=c(156, 241), chr=1)

ipleiotropy(grav, out.hk, lodcolumn=1:2, pheno.col=c(156, 241), chr=1,
            interval=c(53,68))

kbroman/qtlcharts documentation built on May 10, 2023, 6:07 p.m.