def go_term_info(goid=None):
if goid is None:
flash("That page needs a GOID")
return redirect(url_for("homepage"))
go_info = es.go_info(goid)
genes = es.go_to_gene(goid)
if go_info is None or genes is None:
if go_info is None:
flash("You need to implement <code>go_info()</code>")
if genes is None:
flash("You need to implement <code>go_to_gene()</code>")
return redirect(url_for('homepage'))
genes = [(gene, es.gene_name(gene)) for gene in es.go_to_gene(goid)]
local_network = hs.go_network(goid, True, 1)
if local_network is None:
flash("You could maybe implement <code>go_network()</code>")
elif len(local_network['nodes']) > 1000:
flash("Too many nodes to plot a network, forget it")
local_network = {'nodes': [], 'links': []}
return render_template("go_term.html", go_info=go_info, genes=genes,
local_network=local_network)
def experiment(exp):
exp_data = es.experiment()
if exp_data is not None:
genes = [(es.gene_name(g) or 'Implement gene_name()!', v)
for g,v in exp_data[exp]]
fig_dict = hs.plot_experiment(genes)
else:
fig_dict = None
return render_template("experiment.html", exp=exp, fig_dict=fig_dict)
def go_term_info(goid):
go_info = es.go_info(goid)
genes = es.go_to_gene(goid)
if go_info is None or genes is None:
if go_info is None:
flash("You need to implement <code>go_info()</code>")
if genes is None:
flash("You need to implement <code>go_to_gene()</code>")
return redirect(url_for('homepage'))
genes = [(gene, es.gene_name(gene)) for gene in es.go_to_gene(goid)]
return render_template("go_term.html", go_info=go_info, genes=genes)
def gene_info(gene):
gene_name = es.gene_name(gene)
gene_info = es.gene_info(gene)
if gene_name is None or gene_info is None:
if gene_name is None:
flash("You need to implement <code>gene_name()</code>")
if gene_info is None:
flash("You need to implement <code>gene_info()</code>")
return redirect(url_for('homepage'))
go_terms = es.gene_to_go(gene)
go_terms = [(goid, es.go_info(goid)) for goid in go_terms]
return render_template("gene.html", gene_name=gene_name,
gene_info=gene_info, go_terms=go_terms)
def gene_list(exp, top_or_bottom, n=10):
exp_data = es.experiment()
if exp_data is not None:
exp_data = exp_data[exp]
else:
flash("You need to implement <code>experiment()</code> first!")
return redirect(url_for('experiment', exp=exp))
if top_or_bottom == 'top':
genes = sorted(exp_data, key=lambda g: g[1], reverse=True)[:n]
elif top_or_bottom == 'bottom':
genes = sorted(exp_data, key=lambda g: g[1])[:n]
else:
flash("That's not a valid URL")
return redirect(url_for('experiment', exp=exp))
return render_template("gene_list.html", label=top_or_bottom, n=n, exp=exp,
genes=[(g, es.gene_name(g), v) for g,v in genes])
def experiment(exp, n=10):
exp_data = es.experiment()
if exp_data is not None:
genes = [(es.gene_name(g) or 'Implement gene_name()!', v)
for g,v in exp_data[exp]]
fig_dict = hs.plot_experiment(genes)
if fig_dict is None:
flash("You could maybe implement <code>plot_experiment()</code>")
similar_exp = hs.similar_experiments(exp, exp_data, n)
if similar_exp is None:
flash("You could maybe implement <code>similar_experiments()</code>")
similar_exp = []
else:
flash("You should implement <code>experiment()</code> and <code>plot_experiment()</code>")
fig_dict = None
similar_exp = []
return render_template("experiment.html", exp=exp, n=n,
fig_dict=fig_dict, similar_exp=similar_exp)
def gene_info(gene=None):
if gene is None:
flash("That page needs a gene ID")
return redirect(url_for("homepage"))
gene_name = es.gene_name(gene)
gene_info = es.gene_info(gene)
if gene_name is None or gene_info is None:
if gene_name is None:
flash("You need to implement <code>gene_name()</code>")
if gene_info is None:
flash("You need to implement <code>gene_info()</code>")
return redirect(url_for('homepage'))
go_terms = es.gene_to_go(gene)
if go_terms is None:
flash("You need to implement <code>gene_to_go()</code>")
go_terms = []
else:
go_terms = [(goid, es.go_info(goid)) for goid in go_terms]
if go_terms and go_terms[0][1] is None:
flash("You need to implement <code>go_info()</code>")
gene_data = es.gene_data(gene)
if gene_data is None:
flash("You need to implement <code>gene_data()</code>")
hist_data = None
else:
gene_hist,bins = np.histogram(gene_data, bins=np.linspace(-8, 8, 33))
hist_data = zip(bins[:-1], gene_hist)
local_network = hs.go_network(gene, False, 2)
if local_network is None:
flash("You should implement <code>go_network()</code>")
elif len(local_network['nodes']) > 1000:
flash("Too many nodes to plot a network, forget it")
local_network = {'nodes': [], 'links': []}
return render_template("gene.html", gene_name=gene_name,
gene_info=gene_info, go_terms=go_terms,
hist_data=hist_data, local_network=local_network)