INNER CODE UNIT · Python

sumexp

wheaton5/souporcell · consensus.py:379

                sumexp = logsumexp(genotypes)#[math.exp(x) for x in genotypes])
                #//print(genotypes)
                go = []
                for g in genotypes:
                    if math.isnan(g):
                        go.append('NaN')
                    else:
                        go.append(str(int(g)))
                    
                go = ",".join(go)
                
                #if sumexp == 0:
                #    posteriors = [0.333 for x in range(len(genotypes))]
                #else:
                truth = fit['truth'][locus_index[locus]]
                err = fit['err'][locus_index[locus]]
                total_truth_err = logsumexp([truth,err])
                err = np.exp(err-total_truth_err)

View source record →