Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----eval=FALSE---------------------------------------------------------------
# library("SC.MEB")
## ----eval=FALSE---------------------------------------------------------------
# library(mvtnorm)
# library(GiRaF)
# library(SingleCellExperiment)
# set.seed(100)
# G <- 4
# Bet <- 1
# KK <- 5
# p <- 15
# mu <- matrix(c( c(-6, rep(-1.5, 14)),
# rep(0, 15),
# c(6, rep(1.5, 14)),
# c(rep(-1.5, 7), rep(1.5, 7), 6),
# c(rep(1.5, 7), rep(-1.5, 7), -6)), ncol = KK)
# height <- 70
# width <- 70
# n <- height * width # # of cell in each indviduals
## ----eval=FALSE---------------------------------------------------------------
# X <- sampler.mrf(iter = n, sampler = "Gibbs", h = height, w = width, ncolors = KK, nei = G, param = Bet,initialise = FALSE, view = TRUE)
# x <- c(X) + 1
# y <- matrix(0, nrow = n, ncol = p)
#
# for(i in 1:n) { # cell
# mu_i <- mu[, x[i]]
# Sigma_i <- ((x[i]==1)*2 + (x[i]==2)*2.5 + (x[i]==3)*3 +
# (x[i]==4)*3.5 + (x[i]==5)*4)*diag(1, p)*0.3
# y[i, ] <- rmvnorm(1, mu_i, Sigma_i)
# }
#
# pos <- cbind(rep(1:height, width), rep(1:height, each=width))
## ----eval=FALSE---------------------------------------------------------------
# # -------------------------------------------------
# # make SC-MEB metadata used in SC-MEB
# counts <- t(y)
# rownames(counts) <- paste0("gene_", seq_len(p))
# colnames(counts) <- paste0("spot_", seq_len(n))
#
# ## Make array coordinates - filled rectangle
# cdata <- list()
# nrow <- height; ncol <- width
# cdata$row <- rep(seq_len(nrow), each=ncol)
# cdata$col <- rep(seq_len(ncol), nrow)
# cdata <- as.data.frame(do.call(cbind, cdata))
# ## Scale and jitter image coordinates
# #scale.factor <- rnorm(1, 8); n_spots <- n
# #cdata$imagerow <- scale.factor * cdata$row + rnorm(n_spots)
# #cdata$imagecol <- scale.factor * cdata$col + rnorm(n_spots)
# cdata$imagerow <- cdata$row
# cdata$imagecol <- cdata$col
# ## Make SCE
# ## note: scater::runPCA throws warning on our small sim data, so use prcomp
# sce <- SingleCellExperiment(assays=list(counts=counts), colData=cdata)
# reducedDim(sce, "PCA") <- y
# # sce$spatial.cluster <- floor(runif(ncol(sce), 1, 3))
#
# metadata(sce)$SCMEB.data <- list()
# metadata(sce)$SCMEB.data$platform <- "ST"
# metadata(sce)$SCMEB.data$is.enhanced <- FALSE
## ----eval=FALSE---------------------------------------------------------------
# platform = "ST"
# beta_grid = seq(0,4,0.2)
# K_set= 2:10
# parallel=TRUE
# num_core = 3
# PX = TRUE
# maxIter_ICM = 10
# maxIter = 50
## ----eval=FALSE---------------------------------------------------------------
# Adj_sp <- getneighborhood_fast(as.matrix(pos), cutoff = 1.2)
## ----eval=FALSE---------------------------------------------------------------
# Adj_sp <- find_neighbors2(sce, platform = platform)
# Adj_sp[1:10,1:10]
## ----eval=FALSE---------------------------------------------------------------
# fit = SC.MEB(y, Adj_sp, beta_grid = beta_grid, K_set= K_set, parallel=parallel, num_core = num_core, PX = PX, maxIter_ICM=maxIter_ICM, maxIter=maxIter)
# str(fit[,1])
## ----eval=FALSE---------------------------------------------------------------
# selectKPlot(fit, K_set = K_set, criterion = "BIC")
## ----eval=FALSE---------------------------------------------------------------
# selectKPlot(fit, K_set = K_set, criterion = "MBIC")
## ----eval=FALSE---------------------------------------------------------------
# out = selectK(fit, K_set = K_set, criterion = "BIC")
# ClusterPlot(out, pos)
## ----eval=FALSE---------------------------------------------------------------
# ClusterPlot(out, pos) +
# theme_bw() +
# xlab("Row") +
# ylab("Column") +
# labs(title="Spatial clustering")
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.