Select a subset of models from a 'celda_list' object generated by celdaGridSearch.

Convenience function for picking out specific models from a celda_list object. Models can be selected by various parameters, most importantly the K/L parameters (number of cell clusters / number of feature clusters).

subsetCeldaList(celda.list, params)

Arguments

celda.list	celda.list Object of class "celda_list". An object containing celda models returned from `celdaGridSearch()`.
params	List. List of parameters used to subset celda.list.

Value

A new 'celda_list' object containing all models matching the provided criteria in 'params'. If entry in the list matches, the results for the matching model will be returned.

Examples

celda.sim = simulateCells(model="celda_CG", K=5, L=10, G=100)
cgs = celdaGridSearch(celda.sim$counts, model="celda_CG", params.test=list(K=4:6, L=9:11),
                      params.fixed=list(sample.label=celda.sim$sample.label),
                      best.only=TRUE, nchains=1)
#> --------------------------------------------------------------------
#> Starting celdaGridSearch with celda_CG
#> Number of cores: 1
#> --------------------------------------------------------------------
#> Error in {    current.run = c(run.params[i, ])    chain.params = list()    for (j in names(params.test)) {        chain.params[[j]] = current.run[[j]]    }    chain.params$counts = counts    chain.params$seed = all.seeds[ifelse(i%%nchains == 0, nchains,         i%%nchains)]    chain.params$max.iter = max.iter    chain.params$nchain = 1    chain.params$count.checksum = count.checksum    chain.params$verbose = verbose    chain.params$logfile = paste0(logfile.prefix, "_", paste(paste(colnames(run.params),         run.params[i, ], sep = "-"), collapse = "_"), "_Seed-",         chain.params$seed, "_log.txt")    res = do.call(model, c(chain.params, params.fixed))    return(list(res))}: task 1 failed - "unused argument (gamma)"
res.K5.L10 = subsetCeldaList(cgs, params=list(K=5, L=10))
#> Error in isTRUE(class(celda.list)[1] == "celda_list"): object 'cgs' not found