if(require(brassicaData)){
data(raw_napus, package = "brassicaData", envir = environment())
dat <- intens_theta(raw_napus)
dat <- remove_suffix(dat, "_Grn")
dat <- geno_baf_rratio(dat, delthresh = 11)
dat <- segm(dat)
dat <- cnv(dat, dup = 0.03, del = -0.06)
data("synteny_blocks", package = "brassicaData", envir = environment())
dat <- trans_location(dat, synteny_blocks, min1 = 5, min2 = 5, maxdiff = 20)
plot_trans_locations(dat, sb = synteny_blocks$blocks)
}
Run the code above in your browser using DataLab