In kbroman/qtlcharts: Interactive Graphics for QTL Experiments

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

set.seed(69891250)
library(qtl)
library(qtlcharts)
p <- 500
n <- 300
SD <- runif(p, 1, 5)
r <- runif(p, -1, 1)
scat2 <- vector("list", p)
for(i in 1:p)
   scat2[[i]] <- matrix(rnorm(2*n), ncol=2) %*% chol(SD[i]^2*matrix(c(1, r[i], r[i], 1), ncol=2))
scat1 <- cbind(SD=SD, r=r)

scat2scat(scat1, scat2, chartOpts=list(width=1000, height=500))

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