Skip to content

mecfs_bio.analysis.red_blood_cell_volume

Analyze red blood cell volume GWAS

Functions:

run_red_blood_cell_volume_analysis

run_red_blood_cell_volume_analysis()

Analyze GWAS of red blood cell volume from a million veterans Includes: - Apply MAGMA to estimate significant genes - Combine Magma results with GTEx to estimate significant tissues - Apply S-LDSC using standard reference datasets to estimate key cells and tissues

Source code in mecfs_bio/analysis/red_blood_cell_volume.py
def run_red_blood_cell_volume_analysis():
    """
    Analyze GWAS of red blood cell volume from a million veterans
    Includes:
    - Apply MAGMA to estimate significant genes
    - Combine Magma results with GTEx to estimate significant tissues
    - Apply S-LDSC using standard reference datasets to estimate key cells and tissues


    """
    DEFAULT_RUNNER.run(
        [
            MILLION_VETERANS_EUR_RBC_VOLUME_MAGMA_TASKS.inner.bar_plot_task,
            MILLION_VETERANS_EUR_RBC_VOLUME_MAGMA_TASKS.inner.gene_analysis_task,
            MILLION_VETERANS_EUR_RBC_VOLUME_MAGMA_TASKS.inner.labeled_filtered_gene_analysis_task,
        ]
        + MILLION_VETERANS_EUR_RBC_VOL_S_LDSC_TASKS.get_terminal_tasks(),
        incremental_save=True,
    )