How to subset SummarizedExperiment object with multiple experiments?
I have a big SummarizedExperiment object. I want to subset all of the tables based on the information in the colData. I think there should be an efficient way rather than piping filters like:
filter <- colData(dat)[id ,which(colData(dat)$parameter) == "cancer" ]
cancer <- assays(dat)[,filter]
Now, I need to do that for all the values of parameter. There should be a complete function somewhere to automate this process and be memory efficient.
Any ideas?
• 4,439 views
•
link
updated
by
Biostar
•
written
by
Omar Mohamed •
•
written
by
Omar Mohamed •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Problem after subset of single cell integrated dataset in seurat using sctransform
written by NGS_enthusiast •
-
I'm not sure what is going on with my DEXSeq object
written by cftattersfield •
-
Diffbind and DESeq2 comparison for differential sites
written by Ankit
-
Filter out doublet and Clustering
written by synat.keam •
-
SoupX with scDblFinder, matrix error
written by sushisena •
-
How to get the counts matrix correctly for DESeq2 from the counts table resulting from nextflow rna…
written by Josh •
-
Heatmap deseq2
written by bart •
-
Should I cluster and PCA the whole sample or only the genes of interest?
written by Omar Mohamed •
-
DESeq2 uses a lot of memory even for a few samples
written by bounlu
-
Deseq2 Single-Cell analysis with two variables
written by michele.tebaldi.92 •
dat[,dat$parameter == "cancer"]will subset the entire SE. You can then extract assays as you like. The idea of SE is exactly that, so not having the need to filter every assay but one filter call on the entire object subsets all assays, all colData etc. What isidhere? I do not see exactly how filtering for rows of colData makes sense.