tests/test_mapqtl_io.R

# test input/output in mapqtl format

library(qtl)
data(fake.4way)

write.cross(fake.4way, "mapqtl", "fake_4way_mapqtl")

x <- read.cross("mapqtl", "", genfile="fake_4way_mapqtl.loc",
                phefile="fake_4way_mapqtl.qua",
                mapfile="fake_4way_mapqtl_female.map")

x <- replace.map(x, pull.map(fake.4way))

comparecrosses(x, fake.4way)

Try the qtl package in your browser

Any scripts or data that you put into this service are public.

qtl documentation built on Nov. 28, 2023, 1:09 a.m.