空间聚类算法BayesSpace
空间基因表达技术能够在保留空间背景信息的同时，全面测量转录组图谱。然而，现有的分析方法并没有解决技术分辨率有限或有效利用空间信息的问题。
BayesSpace是一种完全贝叶斯统计方法，用于聚类分析来增强空间转录组数据的分辨率，并且可以通过 SingleCellExperiment 对象无缝集成到当前转录组的分析工作流程中。BayesSpace利用多元t分布对基因表达的低维表示进行建模，然后通过Potts模型合并空间信息，这使得相邻点属于同一簇。这种方法借鉴了用于图像分析的成熟的空间统计方法。有效地利用空间信息可以极大地提高空间数据聚类的性能。
readVisium()读入的是Spatial Ranger的输出数据，格式参考： https://support.10xgenomics.com/spatial-gene-expression/software/pipelines/latest/output/overview
BayesSpace需要很少的数据预处理，在这里这一步封装在了spatialPreprocess()这个函数。
spatialPreprocess() 对count矩阵进行了log标准化( log.normalize=T)，并对top n.HVGs 高变基因进行了PCA处理，保留了top n.PCs 个主成分。此外，spatial sequencing platform被作为metadata的信息添加到 SingleCellExperiment 对象用于下游分析。如果不想进行PCA这一步，在进行 spatialPreprocess() 的时候可以设置 skip.PCA=TRUE 。
我们使用 qTune() 和 qPlot() 功能来帮助选择q值，也就是 用于后续分析的cluster数 。（q值直接决定了后面聚出多少个cluster）
spatialCluster() 功能对spots进行聚类，并向 SingleCellExperiment 对象增加预测的cluster标签。通常，我们建议至少10,000次迭代（ nrep=10000 ），但是这个时间通常会很长，所以为了演示，这里只做1000次迭代。（注意：为了结果的一致性，这一步需要设置随机种子。）
在这一步之后，mclust初始化( cluster.init )和贝叶斯cluster分配( spatial.cluster )的结果都已经得到并储存在SingleCellExperiment的 colData 中。
clusterPlot()的结果是基于ggplot2的，因此可以使用ggplot2的语法对结果图进行优化。
spatialEnhance() 功能将增强主成分的分辨率，并将这些PC以及subspot分辨率的预测聚类标签添加到新的SingleCellExperiment对象。和前面的 spatialCluster 一样，为了演示这里的迭代也是使用了1000，但是建议使用10000或更大的迭代值。
增强的 SingleCellExperiment 包括了原始 sce 中父类spot的指标( spot.idx )，也包含了subspot的指标( subspot.idx )。它将offsets添加在原始点坐标上，并提供了增强的聚类坐标( spatial.cluster )。
可视化
BayesSpace对基因表达矩阵的主成分进行运算，因此 spatialEnhanced() 计算的是增强分辨率的PC向量。增强的基因表达并不是直接计算的，而是通过回归算法计算的。对于每个基因，使用每个点的PC向量训练了模型来预测spot-level的基因表达，并使用拟合的模型来预测来自subspot PC的subspot表达。
基因表达增强是使用 enhancedFeature() 函数实现的。BayesSpace默认使用 xgboost 预测基因表达，但也可以通过模型参数使用线性和 Dirichlet狄利克雷回归。使用xgboost时，我们建议将nrounds参数设置为0来自动调整它，尽管这大大增加了运行时间（比预先指定的nrounds慢了约4倍）。
EnnhancedFeatures() 可用于计算所有基因或感兴趣基因集的subspot表达水平。在这里演示了四个常见基因。
默认情况下，计算的是log标准化的expression( logcounts(sce) )，尽管也可以指定其他assays或者arbitrary features矩阵。
Diagnostic measures from each predictive model, such as rmse when using xgboost, are added to the rowData of the enhanced dataset.
参考： https://edward130603.github.io/BayesSpace/articles/BayesSpace.html
扩展阅读：免费iq测试题智力测试 ... 过客软件园qq群号2023 ... 过客软件库qq群232581034 ... citespace聚类太拥挤 ... citespace聚类分析图 ... 填料最小承载比cbr ... 过客园软件网址进入 ... 过客软件库合集入口 ... citespace聚类轮廓值 ...
