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
42
|
\name{catchSalmon}
\alias{catchSalmon}
\alias{catchKallisto}
\title{Process Kallisto or Salmon Output}
\description{Read transcriptwise counts from kallisto or Salmon output for a series of samples and use the bootstrap samples to estimate the mapping uncertainty for each transcript.}
\usage{
catchKallisto(paths, verbose = TRUE)
catchSalmon(paths, verbose = TRUE)
}
\arguments{
\item{paths}{character vector giving paths to the directories created by kallisto.}
\item{verbose}{logical. If \code{TRUE}, summary information is printed as each sample is catch.}
}
\details{
These functions assume that kallisto or Salmon have been run to obtain estimated transcript counts for one or more RNA samples, and that bootstrap samples have also been generated.
These functions catch the counts and use the bootstrap samples to estimate an over-dispersion coefficient for each transcript.
Transcripts that overlap other transcripts and have greater read mapping uncertaintly will have greater over-dispersion coefficients.
The data is then ready for analysis in edgeR.
}
\value{
A list containing components
\item{counts}{matrix of transcript counts}
\item{annotation}{data.frame of transcript information.}
}
\author{Gordon Smyth}
\examples{
\dontrun{
s <- catchSalmon(paths)
dge <- DGEList(counts=s$counts/s$annotation$Overdispersion, genes=s$annotation)
}}
\concept{Transcript expression}
|
