scrna5/6 Jupyter Notebook lamindata

Train a machine learning model on a collection#

Here, we iterate over the artifacts within a collection to train a machine learning model at scale.

import lamindb as ln
import anndata as ad
import numpy as np

💡 connected lamindb: testuser1/test-scrna

ln.transform.stem_uid = "Qr1kIHvK506r"
ln.transform.version = "1"
ln.track()

💡 notebook imports: anndata==0.9.2 lamindb==0.69.2 numpy==1.26.4 torch==2.2.1

💡 saved: Transform(uid='Qr1kIHvK506r5zKv', name='Train a machine learning model on a collection', key='scrna5', version='1', type=notebook, updated_at=2024-03-26 12:03:31 UTC, created_by_id=1)

💡 saved: Run(uid='HMVWOW0j5fjiiCanMw55', transform_id=5, created_by_id=1)

Query our collection:

collection = ln.Collection.filter(
    name="My versioned scRNA-seq collection", version="2"
).one()
collection.describe()
Hide code cell output 
Collection(uid='IXzSUQjHF6tJklaahF5o', name='My versioned scRNA-seq collection', version='2', hash='HNR3VFV60_yqRnUka11E', visibility=1, updated_at=2024-03-26 12:03:09 UTC)

Provenance:
  📔 transform: Transform(uid='ManDYgmftZ8C5zKv', name='Standardize and append a batch of data', key='scrna2', version='1', type='notebook', updated_at=2024-03-26 12:02:46 UTC, created_by_id=1)
  👣 run: Run(uid='0SDPuvwMMvHbz8XZ6Kxl', started_at=2024-03-26 12:02:46 UTC, is_consecutive=True, transform_id=2, created_by_id=1)
  👤 created_by: User(uid='DzTjkKse', handle='testuser1', name='Test User1', updated_at=2024-03-26 12:00:16 UTC)
  ⬇️ input_of (core.Run): ['2024-03-26 12:03:21 UTC']
Features:
  var: FeatureSet(uid='JEqpMZS8zoJ3MS1sRj3e', n=36508, type='number', registry='bionty.Gene', hash='b5NMddLHEyZqn-vSYvBI', updated_at=2024-03-26 12:03:08 UTC, created_by_id=1)
    'MIR1302-2HG', 'FAM138A', 'OR4F5', 'None', 'None', 'None', 'None', 'None', 'None', 'None', 'OR4F29', 'None', 'OR4F16', 'None', 'LINC01409', 'FAM87B', 'LINC01128', 'LINC00115', 'FAM41C', 'None', ...
  obs: FeatureSet(uid='c5eCX4WUUMsPX69Enzc4', n=4, registry='core.Feature', hash='taDIZs0vy7CqVljgeqVZ', updated_at=2024-03-26 12:02:39 UTC, created_by_id=1)
    🔗 donor (12, core.ULabel): 'D496', '621B', 'A29', 'A36', 'A35', '637C', 'A52', 'A37', 'D503', '640C', ...
    🔗 tissue (17, bionty.Tissue): 'blood', 'thoracic lymph node', 'spleen', 'lung', 'mesenteric lymph node', 'lamina propria', 'liver', 'jejunal epithelium', 'omentum', 'bone marrow', ...
    🔗 cell_type (40, bionty.CellType): 'dendritic cell', 'effector memory CD4-positive, alpha-beta T cell, terminally differentiated', 'cytotoxic T cell', 'CD8-positive, CD25-positive, alpha-beta regulatory T cell', 'CD14-positive, CD16-negative classical monocyte', 'CD38-positive naive B cell', 'B cell, CD19-positive', 'CD4-positive, alpha-beta T cell', 'classical monocyte', 'T follicular helper cell', ...
    🔗 assay (3, bionty.ExperimentalFactor): '10x 3' v3', '10x 5' v2', '10x 5' v1'
Labels:
  🏷️ tissues (17, bionty.Tissue): 'blood', 'thoracic lymph node', 'spleen', 'lung', 'mesenteric lymph node', 'lamina propria', 'liver', 'jejunal epithelium', 'omentum', 'bone marrow', ...
  🏷️ cell_types (40, bionty.CellType): 'dendritic cell', 'effector memory CD4-positive, alpha-beta T cell, terminally differentiated', 'cytotoxic T cell', 'CD8-positive, CD25-positive, alpha-beta regulatory T cell', 'CD14-positive, CD16-negative classical monocyte', 'CD38-positive naive B cell', 'B cell, CD19-positive', 'CD4-positive, alpha-beta T cell', 'classical monocyte', 'T follicular helper cell', ...
  🏷️ experimental_factors (3, bionty.ExperimentalFactor): '10x 3' v3', '10x 5' v2', '10x 5' v1'
  🏷️ ulabels (12, core.ULabel): 'D496', '621B', 'A29', 'A36', 'A35', '637C', 'A52', 'A37', 'D503', '640C', ...

Create a map-style dataset#

Let us create a map-style dataset using using mapped(): a MappedCollection. This is what, for example, the PyTorch DataLoader expects as an input.

Under-the-hood, it performs a virtual inner join of the features of the underlying AnnData objects and thus allows to work with very large collections.

You can either perform a virtual inner join:

with collection.mapped(label_keys=["cell_type"], join="inner") as dataset:
    print(len(dataset.var_joint))

749

Or a virtual outer join:

dataset = collection.mapped(label_keys=["cell_type"], join="outer")

len(dataset.var_joint)

36508

This is compatible with a PyTorch DataLoader because it implements __getitem__ over a list of backed AnnData objects. The 5th cell in the collection can be accessed like:

dataset[5]
Hide code cell output 
{'x': array([ 0.   ,  0.   ,  0.   , ...,  0.   ,  0.   , -0.456], dtype=float32),
 '_storage_idx': 0,
 'cell_type': 28}

The labels are encoded into integers:

dataset.encoders
Hide code cell output 
{'cell_type': {'progenitor cell': 0,
  'conventional dendritic cell': 1,
  'mucosal invariant T cell': 2,
  'regulatory T cell': 3,
  'plasmacytoid dendritic cell': 4,
  'naive thymus-derived CD8-positive, alpha-beta T cell': 5,
  'classical monocyte': 6,
  'gamma-delta T cell': 7,
  'dendritic cell, human': 8,
  'macrophage': 9,
  'mast cell': 10,
  'naive thymus-derived CD4-positive, alpha-beta T cell': 11,
  'non-classical monocyte': 12,
  'effector memory CD8-positive, alpha-beta T cell, terminally differentiated': 13,
  'CD16-negative, CD56-bright natural killer cell, human': 14,
  'CD14-positive, CD16-negative classical monocyte': 15,
  'T follicular helper cell': 16,
  'alveolar macrophage': 17,
  'alpha-beta T cell': 18,
  'CD4-positive, alpha-beta T cell': 19,
  'megakaryocyte': 20,
  'plasmablast': 21,
  'lymphocyte': 22,
  'memory B cell': 23,
  'germinal center B cell': 24,
  'naive B cell': 25,
  'group 3 innate lymphoid cell': 26,
  'CD8-positive, alpha-beta memory T cell, CD45RO-positive': 27,
  'cytotoxic T cell': 28,
  'dendritic cell': 29,
  'CD38-positive naive B cell': 30,
  'plasma cell': 31,
  'B cell, CD19-positive': 32,
  'CD4-positive helper T cell': 33,
  'animal cell': 34,
  'effector memory CD4-positive, alpha-beta T cell, terminally differentiated': 35,
  'CD8-positive, CD25-positive, alpha-beta regulatory T cell': 36,
  'effector memory CD4-positive, alpha-beta T cell': 37,
  'CD16-positive, CD56-dim natural killer cell, human': 38,
  'CD8-positive, alpha-beta memory T cell': 39}}

Create a pytorch DataLoader#

Let us use a weighted sampler:

from torch.utils.data import DataLoader, WeightedRandomSampler

# label_key for weight doesn't have to be in labels on init
sampler = WeightedRandomSampler(
    weights=dataset.get_label_weights("cell_type"), num_samples=len(dataset)
)
dataloader = DataLoader(dataset, batch_size=128, sampler=sampler)

We can now iterate through the data loader:

for batch in dataloader:
    pass

Close the connections in MappedCollection:

dataset.close()
In practice, use a context manager
with collection.mapped(label_keys=["cell_type"]) as dataset:
    sampler = WeightedRandomSampler(
        weights=dataset.get_label_weights("cell_type"), num_samples=len(dataset)
    )
    dataloader = DataLoader(dataset, batch_size=128, sampler=sampler)
    for batch in dataloader:
        pass