R/gtfToBed.R

In fischuu/GenomicTools: Collection of Tools for Genomic Data Analysis

gtfToBed <- function(gtf, correctBases=TRUE){
  correction <- 0
  if(correctBases) correction <- 1
  if(sum(class(gtf)=="gtf")==0) stop("This function requires a gtf-object as import (as given by the importGTF() function)")
 bedOut <- data.frame(Chr=gtf[,1],
                      Start=gtf[,4] - correction,
                      End=gtf[,5] - correction,
                      Gene=gtf$gene_id,
                      stringsAsFactors=FALSE)

 colnames(bedOut) <- c("Chr", "Start", "End", "Gene")
 class(bedOut) <- append(class(bedOut), "bed")
 bedOut
}