-
Notifications
You must be signed in to change notification settings - Fork 10
Closed
Description
Hi there,
I have an sce object obtained after doing CATALYST clustering of the flow cytometry data, and I want to use tidySingleCellExperiement package to manipulate the object.
I was trying to plot some channels/markers as below, but unsure which data the function .abundance_counts is drawing from, given the error below.
Would you mind giving me some pointer?
Thank you.
> F37_sce_backboneClustering
# A SingleCellExperiment-tibble abstraction: 4,955,453 × 17
# Features=270 | Assays=exprs
.cell sample_id cell_count_prePeacoQC cell_count_postPeacoQC manualInspect_strang…¹
<chr> <fct> <fct> <fct> <fct>
1 CATALYST28meta16_12_1 F37_10D_CD4… 50819 50785 Expected
2 CATALYST28meta16_3_2 F37_10D_CD4… 50819 50785 Expected
3 CATALYST28meta16_3_3 F37_10D_CD4… 50819 50785 Expected
4 CATALYST28meta16_2_4 F37_10D_CD4… 50819 50785 Expected
5 CATALYST28meta16_3_5 F37_10D_CD4… 50819 50785 Expected
6 CATALYST28meta16_2_6 F37_10D_CD4… 50819 50785 Expected
7 CATALYST28meta16_8_7 F37_10D_CD4… 50819 50785 Expected
8 CATALYST28meta16_3_8 F37_10D_CD4… 50819 50785 Expected
9 CATALYST28meta16_3_9 F37_10D_CD4… 50819 50785 Expected
10 CATALYST28meta16_3_10 F37_10D_CD4… 50819 50785 Expected
# ℹ 4,955,443 more rows
# ℹ abbreviated name: ¹manualInspect_strangeFCS
# ℹ 12 more variables: preFilter_liveCD3CD8_count <fct>, flowJo_prePeacoQC_gMFI_PE <fct>,
# flowCore_postPeacoQC_medMFI_PE <fct>, cluster_id <fct>, CATALYST28meta16 <fct>, nColID <fct>,
# CATALYST28meta16_nColID <chr>, CATALYST28meta16_uniqueLetterID <fct>, UMAP1 <dbl>, UMAP2 <dbl>,
# TSNE1 <dbl>, TSNE2 <dbl>
# ℹ Use `print(n = ...)` to see more rows
> F37_sce_backboneClustering |>
+ join_features(features=c("CD4", "CD8")) |>
+ ggplot(aes(CATALYST28meta16, .abundance_counts + 1, fill=CATALYST28meta16)) +
+ geom_boxplot(outlier.shape=NA) +
+ scale_y_log10() +
+ custom_theme
tidySingleCellExperiment says: This operation lead to duplicated cell names. A data frame is returned for independent data analysis.
Error in `geom_boxplot()`:
! Problem while computing aesthetics.
ℹ Error occurred in the 1st layer.
Caused by error in `FUN()`:
! object '.abundance_counts' not found
Backtrace:
1. base (local) `<fn>`(x)
2. ggplot2:::print.ggplot(x)
4. ggplot2:::ggplot_build.ggplot(x)
5. ggplot2:::by_layer(...)
12. ggplot2 (local) f(l = layers[[i]], d = data[[i]])
13. l$compute_aesthetics(d, plot)
14. ggplot2 (local) compute_aesthetics(..., self = self)
15. base::lapply(aesthetics, eval_tidy, data = data, env = env)
16. rlang (local) FUN(X[[i]], ...)
Metadata
Metadata
Assignees
Labels
No labels