Guide Claude through ingesting TCGA sample sheets, expression archives, and clinical carts into omicverse, initialising survival metadata, and exporting annotated AnnData files.
This skill automates the preprocessing of TCGA (The Cancer Genome Atlas) data using the Omicverse package, providing a streamlined method for integrating and analyzing expression and clinical data. It specifically handles the ingestion of TCGA sample sheets, expression archives, and clinical carts, while also initializing survival metadata and exporting annotated AnnData files for downstream analysis. This solution is designed to facilitate the early stages of bioinformatics workflows, saving time by automating data processing tasks that would otherwise be tedious and error-prone when performed manually. It allows users to focus on higher-level analysis and interpretation rather than the intricacies of data preprocessing.
The skill offers a comprehensive set of features for TCGA data, including initialization of raw count, FPKM, and TPM matrices, automatic attachment of clinical metadata, and survival analysis. Users can quickly generate survival curves for individual genes or process all genes for survival analysis. It also allows exporting enriched AnnData files for further analysis and sharing with collaborators. This tool is particularly valuable for bioinformaticians, cancer researchers, and data scientists working with TCGA datasets. The integration with Omicverse ensures compatibility with other analysis tools and easy integration into broader omics workflows.
You will need the `gdc_sample_sheet.<date>.tsv` from the TCGA Sample Sheet export, the decompressed `gdc_download_xxxxx` directory containing expression archives, and the `clinical.cart.<date>` directory with clinical XML/JSON files.
You can call `aml_tcga.survial_analysis_all()` to perform survival analysis across all genes in the dataset. Please note that this may take some time depending on the dataset size.
Yes, you can reload a saved AnnData file by calling `aml_tcga.adata_read(<path>)` to reconstruct the class with the same paths and data.
If clinical dates are missing, survival plots cannot be generated. You should check the contents of the `clinical_cart` to ensure all necessary data is present.
Ensure that the case IDs match between the sample sheet and the expression files. If there is an issue with the IDs, you may need to re-download the data to ensure they align.
Quick Setup:
.claude/skills/Repository
Starlitnightly/omicverse