Description Usage Arguments Details Value Author(s) Examples
View source: R/QFeatures-join.R
This function applies a full-join type of operation on 2 or more
assays in a QFeatures
instance.
1 | joinAssays(x, i, name = "joinedAssay")
|
x |
An instance of class QFeatures. |
i |
The indices or names of al least two assays to be joined. |
name |
A |
The rows to be joined are chosen based on the rownames of the respective assays. It is the user's responsability to make sure these are meaningful, such as for example refering to unique peptide sequences or proteins.
The join operation acts along the rows and expects the samples
(columns) of the assays to be disjoint, i.e. the assays mustn't
share any samples. Rows that aren't present in an assay are set to
NA
when merged.
The rowData
slots are also joined. However, only columns that
are shared and that have the same values for matching columns/rows
are retained. For example of a feature variable A
in sample S1
contains value a1
and variable A
in sample S2
in a different
assay contains a2
, then the feature variable A
is dropped in
the merged assay.
The joined assay is linked to its parent assays through an AssayLink
object. The link between the child assay and the parent assays is based on
the assay row names, just like the procedure for joining the parent assays.
A QFeatures
object with an additional assay.
Laurent Gatto
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 | ## -----------------------------------------------
## An example QFeatures with 3 assays to be joined
## -----------------------------------------------
data(feat2)
feat2
feat2 <- joinAssays(feat2, 1:3)
## Individual assays to be joined, each with 4 samples and a
## variable number of rows.
assay(feat2[[1]])
assay(feat2[[2]])
assay(feat2[[3]])
## The joined assay contains 14 rows (corresponding to the union
## of those in the initial assays) and 12 samples
assay(feat2[["joinedAssay"]])
## The individual rowData to be joined.
rowData(feat2[[1]])
rowData(feat2[[2]])
rowData(feat2[[3]])
## Only the 'Prot' variable is retained because it is shared among
## all assays and the values and coherent across samples (the
## value of 'Prot' for row 'j' is always 'Pj'). The variable 'y' is
## missing in 'assay1' and while variable 'x' is present is all
## assays, the values for the shared rows are different.
rowData(feat2[["joinedAssay"]])
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.