R/readSNPTable.R

Defines functions readSNPTable

Documented in readSNPTable

readSNPTable <-
function(FileName, seqLengths, format="text"){

  message("reading SNP table")
  
  if(format=="GRanges"){
    snp <- readRDS(FileName)
  }
  else if(format=="text"){
    dat <- scan(FileName, what=list(chr=character(0), pos=integer(0)), sep="\t")
    snp <- GRanges(seqnames=dat$chr, ranges=IRanges(start=dat$pos, width=1), strand='*', seqlengths=seqLengths)
  }
  else{
    stop("unknown format")
  }

  # order by chromosome and CpG position
  snp <- snp[order(as.vector(seqnames(snp)), start(snp))]
  snp
  
}

Try the MethylSeekR package in your browser

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

MethylSeekR documentation built on Nov. 8, 2020, 6:57 p.m.