inst/auxiliary/makeChlens.R

chdata <- read.table("chlens.txt", sep="\t") # chromnumber, length
chnames <- paste("chr", c(1:22, "X", "Y"), sep="")
chdata$Chrom <- factor(chnames[chdata$V1], levels = chnames)
chlens <- as.vector(tapply(chdata$V2, list(chdata$Chrom), max))
names(chlens) <- chnames

save(chlens, file = "../../Data/sysdata.rda")

Try the CloneSeeker package in your browser

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

CloneSeeker documentation built on July 1, 2022, 3 a.m.