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()