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
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
\name{SeqExpressionSet-class}
\docType{class}
\alias{SeqExpressionSet-class}
\alias{offst}
\alias{offst<-}
\alias{normCounts}
\alias{normCounts<-}
\alias{initialize,SeqExpressionSet-method}
\alias{exprs,SeqExpressionSet-method}
\alias{exprs<-,SeqExpressionSet,ANY-method}
\alias{offst,SeqExpressionSet-method}
\alias{offst<-,SeqExpressionSet,ANY-method}
\alias{offst<-,SeqExpressionSet-method}
\alias{counts,SeqExpressionSet-method}
\alias{counts<-,SeqExpressionSet-method}
\alias{normCounts,SeqExpressionSet-method}
\alias{normCounts<-,SeqExpressionSet,ANY-method}
\alias{normCounts<-,SeqExpressionSet-method}
\title{"SeqExpressionSet" class for collections of short reads}
\description{
This class represents a collection of digital expression data (usually counts from RNA-Seq technology) along with sample information.
}
\section{Objects from the Class}{
Objects of this class can be created from a call to the
\code{\link{newSeqExpressionSet}} constructor.
}
\section{Extends}{
Class \code{eSet}, directly.
Class \code{VersionedBiobase}, by class \code{eSet}, distance 2.
Class \code{Versioned}, by class \code{eSet}, distance 3.
}
\section{Slots}{
Inherited from \code{eSet}:
\describe{
\item{\code{assayData}}{Contains matrices with equal dimensions, and with
column number equal to \code{nrow(phenoData)}.\code{assayData} must
contain a matrix \code{counts} with rows represening features
(e.g., genes) and columns representing samples.
The optional matrices \code{normalizedCounts} and \code{offset} can be added to represent a normalization in terms of pseudo-counts or offset, respectively, to be used for subsequent analyses. See the vignette for details.
Class: \code{\link{AssayData-class}}.}
\item{\code{phenoData}}{Sample information. For compatibility with DESeq, there should be at least the column \code{conditions}. See \code{\link{eSet}} for details.}
\item{\code{featureData}}{Feature information. It is recomended to include at least length and GC-content information. This slot is used for \code{\link{withinLaneNormalization}}. See \code{\link{eSet}} for details.}
\item{\code{experimentData}}{See \code{\link{eSet}}}
\item{\code{annotation}}{See \code{\link{eSet}}}
\item{\code{protocolData}}{See \code{link{eSet}}}
}
}
\section{Methods}{
See \code{\link{eSet}} for inherited methods. Additional methods:
\describe{
\item{counts}{\code{signature(object="SeqExpressionSet")}: returns the \code{counts} matrix. }
\item{counts<-}{\code{signature(object = "SeqExpressionSet")}: method to replace the \code{counts} matrix.}
\item{normCounts}{\code{signature(object="SeqExpressionSet")}: returns the \code{normalizedCounts} matrix. }
\item{normCounts<-}{\code{signature(object = "SeqExpressionSet")}: method to replace the \code{normalizedCounts} matrix.}
\item{offst}{\code{signature(object = "SeqExpressionSet")}: returns the \code{offset} matrix.}
\item{offst<-}{\code{signature(object = "SeqExpressionSet")}: method to replace the \code{offset} slot.}
\item{boxplot}{\code{signature(x = "SeqExpressionSet")}: produces a boxplot of the log counts.}
\item{meanVarPlot}{\code{signature(x = "SeqExpressionSet")}: produces a \code{\link{smoothScatter}} plot of the mean variance relation. See \code{\link{meanVarPlot}} for details.}
\item{biasPlot}{\code{signature(x = "SeqExpressionSet", y = "character")}: produces a plot of the \code{\link{lowess}} regression of the counts on some covariate of interest (usually GC-content or length). See \code{\link{biasPlot}} for details.}
\item{wihtinLaneNormalization}{\code{signature(x = "SeqExpressionSet", y = "missing")}: within lane normalization for GC-content (or other lane specific) bias. See \code{\link{withinLaneNormalization}} for details.}
\item{betweenLaneNormalization}{\code{signature(x = "SeqExpressionSet")}: between lane normalization for sequencing depth and possibly other distributional differences between lanes. See \code{\link{betweenLaneNormalization}} for details.}
}
}
\author{Davide Risso <risso.davide@gmail.com>}
\seealso{
\code{\linkS4class{eSet}}, \code{\link{newSeqExpressionSet}}, \code{\link{biasPlot}}, \code{\link{withinLaneNormalization}}, \code{\link{betweenLaneNormalization}}
}
\examples{
showMethods(class="SeqExpressionSet", where=getNamespace("EDASeq"))
counts <- matrix(data=0, nrow=100, ncol=4)
for(i in 1:4) {
counts[,i] <- rpois(100,lambda=50)
}
cond <- c(rep("A", 2), rep("B", 2))
data <- newSeqExpressionSet(counts, phenoData=AnnotatedDataFrame(data.frame(conditions=cond)))
head(counts(data))
boxplot(data, col=as.numeric(pData(data)[,1])+1)
}
\keyword{classes}
|
