INNER CODE UNIT · Python
myopen
wheaton5/souporcell · consensus.py:22
def myopen(fname): return open(fname, 'rb') if fname.endswith('.gz') else open(fname)
cell_genotype_consensus = """
data {
int<lower=0> cells; // number of cells
int<lower=0> loci; // number of loci
int<lower=0> msoup; //number of loci used to estimate soup
int<lower=0> k; // number of clusters
int<lower=1,upper=2> ploidy;
int cluster_allele_counts[loci, k, 2];
int cluster_allele_counts_soup[msoup,k,2];
int cluster_num_cells[k];
real average_allele_expression[loci, 2];
real average_allele_expression_soup[msoup,2];
}
transformed data {
real<upper=0> neg_log_3;
real<upper=0> fp_prior;