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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/rowSubset.R
\docType{methods}
\name{rowSubset}
\alias{rowSubset}
\alias{rowSubset<-}
\alias{rowSubset,SingleCellExperiment-method}
\alias{rowSubset<-,SingleCellExperiment-method}
\title{Get or set the row subset}
\usage{
\S4method{rowSubset}{SingleCellExperiment}(x, field = "subset", onAbsence = "none")
\S4method{rowSubset}{SingleCellExperiment}(x, field = "subset", ...) <- value
}
\arguments{
\item{x}{A \linkS4class{SingleCellExperiment} object.}
\item{field}{String containing the name of the field in the \code{\link{rowData}} to get or set subsetting data.}
\item{onAbsence}{String indicating an additional action to take when labels are absent:
nothing (\code{"none"}), a warning (\code{"warn"}) or an error (\code{"error"}).}
\item{...}{Additional arguments, currently ignored.}
\item{value}{Any character, logical or numeric vector specifying rows of \code{x} to include in the subset.
Alternatively \code{NULL}, in which case existing subsetting information is removed.}
}
\value{
For \code{rowSubset}, a logical vector is returned specifying the rows to retain in the subset of interest.
If no subset is available, a \code{NULL} is returned (and/or a warning or error, depending on \code{onAbsence}).
For \code{rowSubset<-}, a modified \code{x} is returned a subsetting vector in its \code{\link{rowData}}.
}
\description{
Get or set the row subset in an instance of a \linkS4class{SingleCellExperiment} class.
This is assumed to specify some interesting subset of genes to be favored in downstream analyses.
}
\details{
A frequent task in single-cell data analyses is to focus on a subset of genes of interest,
e.g., highly variable genes, derived marker genes for clusters, known markers for cell types.
A related task is to filter out uninteresting genes such as ribosomal protein genes or mitochondrial transcripts,
in which case we want to subset to exclude those genes.
These functions store a set of genes of interest inside a \linkS4class{SingleCellExperiment}
for later retrieval and use in downstream functions.
Character and numeric \code{value} are converted to logical vectors that are parallel to the rows of \code{x},
allowing them to be added to the \code{\link{rowData}} for synchronized row-level operations.
For developers, \code{onAbsence} is provided to make it easier to mandate that \code{x} actually has labels.
This avoids silent \code{NULL} values that flow to the rest of the function and make debugging difficult.
}
\examples{
example(SingleCellExperiment, echo=FALSE) # Using the class example
rowSubset(sce, "hvgs") <- 1:10
rowSubset(sce, "hvgs")
rowSubset(sce) <- rbinom(nrow(sce), 1, 0.5)==1
rowSubset(sce)
}
\seealso{
\linkS4class{SingleCellExperiment}, for the underlying class definition.
}
\author{
Aaron Lun
}
|
