In [1]:

    

%matplotlib notebook


    

In [2]:

    

import manhattan_plot
import random


    

In [3]:

    

j = {}
for c in range(1, 5):
    j['chr' + str(c)] = []
    for i in range(100, int((1/c)*5e5)):
        j['chr' + str(c)].append((i, random.random()))
j['chrX'] = []
for i in range(100, 10000):
    j['chrX'].append((i, random.random()))
j['chrY'] = []
for i in range(100, 5000):
    j['chrY'].append((i, random.random()))
j['chrMT'] = []
for i in range(100, 3000):
    j['chrMT'].append((i, random.random()))


    

In [4]:

    

# Create a significant area
l = len(j['chr1'])
first = True
for i in range(l, int(l+l*0.05)):
    if first:
        rando = 1/random.randrange(1e4, 1e9)
        first = False
    else:
        rando = random.random()
        first = True
    j['chr1'].append((i, rando))
l = len(j['chr1'])
for i in range(l, int(l+l*0.3)):
    j['chr1'].append((i, random.random()))


    

In [10]:

    

plt = manhattan_plot.plot(j, title='Test Plot', sig_line=0.05)
plt.show()


    

In [9]:

    

plt = manhattan_plot.plot(j, title='As a line graph', sig_line=0.05, line_graph=True)
plt.show()


    

In [11]:

    

plt = manhattan_plot.plot(j, title='Raw p-values', sig_line=0.05, log_scale=False)
plt.show()


    

In [ ]: