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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/splitSE.R
\name{splitSE}
\alias{splitSE}
\title{Split SummarizedExperiment by gene categories}
\usage{
splitSE(se, splitDf, assayName)
}
\arguments{
\item{se}{A \code{SummarizedExperiment} object.}
\item{splitDf}{A \code{data.frame} with feature IDs. Each column
represents a separate feature type, and the features in a given
row are considered representatives of the same feature (and will
be represented as one feature in the output object).}
\item{assayName}{A character scalar, indicating the assay of \code{se}
that will be split. Must be one of \code{assayNames(se)}.}
}
\value{
A \code{SummarizedExperiment} object with the same columns
as the input object, and the same number of assays as the number
of columns in \code{splitDf}. The assays will be named by the
column names of \code{splitDf}. The \code{colData} and \code{metadata}
of the input \code{SummarizedExperiment} object are copied to the
output object. The row names are set to the feature IDs in the
first column of \code{splitDf}.
}
\description{
Construct a new SummarizedExperiment by splitting one of the assays
into a list of assays, each of which contains features of a given
'type'. It is assumed that there is a one-to-one correspondence
between feature sets of different types; for example, these can be
spliced and unspliced variants of the same transcripts. The type of
each feature in the original SummarizedExperiment, and the correspondence
between the features of different types, are given in a \code{data.frame}.
}
\examples{
se <- SummarizedExperiment::SummarizedExperiment(
assays = S4Vectors::SimpleList(
counts = as(matrix(1:15, nrow = 5), "sparseMatrix"),
logcounts = log2(matrix(1:15, nrow = 5))
),
colData = S4Vectors::DataFrame(sID = paste0("S", 1:3),
condition = c("A", "A", "B")),
metadata = list(md1 = "annotation")
)
rownames(se) <- paste0("G", 1:5)
colnames(se) <- paste0("P", 1:3)
splitDf <- data.frame(spliced = c("G1", "G2", "G6"),
unspliced = c("G3", "G5", "G4"),
stringsAsFactors = FALSE)
splse <- splitSE(se = se, splitDf = splitDf, assayName = "counts")
}
|
