This document serves to test the R/qtlcharts function iplotCurves. This one has a plot with just curves.

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

Reduce size of data.

phe <- grav$pheno[,seq(1, nphe(grav), by=5)]
phe[14,ncol(phe)] <- NA
times <- as.numeric(sub("T", "", colnames(phe)))/60
iplotCurves(phe, times,
            chartOpts=list(curves_xlab="Time", curves_ylab="Root tip angle"))


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