# #——包括= FALSE --------------------------------------------------------- knitr: opts_chunk美元集(崩溃= TRUE,发表评论 = "#>" ) ## ----------------------------------------------------------------------------- suppressPackageStartupMessages({库(waddR)图书馆(SingleCellExperiment)库(“BiocFileCache ") }) ## ----------------------------------------------------------------------------- url。base <- "https://github.com/goncalves-lab/waddR-data/blob/master/data/" sce.blood.url <- paste0(url。Base, "data_blood.rda?raw=true") sce.decidua.url <- paste0(url。base, "data_decidua.rda?raw=true") getCachedPath <- function(url, rname){bfc <- BiocFileCache() # fire up cache res <- bfcquery(bfc, url, field="fpath", exact= true) if (bfccount(res) == 0L) cachedFilePath <- bfcadd(bfc, rname=rname, fpath=url) else cachedFilePath <- bfcpath(bfc, res[["rid"]]) cachedFilePath} load(getCachedPath(cfc, res[["rid"]])url,“data_blood”)负载(getCachedPath (data_decidua sce.decidua.url。 ")) ## ----------------------------------------------------------------------------- set.seed (28) randgenes < -样本(rownames (data_blood), 1000年,取代= FALSE)预计。血< - data_blood [randgenes,] sce.decidua < - data_decidua [randgenes , ] ## ----------------------------------------------------------------------------- 化验(sce.blood)化验(sce.decidua ) ## ----------------------------------------------------------------------------- res < wasserstein.sc (sce。血液、sce.decidua方法= " t ", permnum = 1000,种子= 24)头(res, n = 10 l) # #——session-info ------------------------------------------------------------- sessionInfo ()