In kbroman/qtlcharts: Interactive Graphics for QTL Experiments

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

set.seed(69891250)
library(qtl)
library(qtlcharts)
data(grav)
grav$pheno <- grav$pheno[,seq(1, nphe(grav), by=5)]

rowmedian <- apply(grav$pheno, 1, median, na.rm=TRUE)
grp <- as.numeric(rowmedian < median(unlist(grav$pheno), na.rm=TRUE))+1
iplot(grav$pheno$T0, grav$pheno$T60,
      group=grp,
      indID=paste0("RIL", 1:nind(grav)),
      chartOpts=list(xlab="Angle at time 0", ylab="Angle at 60 min"))

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