def benjamini(self):
genes = self.values()
p = {}
for id, gene in genes:
p[id] = gene.pvalue
benjamini_pvalues = calc_benjamini_hochberg_corrections(p.values(), len(p))
for index, in enumerate(benjamini_pvalues):
gene = genes[index]
gene.benjamini = benjamini_pvalues[0]
self[gene.id] = gene
def benjamini(self):
genes = self.values()
p = {}
for id, gene in genes:
p[id] = gene.pvalue
benjamini_pvalues = calc_benjamini_hochberg_corrections(
p.values(), len(p))
for index, in enumerate(benjamini_pvalues):
gene = genes[index]
gene.benjamini = benjamini_pvalues[0]
self[gene.id] = gene
def calc_benjamini(signature):
"""Function which preforms the actual Benjamini correction."""
transcripts = signature.transcripts.all()
p = [] # seq_id - p-values mapping.
t = [] #{} # seq_id - transcript mapping.
for transcript in transcripts.order_by('pvalue'):
t.append(transcript)
p.append(transcript.pvalue)
benjamini_pvalues = calc_benjamini_hochberg_corrections(p, len(p))
for index, benjamini_pvalue in enumerate(benjamini_pvalues):
#print index, benjamini_pvalue
t[index].benjamini = benjamini_pvalue[1] #expression__
t[index].save()
def calc_benjamini(signature):
"""Function which preforms the actual Benjamini correction."""
transcripts = signature.transcripts.all()
p = [] # seq_id - p-values mapping.
t = [] #{} # seq_id - transcript mapping.
for transcript in transcripts.order_by('pvalue'):
t.append(transcript)
p.append(transcript.pvalue)
benjamini_pvalues = calc_benjamini_hochberg_corrections(p, len(p))
for index, benjamini_pvalue in enumerate(benjamini_pvalues):
#print index, benjamini_pvalue
t[index].benjamini = benjamini_pvalue[1] #expression__
t[index].save()
