findParams {ROSeq} | R Documentation

## Finds the optimal values of parameters a and b that model the probability distribution of ranks, by Maximising the Log-Likelihood

### Description

Takes in as input the read count data corresponding to one sub-population and the typical gene statistics. Then it splits the entire range into equally sized bins of size k * σ, where k is a scalar with a default value of 0.05, and σ is the standard deviation of the pulled expression estimates across the cell-groups. Each of these bins corresponds to a rank. Therefore, for each group, cell frequency for each bin maps to a rank. These frequencies are normalized group-wise by dividing by the total cell count within a concerned group.

### Usage

    findParams(ds, geneStats)

### Arguments

`ds` | The (normalized and filtered) read count data corresponding to a sub-population
`geneStats` | A vector containing 7 values corresponding to the gene data (maximum, minimum, mean, standard deviation, upper multiple of standard deviation, lower multiple of standard deviation and log_2(fold change))

### Value

results A vector containing 5 values (a, b, A, number of bins, R2)

[Package ROSeq version 1.0.0 Index]