iplotMScanone within R Markdown

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

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

Reduce size of data.

grav$pheno <- grav$pheno[,seq(1, nphe(grav), by=5)]
times <- as.numeric(sub("T", "", phenames(grav)))/60

Run scanone.

grav <- calc.genoprob(grav, step=1)
out.hk <- scanone(grav, pheno.col=1:nphe(grav), method="hk")
iplotMScanone(out.hk, grav, chartOpts=list(eff_ylab="QTL effect"))


Try the qtlcharts package in your browser

Any scripts or data that you put into this service are public.

qtlcharts documentation built on Jan. 8, 2022, 1:06 a.m.