Skip to content

Commit

Permalink
[pre-commit.ci] auto fixes from pre-commit.com hooks
Browse files Browse the repository at this point in the history
for more information, see https://pre-commit.ci
  • Loading branch information
pre-commit-ci[bot] committed Mar 29, 2023
1 parent ce30a82 commit 8ad2e8b
Showing 1 changed file with 9 additions and 7 deletions.
16 changes: 9 additions & 7 deletions scirpy/datasets/_processing_scripts/maynard2020.py
Original file line number Diff line number Diff line change
Expand Up @@ -20,23 +20,23 @@
from pathlib import Path

import numpy as np
import numpy.testing as npt
import pandas as pd
import pandas.testing as pdt
import scanpy as sc
import scipy.sparse as sp
import scirpy as ir
from mudata import MuData

import scirpy as ir

DATASET_DIR = Path("/data/datasets/Maynard_Bivona_2020_NSCLC/")
# -

# The dataset has been downloaded from ENA and then processed using BraCer, TraCeR and the nf-core RNA-seq pipeline
# using `salmon`.
# using `salmon`.
#
# We previously processed this dataset with STAR + featureCounts, but following up the discussion
# on nf-core, featureCounts is not state-of-the art any more for estimating transcript abundances.
# We, therefore, switched to the nf-core RNA-seq pipeline and Salmon.
# on nf-core, featureCounts is not state-of-the art any more for estimating transcript abundances.
# We, therefore, switched to the nf-core RNA-seq pipeline and Salmon.

with open(
"/data/genomes/hg38/annotation/gencode/gencode.v33.primary_assembly.annotation.gtf",
Expand Down Expand Up @@ -201,7 +201,7 @@ def read_salmon(path):

adata = sc.AnnData(
var=res[0]["var"],
X=tpm_mat[sample_id_mask, :],
X=tpm_mat[sample_id_mask, :],
obs=meta.loc[sample_ids[sample_id_mask], :],
)

Expand Down Expand Up @@ -245,7 +245,9 @@ def read_salmon(path):
mdata_vis, groupby="airr:receptor_subtype", target_col="gex:patient"
)

_ = ir.pl.group_abundance(mdata_vis, groupby="airr:chain_pairing", target_col="gex:patient")
_ = ir.pl.group_abundance(
mdata_vis, groupby="airr:chain_pairing", target_col="gex:patient"
)

# ## save MuData

Expand Down

0 comments on commit 8ad2e8b

Please sign in to comment.