diff options
Diffstat (limited to 'topics')
| -rw-r--r-- | topics/lmms/rqtl2/using-rqtl2-lmdb-adapter.gmi | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/topics/lmms/rqtl2/using-rqtl2-lmdb-adapter.gmi b/topics/lmms/rqtl2/using-rqtl2-lmdb-adapter.gmi index 2fa0491..bea482d 100644 --- a/topics/lmms/rqtl2/using-rqtl2-lmdb-adapter.gmi +++ b/topics/lmms/rqtl2/using-rqtl2-lmdb-adapter.gmi @@ -60,7 +60,24 @@ guix shell r r-thor r-rjson r-qtl2 -- \ # Example Rscript https://github.com/genenetwork/genenetwork3/pull/235/files ./lmdb_path + + ``` +### using this with rqtl2 Example + +```r +cross <- read_lmdb_cross(LMDB_DB_PATH) +summary(cross) +cat("Is this cross okay", check_cross2(cross), "\n") +warnings() # enable warnings for the debug purposes only! +pr <- calc_genoprob(cross) +out <- scan1(pr, cross$pheno, cores=4) +par(mar=c(5.1, 4.1, 1.1, 1.1)) +ymx <- maxlod(out) +plot(out, cross$gmap, lodcolumn=1, col="slateblue") # test generating of qtl plots + +``` + ### References => https://kbroman.org/qtl2/assets/vignettes/developer_guide.html |