<jupyter_start><jupyter_text>ATAC<jupyter_code>#load pfm
atac_sum_6mer=np.loadtxt("atac.gm12878.0.computed.6mer.sum")
atac_sum_24mer=np.loadtxt("atac.gm12878.0.computed.24mer.sum")
atac_first_6mer=np.loadtxt("atac.gm12878.0.computed.6mer.first")
atac_first_24mer=np.loadtxt("atac.gm12878.0.computed.24mer.first")
atac_center_6mer=np.loadtxt("atac.gm12878.0.computed.6mer.center")
atac_center_24mer=np.loadtxt("atac.gm12878.0.computed.24mer.center")
#get pssm from pfm
atac_sum_6mer_pssm=np.log2(atac_sum_6mer/background)
atac_sum_24mer_pssm=np.log2(atac_sum_24mer/background)
atac_first_6mer_pssm=np.log2(atac_first_6mer/background)
atac_first_24mer_pssm=np.log2(atac_first_24mer/background)
atac_center_6mer_pssm=np.log2(atac_center_6mer/background)
atac_center_24mer_pssm=np.log2(atac_center_24mer/background)
atac_sum_6mer
atac_sum_6mer_pssm
f,axes=plt.subplots(4,3,dpi=80,figsize=(12,10))
matplotlib.pyplot.subplots_adjust(wspace=0.4,hspace=0.4)

axes[0,0]=plot_bases_on_ax(atac_sum_6mer.T,axes[0,0],show_ticks=True)
axes[0,0].set_title("6-mer, sum, PFM")
axes[1,0]=plot_bases_on_ax(atac_sum_6mer_pssm.T,axes[1,0],show_ticks=True)
axes[1,0].set_title("6-mer, sum, PSSM")

axes[2,0]=plot_bases_on_ax(atac_sum_24mer.T,axes[2,0],show_ticks=True)
axes[2,0].set_title("24-mer, sum, PFM")
axes[3,0]=plot_bases_on_ax(atac_sum_24mer_pssm.T,axes[3,0],show_ticks=True)
axes[3,0].set_title("24-mer, sum, PSSM")

axes[0,1]=plot_bases_on_ax(atac_first_6mer.T,axes[0,1],show_ticks=True)
axes[0,1].set_title("6-mer, first, PFM")
axes[1,1]=plot_bases_on_ax(atac_first_6mer_pssm.T,axes[1,1],show_ticks=True)
axes[1,1].set_title("6-mer, first, PSSM")

axes[2,1]=plot_bases_on_ax(atac_first_24mer.T,axes[2,1],show_ticks=True)
axes[2,1].set_title("24-mer, first, PFM")
axes[3,1]=plot_bases_on_ax(atac_first_24mer_pssm.T,axes[3,1],show_ticks=True)
axes[3,1].set_title("24-mer, first, PSSM")



axes[0,2]=plot_bases_on_ax(atac_center_6mer.T,axes[0,2],show_ticks=True)
axes[0,2].set_title("6-mer, center, PFM")
axes[1,2]=plot_bases_on_ax(atac_center_6mer_pssm.T,axes[1,2],show_ticks=True)
axes[1,2].set_title("6-mer, center, PSSM")

axes[2,2]=plot_bases_on_ax(atac_center_24mer.T,axes[2,2],show_ticks=True)
axes[2,2].set_title("24-mer, center, PFM")
axes[3,2]=plot_bases_on_ax(atac_center_24mer_pssm.T,axes[3,2],show_ticks=True)
axes[3,2].set_title("24-mer, center, PSSM")

plt.suptitle("ATAC, GM12878, fold 0")<jupyter_output><empty_output><jupyter_text>DNASE GM12878<jupyter_code>#load pfm
dnase_sum_6mer=np.loadtxt("dnase.gm12878.0.computed.6mer.sum")
dnase_sum_24mer=np.loadtxt("dnase.gm12878.0.computed.24mer.sum")
dnase_first_6mer=np.loadtxt("dnase.gm12878.0.computed.6mer.first")
dnase_first_24mer=np.loadtxt("dnase.gm12878.0.computed.24mer.first")
dnase_center_6mer=np.loadtxt("dnase.gm12878.0.computed.6mer.center")
dnase_center_24mer=np.loadtxt("dnase.gm12878.0.computed.24mer.center")
#get pssm from pfm
dnase_sum_6mer_pssm=np.log2(dnase_sum_6mer/background)
dnase_sum_24mer_pssm=np.log2(dnase_sum_24mer/background)
dnase_first_6mer_pssm=np.log2(dnase_first_6mer/background)
dnase_first_24mer_pssm=np.log2(dnase_first_24mer/background)
dnase_center_6mer_pssm=np.log2(dnase_center_6mer/background)
dnase_center_24mer_pssm=np.log2(dnase_center_24mer/background)
f,axes=plt.subplots(4,3,dpi=80,figsize=(12,10))
matplotlib.pyplot.subplots_adjust(wspace=0.4,hspace=0.4)

axes[0,0]=plot_bases_on_ax(dnase_sum_6mer.T,axes[0,0],show_ticks=True)
axes[0,0].set_title("6-mer, sum, PFM")
axes[1,0]=plot_bases_on_ax(dnase_sum_6mer_pssm.T,axes[1,0],show_ticks=True)
axes[1,0].set_title("6-mer, sum, PSSM")

axes[2,0]=plot_bases_on_ax(dnase_sum_24mer.T,axes[2,0],show_ticks=True)
axes[2,0].set_title("24-mer, sum, PFM")
axes[3,0]=plot_bases_on_ax(dnase_sum_24mer_pssm.T,axes[3,0],show_ticks=True)
axes[3,0].set_title("24-mer, sum, PSSM")

axes[0,1]=plot_bases_on_ax(dnase_first_6mer.T,axes[0,1],show_ticks=True)
axes[0,1].set_title("6-mer, first, PFM")
axes[1,1]=plot_bases_on_ax(dnase_first_6mer_pssm.T,axes[1,1],show_ticks=True)
axes[1,1].set_title("6-mer, first, PSSM")

axes[2,1]=plot_bases_on_ax(dnase_first_24mer.T,axes[2,1],show_ticks=True)
axes[2,1].set_title("24-mer, first, PFM")
axes[3,1]=plot_bases_on_ax(dnase_first_24mer_pssm.T,axes[3,1],show_ticks=True)
axes[3,1].set_title("24-mer, first, PSSM")



axes[0,2]=plot_bases_on_ax(dnase_center_6mer.T,axes[0,2],show_ticks=True)
axes[0,2].set_title("6-mer, center, PFM")
axes[1,2]=plot_bases_on_ax(dnase_center_6mer_pssm.T,axes[1,2],show_ticks=True)
axes[1,2].set_title("6-mer, center, PSSM")

axes[2,2]=plot_bases_on_ax(dnase_center_24mer.T,axes[2,2],show_ticks=True)
axes[2,2].set_title("24-mer, center, PFM")
axes[3,2]=plot_bases_on_ax(dnase_center_24mer_pssm.T,axes[3,2],show_ticks=True)
axes[3,2].set_title("24-mer, center, PSSM")

plt.suptitle("dnase, GM12878, fold 0")<jupyter_output><empty_output><jupyter_text>DNASE K562<jupyter_code>#load pfm
dnase_sum_6mer=np.loadtxt("dnase.k562.0.computed.6mer.sum")
dnase_sum_24mer=np.loadtxt("dnase.k562.0.computed.24mer.sum")
dnase_first_6mer=np.loadtxt("dnase.k562.0.computed.6mer.first")
dnase_first_24mer=np.loadtxt("dnase.k562.0.computed.24mer.first")
dnase_center_6mer=np.loadtxt("dnase.k562.0.computed.6mer.center")
dnase_center_24mer=np.loadtxt("dnase.k562.0.computed.24mer.center")
#get pssm from pfm
dnase_sum_6mer_pssm=np.log2(dnase_sum_6mer/background)
dnase_sum_24mer_pssm=np.log2(dnase_sum_24mer/background)
dnase_first_6mer_pssm=np.log2(dnase_first_6mer/background)
dnase_first_24mer_pssm=np.log2(dnase_first_24mer/background)
dnase_center_6mer_pssm=np.log2(dnase_center_6mer/background)
dnase_center_24mer_pssm=np.log2(dnase_center_24mer/background)
f,axes=plt.subplots(4,3,dpi=80,figsize=(12,10))
matplotlib.pyplot.subplots_adjust(wspace=0.4,hspace=0.4)

axes[0,0]=plot_bases_on_ax(dnase_sum_6mer.T,axes[0,0],show_ticks=True)
axes[0,0].set_title("6-mer, sum, PFM")
axes[1,0]=plot_bases_on_ax(dnase_sum_6mer_pssm.T,axes[1,0],show_ticks=True)
axes[1,0].set_title("6-mer, sum, PSSM")

axes[2,0]=plot_bases_on_ax(dnase_sum_24mer.T,axes[2,0],show_ticks=True)
axes[2,0].set_title("24-mer, sum, PFM")
axes[3,0]=plot_bases_on_ax(dnase_sum_24mer_pssm.T,axes[3,0],show_ticks=True)
axes[3,0].set_title("24-mer, sum, PSSM")

axes[0,1]=plot_bases_on_ax(dnase_first_6mer.T,axes[0,1],show_ticks=True)
axes[0,1].set_title("6-mer, first, PFM")
axes[1,1]=plot_bases_on_ax(dnase_first_6mer_pssm.T,axes[1,1],show_ticks=True)
axes[1,1].set_title("6-mer, first, PSSM")

axes[2,1]=plot_bases_on_ax(dnase_first_24mer.T,axes[2,1],show_ticks=True)
axes[2,1].set_title("24-mer, first, PFM")
axes[3,1]=plot_bases_on_ax(dnase_first_24mer_pssm.T,axes[3,1],show_ticks=True)
axes[3,1].set_title("24-mer, first, PSSM")



axes[0,2]=plot_bases_on_ax(dnase_center_6mer.T,axes[0,2],show_ticks=True)
axes[0,2].set_title("6-mer, center, PFM")
axes[1,2]=plot_bases_on_ax(dnase_center_6mer_pssm.T,axes[1,2],show_ticks=True)
axes[1,2].set_title("6-mer, center, PSSM")

axes[2,2]=plot_bases_on_ax(dnase_center_24mer.T,axes[2,2],show_ticks=True)
axes[2,2].set_title("24-mer, center, PFM")
axes[3,2]=plot_bases_on_ax(dnase_center_24mer_pssm.T,axes[3,2],show_ticks=True)
axes[3,2].set_title("24-mer, center, PSSM")

plt.suptitle("dnase, k562, fold 0")<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.17560099474993093, 'fraction_numerical': 0.08400110527770102, 'mean_word_length': 13.05631067961165, 'pattern_counts': {'":': 0, '<': 13, '<?xml version=': 0, '>': 13, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '3', 'dir': 'jupyter-structured-clean-dedup', 'id': '517273', 'license': 'MIT', 'n_tokens_mistral': 3908, 'n_tokens_neox': 3619, 'n_words': 218, 'path': 'enzymatic_bias/get_pwm_from_bpnet/.ipynb_checkpoints/BPNET_Model_PWM_from_Predictions-checkpoint.ipynb', 'repo_name': 'kundajelab/bias_correction', 'size': '7238'}