forked from lima1/PureCN
-
Notifications
You must be signed in to change notification settings - Fork 0
/
readCoverageFile.Rd
41 lines (36 loc) · 1.16 KB
/
readCoverageFile.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/readCoverageFile.R
\name{readCoverageFile}
\alias{readCoverageFile}
\title{Read coverage file}
\usage{
readCoverageFile(file, format, zero = NULL, read.length = 100)
}
\arguments{
\item{file}{Target coverage file.}
\item{format}{File format. If missing, derived from the file
extension. Currently GATK3 DepthofCoverage, GATK4 CollectFragmentCounts
(hdf5), and CNVkit formats supported.}
\item{zero}{Start position is 0-based. Default is \code{FALSE}
for GATK, \code{TRUE} for BED file based intervals.}
\item{read.length}{For output formats which do not provide both counts
and total coverages, approximate them using the specified read length.}
}
\value{
A \code{data.frame} with the parsed coverage information.
}
\description{
Read coverage file produced by external tools like The Genome Analysis
Toolkit or by \code{\link{calculateBamCoverageByInterval}}.
}
\examples{
tumor.coverage.file <- system.file("extdata", "example_tumor.txt.gz",
package = "PureCN")
coverage <- readCoverageFile(tumor.coverage.file)
}
\seealso{
\code{\link{calculateBamCoverageByInterval}}
}
\author{
Markus Riester
}