@@ -14,13 +14,35 @@ library(tidySingleCellExperiment)
1414Load the data
1515
1616``` r
17- hca_metadata = get_metadata()
17+ get_metadata()
1818```
1919
20+ ## # Source: table<metadata> [?? x 48]
21+ ## # Database: sqlite 3.39.2 [/vast/projects/RCP/human_cell_atlas/metadata.sqlite]
22+ ## .cell .sample .samp…¹ assay assay…² cell_…³ cell_…⁴ devel…⁵ devel…⁶ disease
23+ ## <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr>
24+ ## 1 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
25+ ## 2 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
26+ ## 3 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… lumina… CL:000… 31-yea… HsapDv… normal
27+ ## 4 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… lumina… CL:000… 31-yea… HsapDv… normal
28+ ## 5 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
29+ ## 6 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
30+ ## 7 AAACCT… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
31+ ## 8 AAACGG… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
32+ ## 9 AAACGG… 5f20d7… D17PrP… 10x … EFO:00… lumina… CL:000… 31-yea… HsapDv… normal
33+ ## 10 AAACGG… 5f20d7… D17PrP… 10x … EFO:00… basal … CL:000… 31-yea… HsapDv… normal
34+ ## # … with more rows, 38 more variables: disease_ontology_term_id <chr>,
35+ ## # ethnicity <chr>, ethnicity_ontology_term_id <chr>, file_id <chr>,
36+ ## # is_primary_data.x <chr>, organism <chr>, organism_ontology_term_id <chr>,
37+ ## # sample_placeholder <chr>, sex <chr>, sex_ontology_term_id <chr>,
38+ ## # tissue <chr>, tissue_ontology_term_id <chr>, dataset_id <chr>,
39+ ## # collection_id <chr>, cell_count <int>, dataset_deployments <chr>,
40+ ## # is_primary_data.y <chr>, is_valid <int>, linked_genesets <int>, …
41+
2042Explore the HCA content
2143
2244``` r
23- hca_metadata | >
45+ get_metadata() | >
2446 distinct(tissue , file_id ) | >
2547 count(tissue ) | >
2648 arrange(desc(n ))
@@ -48,33 +70,35 @@ Query
4870``` r
4971options(" restore_SingleCellExperiment_show" = TRUE )
5072sce =
51- hca_metadata | >
73+ get_metadata() | >
5274 filter(
53- ethnicity == " African" &
75+ ethnicity == " African" &
5476 assay %LIKE % " %10x%" &
5577 tissue == " lung parenchyma" &
5678 cell_type %LIKE % " %CD4%"
5779 ) | >
80+
5881 get_SingleCellExperiment()
5982```
6083
61- ## Reading 9 files.
84+ ## Reading 1 files.
6285
63- ## .........
86+ ## .
6487
6588``` r
6689sce
6790```
6891
6992 ## class: SingleCellExperiment
70- ## dim: 60661 1571
71- ## metadata(63): X_normalization batch_condition ... schema_version title
93+ ## dim: 28024 1571
94+ ## metadata(0):
7295 ## assays(1): X
73- ## rownames(60661 ): TSPAN6 TNMD ... RP11-175I6.6 PRSS43P
96+ ## rownames(28024 ): TSPAN6 TNMD ... RP11-107E5.4 RP11-299P2.2
7497 ## rowData names(0):
75- ## colnames(1571):
76- ## ACAGCCGGTCCGTTAA_F02526_4fc10a6b85e5fa688b253db4e0db8ba0
77- ## GGGAATGAGCCCAGCT_F02526_4fc10a6b85e5fa688b253db4e0db8ba0 ...
78- ## TCGTAGAGTATATGGA_SC29_948efa032d715a169f9edae65c696836
79- ## GACGTTACACTGTGTA_SC29_948efa032d715a169f9edae65c696836
80- ## colData names(45): .sample .sample_name ... created_at.y updated_at.y
98+ ## colnames(1571): ACAGCCGGTCCGTTAA_F02526 GGGAATGAGCCCAGCT_F02526 ...
99+ ## TACAACGTCAGCATTG_SC84 CATTCGCTCAATACCG_F02526
100+ ## colData names(47): .sample .sample_name ... file_id_db
101+ ## stringr..str_remove.stringr..str_remove..cell...sample...._...
102+ ## reducedDimNames(0):
103+ ## mainExpName: NULL
104+ ## altExpNames(0):
0 commit comments