## ----------------------------------------------------------------------------- library(destiny) data(guo_norm) ## ----------------------------------------------------------------------------- sigmas <- find_sigmas(guo_norm, verbose = FALSE) optimal_sigma(sigmas) ## ----------------------------------------------------------------------------- palette(cube_helix(6)) plots <- lapply( list('local', 5, round(optimal_sigma(sigmas), 2), 100), function(sigma) plot( DiffusionMap(guo_norm, sigma), 1:2, main = paste('σ =', sigma), col_by = 'num_cells', draw_legend = FALSE)) do.call(gridExtra::grid.arrange, c(plots, ncol = 2))