|
1173 | 1173 | "name": "stderr", |
1174 | 1174 | "output_type": "stream", |
1175 | 1175 | "text": [ |
1176 | | - "\r", |
| 1176 | + "\r\n", |
1177 | 1177 | "0.00B [00:00, ?B/s]" |
1178 | 1178 | ] |
1179 | 1179 | }, |
1180 | 1180 | { |
1181 | 1181 | "name": "stderr", |
1182 | 1182 | "output_type": "stream", |
1183 | 1183 | "text": [ |
1184 | | - "\r", |
| 1184 | + "\r\n", |
1185 | 1185 | "56.5kB [00:00, 168MB/s]" |
1186 | 1186 | ] |
1187 | 1187 | }, |
|
1204 | 1204 | "name": "stderr", |
1205 | 1205 | "output_type": "stream", |
1206 | 1206 | "text": [ |
1207 | | - "\r", |
| 1207 | + "\r\n", |
1208 | 1208 | "0.00B [00:00, ?B/s]" |
1209 | 1209 | ] |
1210 | 1210 | }, |
1211 | 1211 | { |
1212 | 1212 | "name": "stderr", |
1213 | 1213 | "output_type": "stream", |
1214 | 1214 | "text": [ |
1215 | | - "\r", |
| 1215 | + "\r\n", |
1216 | 1216 | "449kB [00:00, 4.59MB/s]" |
1217 | 1217 | ] |
1218 | 1218 | }, |
1219 | 1219 | { |
1220 | 1220 | "name": "stderr", |
1221 | 1221 | "output_type": "stream", |
1222 | 1222 | "text": [ |
1223 | | - "\r", |
| 1223 | + "\r\n", |
1224 | 1224 | "1.78MB [00:00, 6.86MB/s]" |
1225 | 1225 | ] |
1226 | 1226 | }, |
1227 | 1227 | { |
1228 | 1228 | "name": "stderr", |
1229 | 1229 | "output_type": "stream", |
1230 | 1230 | "text": [ |
1231 | | - "\r", |
| 1231 | + "\r\n", |
1232 | 1232 | "4.76MB [00:00, 15.6MB/s]" |
1233 | 1233 | ] |
1234 | 1234 | }, |
1235 | 1235 | { |
1236 | 1236 | "name": "stderr", |
1237 | 1237 | "output_type": "stream", |
1238 | 1238 | "text": [ |
1239 | | - "\r", |
| 1239 | + "\r\n", |
1240 | 1240 | "5.91MB [00:00, 15.9MB/s]" |
1241 | 1241 | ] |
1242 | 1242 | }, |
|
1398 | 1398 | ], |
1399 | 1399 | "source": [ |
1400 | 1400 | "# Query Omnipath and get PanglaoDB\n", |
1401 | | - "markers = dc.get_resource(name=\"PanglaoDB\", organism=\"human\")\n", |
| 1401 | + "markers = dc.op.resource(name=\"PanglaoDB\", organism=\"human\")\n", |
1402 | 1402 | "# Keep canonical cell type markers alone\n", |
1403 | 1403 | "markers = markers[markers[\"canonical_marker\"]]\n", |
1404 | 1404 | "\n", |
|
1424 | 1424 | "name": "stderr", |
1425 | 1425 | "output_type": "stream", |
1426 | 1426 | "text": [ |
1427 | | - "\r", |
| 1427 | + "\r\n", |
1428 | 1428 | " 0%| | 0/2 [00:00<?, ?it/s]" |
1429 | 1429 | ] |
1430 | 1430 | }, |
1431 | 1431 | { |
1432 | 1432 | "name": "stderr", |
1433 | 1433 | "output_type": "stream", |
1434 | 1434 | "text": [ |
1435 | | - "\r", |
| 1435 | + "\r\n", |
1436 | 1436 | " 50%|█████ | 1/2 [00:05<00:05, 5.44s/it]" |
1437 | 1437 | ] |
1438 | 1438 | }, |
1439 | 1439 | { |
1440 | 1440 | "name": "stderr", |
1441 | 1441 | "output_type": "stream", |
1442 | 1442 | "text": [ |
1443 | | - "\r", |
| 1443 | + "\r\n", |
1444 | 1444 | "100%|██████████| 2/2 [00:09<00:00, 4.45s/it]" |
1445 | 1445 | ] |
1446 | 1446 | }, |
1447 | 1447 | { |
1448 | 1448 | "name": "stderr", |
1449 | 1449 | "output_type": "stream", |
1450 | 1450 | "text": [ |
1451 | | - "\r", |
| 1451 | + "\r\n", |
1452 | 1452 | "100%|██████████| 2/2 [00:09<00:00, 4.60s/it]" |
1453 | 1453 | ] |
1454 | 1454 | }, |
|
1461 | 1461 | } |
1462 | 1462 | ], |
1463 | 1463 | "source": [ |
1464 | | - "dc.run_mlm(mat=adata, net=markers, weight=None, source=\"cell_type\", target=\"genesymbol\", verbose=True, use_raw=False)" |
| 1464 | + "dc.mt.mlm(adata, net=markers.rename(columns=dict(cell_type=\"source\", genesymbol=\"target\")), verbose=True)" |
1465 | 1465 | ] |
1466 | 1466 | }, |
1467 | 1467 | { |
1468 | 1468 | "cell_type": "markdown", |
1469 | 1469 | "id": "2e01e50a-238f-459f-a9f6-f25880b78a2c", |
1470 | 1470 | "metadata": {}, |
1471 | 1471 | "source": [ |
1472 | | - "The obtained results are stored in the .obsm key, with `mlm_estimate` representing coefficient t-values:" |
| 1472 | + "The obtained results are stored in `.obsm`, with `\"score_mlm\"` representing coefficient t-values, and `\"padj_mlm\"` representing adjusted p-values.:" |
1473 | 1473 | ] |
1474 | 1474 | }, |
1475 | 1475 | { |
|
1714 | 1714 | } |
1715 | 1715 | ], |
1716 | 1716 | "source": [ |
1717 | | - "adata.obsm[\"mlm_estimate\"].head()" |
| 1717 | + "adata.obsm[\"score_mlm\"].head()" |
1718 | 1718 | ] |
1719 | 1719 | }, |
1720 | 1720 | { |
|
1749 | 1749 | } |
1750 | 1750 | ], |
1751 | 1751 | "source": [ |
1752 | | - "acts = dc.get_acts(adata=adata, obsm_key=\"mlm_estimate\")\n", |
| 1752 | + "acts = dc.pp.get_obsm(adata, \"score_mlm\")\n", |
1753 | 1753 | "sc.pl.umap(\n", |
1754 | 1754 | " acts,\n", |
1755 | 1755 | " color=[\n", |
|
1780 | 1780 | "metadata": {}, |
1781 | 1781 | "outputs": [], |
1782 | 1782 | "source": [ |
1783 | | - "mean_enr = dc.summarize_acts(acts, groupby=\"leiden_res0_5\", min_std=1.0)\n", |
1784 | | - "annotation_dict = dc.assign_groups(mean_enr)\n", |
| 1783 | + "enr = dc.tl.rankby_group(acts, groupby=\"leiden_res0_5\")\n", |
| 1784 | + "annotation_dict = enr[enr[\"stat\"] > 0].groupby(\"group\").head(1).set_index(\"group\")[\"name\"].to_dict()\n", |
1785 | 1785 | "adata.obs[\"dc_anno\"] = [annotation_dict[clust] for clust in adata.obs[\"leiden_res0_5\"]]" |
1786 | 1786 | ] |
1787 | 1787 | }, |
|
0 commit comments