def raw_value_comparison(coh, plot=False):
    "Return the average differences in raw copy number values between the\n    gene-level calls in hg19 and hg38 for each gene for a given tumor type \n    'coh.' If plot=True, plot the genes' differences in a histogram."
    (df_38, df_19) = load_data(coh, thresh=False)
    df_s = (df_38 - df_19)
    avg_diff = {g: np.average(df_s[g]) for g in df_s.columns.get_level_values('Gene Symbol')}
    results = []
    std = np.std([avg_diff[x] for x in avg_diff])
    for g in avg_diff:
        if (avg_diff[g] > (4 * std)):
            results.append([coh, 'Pos', g, avg_diff[g]])
        elif (avg_diff[g] < ((- 4) * std)):
            results.append([coh, 'Neg', g, avg_diff[g]])
    if plot:
        plt.hist([avg_diff[x] for x in avg_diff], bins=1000)
        plt.title(coh, fontsize=16)
        plt.xlabel('Average CN Difference Between Alignments', fontsize=14)
        plt.ylabel('Genes', fontsize=14)
        sns.despine()
        plt.savefig((('./genehists/' + coh) + '_genehist.pdf'))
        plt.savefig((('./genehists/' + coh) + '_genehist.png'))
        plt.clf()
    return results
Return the average differences in raw copy number values between the
gene-level calls in hg19 and hg38 for each gene for a given tumor type 
'coh.' If plot=True, plot the genes' differences in a histogram.
def raw_value_comparison(coh, plot=False):
    "Return the average differences in raw copy number values between the\n    gene-level calls in hg19 and hg38 for each gene for a given tumor type \n    'coh.' If plot=True, plot the genes' differences in a histogram."
    (df_38, df_19) = load_data(coh, thresh=False)
    df_s = (df_38 - df_19)
    avg_diff = {g: np.average(df_s[g]) for g in df_s.columns.get_level_values('Gene Symbol')}
    results = []
    std = np.std([avg_diff[x] for x in avg_diff])
    for g in avg_diff:
        if (avg_diff[g] > (4 * std)):
            results.append([coh, 'Pos', g, avg_diff[g]])
        elif (avg_diff[g] < ((- 4) * std)):
            results.append([coh, 'Neg', g, avg_diff[g]])
    if plot:
        plt.hist([avg_diff[x] for x in avg_diff], bins=1000)
        plt.title(coh, fontsize=16)
        plt.xlabel('Average CN Difference Between Alignments', fontsize=14)
        plt.ylabel('Genes', fontsize=14)
        sns.despine()
        plt.savefig((('./genehists/' + coh) + '_genehist.pdf'))
        plt.savefig((('./genehists/' + coh) + '_genehist.png'))
        plt.clf()
    return results