LSHTM_analysis/scripts/temp_results.py

34 lines
1.6 KiB
Python
Executable file

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Thu Nov 12 15:20:32 2020
@author: tanu
"""
print('============================================'
, '\nQuick summary output for', drug, 'and' , gene.lower()
, '\n============================================'
, '\nTotal samples:', total_samples
, '\nNo. of Sus and Res', drug, 'samples:\n', meta_data[drug].value_counts()
, '\n'
, '\nPercentage of Sus and Res', drug, 'samples:\n', meta_data[drug].value_counts(normalize = True)*100
, '\n'
, '\nTotal no.of unique dr muts:', mut_grouped[dr_muts_col]
, '\nTotal no.of unique other muts:', mut_grouped[other_muts_col]
, '\nTotal no.of unique missense muts:', gene_LF1['mutationinformation'].nunique()
, '\nTotal no.of unique positions associated with missense muts:',gene_LF1['position'].nunique()
, '\nTotal no. of samples with missense muts:', len(gene_LF1)
, '\nTotal no. of unique samples with missense muts:', gene_LF1['id'].nunique()
, '\n'
, '\nTotal no.of samples with common_ids:', nu_common_ids['id']
, '\nTotal no.of samples with ambiguous muts:', len(inspect)
#, '\nTotal no.of unique ambiguous muts:', len(common_muts)
, '\nTotal no.of unique ambiguous muts:', inspect['mutation'].nunique()
, '\n============================================================='
, '\n\n\n')
print('\n----------------------------------'
, '\nWriting subsetted gwas data:', gene
, '\n----------------------------------'
, '\nFile', check_file
, '\nDim:', meta_data.shape)