inst/scripts/createGenomicRanges.R

library(GenomicRanges)
a <- read.csv2("~/../Desktop/Libro1.csv", as.is = TRUE)

colnames(a)[c(1, 4:6)] <- c("scoreInvHap.name", "Inv.freq", "Haplotypes", "Num.SNPs")
a[, "Inv.freq"] <- as.numeric(a[, "Inv.freq"])

inversionGR <- GRanges(gsub(",", "", a$Coordinates))
mcols(inversionGR) <- a[, -3]
names(inversionGR) <- a[, 1]
devtools::use_data(inversionGR)

Try the scoreInvHap package in your browser

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

scoreInvHap documentation built on Feb. 6, 2021, 2 a.m.