fig = plt.figure(figsize=(10, 2.5))
NROW = 1
NCOL = 4
axs=[]
for item in range(0, NROW):
for col in range(0, NCOL):
axs.append(plt.subplot2grid((NROW, NCOL), (item, col)))
count=0
# flags for axis control
yaxis={'PP64':1, 'CPD':1}
xlabel={1000:0, 500:0, 100:1}
flank=200
# for proximal TFBS-DHS
pp='CPD'
target='tfbs-bs-seperated'
title='Low (n = 4,878)'
# set axisParam [title, xaxis, xlabel, Lyaxis, Lylabel, Ryaxis, Rylabel ]
axisParam = [title, 1, 0, 1, 1, 0, 0 ]
# plot values
cls.plot_TFBS_xrseq_combined(axs[count], pp, [ 'NHF1', 'CSB' ], "combined", flank, "allTFs_0", axisParam, target)
count+=1
title='Low-medium (n = 4,547)'
# set axisParam [title, xaxis, xlabel, Lyaxis, Lylabel, Ryaxis, Rylabel ]
axisParam = [title, 1, 0, 0, 0, 0, 0 ]
# plot values
cls.plot_TFBS_xrseq_combined(axs[count], pp, [ 'NHF1', 'CSB' ], "combined", flank, "allTFs_1", axisParam, target)
count+=1
title='Medium-high (n = 4,580)'
# set axisParam [title, xaxis, xlabel, Lyaxis, Lylabel, Ryaxis, Rylabel ]
axisParam = [title, 1, 0, 0, 0, 0, 0 ]
# plot values
cls.plot_TFBS_xrseq_combined(axs[count], pp, [ 'NHF1', 'CSB' ], "combined", flank, "allTFs_2", axisParam, target)
count+=1
title='High (n = 4,560)'
# set axisParam [title, xaxis, xlabel, Lyaxis, Lylabel, Ryaxis, Rylabel ]
axisParam = [title, 1, 0, 0, 0, 1, 1 ]
# plot values
cls.plot_TFBS_xrseq_combined(axs[count], pp, [ 'NHF1', 'CSB' ], "combined", flank, "allTFs_3", axisParam, target)
count+=1
plt.tight_layout(pad=0, w_pad=0)
plt.show()