|
229 | 229 | }
|
230 | 230 | ],
|
231 | 231 | "source": [
|
232 |
| - "sc.pp.filter_cells(adata, min_genes=200) # this does nothing, in this specific case\n", |
| 232 | + "sc.pp.filter_cells(adata, min_genes=200) # this does nothing, in this specific case\n", |
233 | 233 | "sc.pp.filter_genes(adata, min_cells=3)\n",
|
234 | 234 | "adata"
|
235 | 235 | ]
|
|
389 | 389 | "metadata": {},
|
390 | 390 | "outputs": [],
|
391 | 391 | "source": [
|
392 |
| - "adata = adata[(adata.obs.n_genes_by_counts < 2500) & (adata.obs.n_genes_by_counts > 200) & (adata.obs.pct_counts_mt < 5), :].copy()\n", |
| 392 | + "adata = adata[\n", |
| 393 | + " (adata.obs.n_genes_by_counts < 2500)\n", |
| 394 | + " & (adata.obs.n_genes_by_counts > 200)\n", |
| 395 | + " & (adata.obs.pct_counts_mt < 5),\n", |
| 396 | + " :,\n", |
| 397 | + "].copy()\n", |
393 | 398 | "adata.layers[\"counts\"] = adata.X.copy()"
|
394 | 399 | ]
|
395 | 400 | },
|
|
484 | 489 | }
|
485 | 490 | ],
|
486 | 491 | "source": [
|
487 |
| - "sc.pp.highly_variable_genes(adata, layer=\"counts\", n_top_genes=2000, min_mean=0.0125, max_mean=3, min_disp=0.5, flavor='seurat_v3')" |
| 492 | + "sc.pp.highly_variable_genes(\n", |
| 493 | + " adata,\n", |
| 494 | + " layer=\"counts\",\n", |
| 495 | + " n_top_genes=2000,\n", |
| 496 | + " min_mean=0.0125,\n", |
| 497 | + " max_mean=3,\n", |
| 498 | + " min_disp=0.5,\n", |
| 499 | + " flavor=\"seurat_v3\",\n", |
| 500 | + ")" |
488 | 501 | ]
|
489 | 502 | },
|
490 | 503 | {
|
|
544 | 557 | }
|
545 | 558 | ],
|
546 | 559 | "source": [
|
547 |
| - "adata.layers[\"normalized\"] = adata.X.copy()\n", |
| 560 | + "adata.layers[\"normalized\"] = adata.X.copy()\n", |
548 | 561 | "adata_subset_hvg = adata[:, adata.var[\"highly_variable\"]].copy()\n",
|
549 | 562 | "sc.pp.regress_out(adata_subset_hvg, [\"total_counts\", \"pct_counts_mt\"])\n",
|
550 | 563 | "sc.pp.scale(adata_subset_hvg, max_value=10)"
|
|
1507 | 1520 | }
|
1508 | 1521 | ],
|
1509 | 1522 | "source": [
|
1510 |
| - "sc.tl.rank_genes_groups(adata_subset_hvg, \"leiden\", groups=[\"0\"], reference=\"1\", method=\"wilcoxon\")\n", |
| 1523 | + "sc.tl.rank_genes_groups(\n", |
| 1524 | + " adata_subset_hvg, \"leiden\", groups=[\"0\"], reference=\"1\", method=\"wilcoxon\"\n", |
| 1525 | + ")\n", |
1511 | 1526 | "sc.pl.rank_genes_groups(adata_subset_hvg, groups=[\"0\"], n_genes=20)"
|
1512 | 1527 | ]
|
1513 | 1528 | },
|
|
0 commit comments