Esempio n. 1
0
def plot_matrix_muts(reactivity_type, max_val = 0.025):
    matrix_data = cotrans_matrix_data(reactivity_type, max_val = max_val)
    plt.matshow(matrix_data)
    ax = plt.gca()
    ax.grid(color='grey',linestyle='-',linewidth='0.5')
    ax.xaxis.set_ticks_position('bottom')
    plt.xlabel("Nucleotide (nt)")
    plt.ylabel("RNA Length (nt)")
    cbar = plt.colorbar(fraction=0.046, pad=0.04)
    cbar.set_label(reactivity_type)
    return plt
Esempio n. 2
0
def plot_cotrans_untreated():
    matrix_data = cotrans_matrix_data('untreated_counts', max_val = 5000)
    plt.matshow(matrix_data,cmap='jet')
    ax = plt.gca()
    ax.grid(color='grey',linestyle='-',linewidth='0.5')
    ax.xaxis.set_ticks_position('bottom')
    plt.xlabel("Nucleotide (nt)")
    plt.ylabel("RNA Length (nt)")
    cbar = plt.colorbar(fraction=0.046, pad=0.04)
    cbar.set_label('Counts')
    plt.title('Untreated Channel Counts')
    return plt