Estimate simple simulation parameters

Estimate simulation parameters for the simple simulation from a real dataset.

simpleEstimate(counts, params = newSimpleParams())

# S3 method for SingleCellExperiment
simpleEstimate(counts, params = newSimpleParams())

# S3 method for matrix
simpleEstimate(counts, params = newSimpleParams())

Arguments

counts	either a counts matrix or a SingleCellExperiment object containing count data to estimate parameters from.
params	SimpleParams object to store estimated values in.

Value

SimpleParams object containing the estimated parameters.

Details

The nGenes and nCells parameters are taken from the size of the input data. The mean parameters are estimated by fitting a gamma distribution to the library size normalised mean expression level using fitdist. See SimpleParams for more details on the parameters.

Examples

# Load example data
library(scater)
set.seed(1)
sce <- mockSCE()

params <- simpleEstimate(sce)
params
#> A Params object of class SimpleParams 
#> Parameters can be (estimable) or [not estimable], 'Default' or  'NOT DEFAULT' 
#> Secondary parameters are usually set during simulation
#> 
#> Global: 
#> (GENES)  (CELLS)   [Seed] 
#>    2000      200   121054 
#> 
#> 3 additional parameters 
#> 
#> Mean: 
#>              (RATE)              (SHAPE) 
#> 0.00315146081285586    0.562634238245172 
#> 
#> Counts: 
#> [Dispersion] 
#>          0.1 
#>