Estimate simulation parameters for the mfa simulation from a real dataset.
Usage
mfaEstimate(counts, params = newMFAParams())
# S3 method for class 'SingleCellExperiment'
mfaEstimate(counts, params = newMFAParams())
# S3 method for class 'matrix'
mfaEstimate(counts, params = newMFAParams())
Details
The nGenes
and nCells
parameters are taken from the size of the
input data. The dropout lambda parameter is estimate using
empirical_lambda
. See MFAParams
for more
details on the parameters.
Examples
# Load example data
if (requireNamespace("mfa", quietly = TRUE)) {
library(mfa)
synth <- create_synthetic(
C = 20, G = 5, zero_negative = TRUE,
model_dropout = TRUE
)
params <- mfaEstimate(synth$X)
params
}
#> A Params object of class MFAParams
#> Parameters can be (estimable) or [not estimable], 'Default' or 'NOT DEFAULT'
#> Secondary parameters are usually set during simulation
#>
#> Global:
#> (GENES) (CELLS) [SEED]
#> 20 5 875438
#>
#> 4 additional parameters
#>
#> Transient:
#> [Proportion]
#> 0
#>
#> Negative:
#> [Zero]
#> TRUE
#>
#> Dropout:
#> [Present] (LAMBDA)
#> FALSE 0.270058531928805
#>