Kersplat_simulation | R Documentation |
This function is used to simulate datasets from learned parameters by kersplatSimulate
function in Splatter package.
Kersplat_simulation(
parameters,
other_prior = NULL,
return_format,
verbose = FALSE,
seed
)
parameters |
A object generated by |
other_prior |
A list with names of certain parameters. Some methods need
extra parameters to execute the estimation step, so you must input them. In
simulation step, the number of cells, genes, groups, batches, the percent of
DEGs and other variables are usually customed, so before simulating a dataset
you must point it out. See |
return_format |
A character. Alternative choices: list, SingleCellExperiment,
Seurat, h5ad. If you select |
verbose |
Logical. Whether to return messages or not. |
seed |
A random seed. |
In addtion to simulate datasets with default parameters, users can set two extra parameters before executing simulation step within Simple method:
nCells. Just only type other_prior = list(nCells = n)
when you want to
simulate a dataset with n cells.
nGenes. Like nCells
, users just only type other_prior = list(nGenes = m)
to get the dataset with m genes.
Zappia L, Phipson B, Oshlack A. Splatter: simulation of single-cell RNA sequencing data. Genome biology, 2017, 18(1): 1-15. https://doi.org/10.1186/s13059-017-1305-0
Bioconductor URL: https://bioconductor.org/packages/release/bioc/html/splatter.html
Github URL: https://github.com/Oshlack/splatter
## Not run:
ref_data <- simmethods::data
# Estimate parameters
estimate_result <- Kersplat_estimation(ref_data = ref_data,
verbose = TRUE,
seed = 111)
# (1) Simulate a dataset with default parameters
simulate_result <- Kersplat_simulation(parameters = estimate_result[["estimate_result"]],
return_format = "list",
verbose = TRUE,
seed = 111)
counts <- simulate_result[["simulate_result"]][["count_data"]]
dim(counts)
# (2) Simulate a dataset with customed number of cells and genes
simulate_result <- Kersplat_simulation(parameters = estimate_result[["estimate_result"]],
return_format = "list",
other_prior = list(nCells = 500,
nGenes = 3000),
verbose = TRUE,
seed = 111)
counts <- simulate_result[["simulate_result"]][["count_data"]]
dim(counts)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.