INNER CODE UNIT · Python

total_truth_err

wheaton5/souporcell · consensus.py:395

                total_truth_err = logsumexp([truth,err])
                err = np.exp(err-total_truth_err)
                #print(err)
                logpost = [x - sumexp for x in genotypes]
                posteriors = np.exp(logpost)
                gn = []
                for g in logpost:
                    if math.isnan(g):
                        gn.append("NaN")
                    else:
                        gn.append(str(int(g)))
                gn = ",".join(gn)

                largest = np.argmax(posteriors)
                if err > 0.5:
                    newrec.FILTER = ['BACKGROUND']
                #print(posteriors)
                if len(posteriors) == 3:

View source record →