tests/test_tidyIO.R

library(qtl)
data(hyper)

# write to tidy format
write.cross(hyper, "tidy", "hyper_tidy")

# read back in
x <- read.cross("tidy", "", genfile="hyper_tidy_gen.csv",
                mapfile="hyper_tidy_map.csv", phefile="hyper_tidy_phe.csv",
                genotypes=c("BB", "BA", "AA"))

# compare results
comparecrosses(x, hyper)
kbroman/qtl documentation built on Aug. 26, 2024, 12:33 a.m.