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

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=32, method="hk")
iplotScanone(out.hk, grav, lodcolumn=1, pheno.col=32)


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