我有一份多组学数据清单,其中包括三种基因/蛋白质/细胞的矩阵,以及列中的个体样本。这三个矩阵共享相同的样本。我使用多个CO惯性分析将每个样本投射到CIA空间(每个样本由基于转录、蛋白质模拟和细胞数据的三个点表示)。现在,我在我原来的空间结构中有了一个新的观点,我想把它投射到中央情报局的坐标系中。
我怎么才能得到新的中情局坐标?
为了重现性,我以NCI60数据为例。
library(omicade4)
data(NCI60_4arrays)
NCI60_4arrays.sub <- lapply(NCI60_4arrays, function(m) { return(m[,1:
(ncol(m)-1)])})
mcoin <- mcia(NCI60_4arrays.sub, cia.nf=4)
用于提取CIA空间中的样本坐标:
sample.cia.coordinates <- mcoin$mcoa$Tli
我怎么能投射一个新的数据点,并得到它的中央情报局坐标?
例如,一个新的数据点:
new.point <- lapply(NCI60_4arrays, function(m) { return(m[,ncol(m)])})
发布于 2019-02-06 15:28:15
当您想要理解分析背后的计算时,一种选择是创建一个低容量的示例分析,例如:
sample <- lapply(NCI60_4arrays[1:2], function(m) { return(m[1:7,1:5])})
sample_mcoa <- mcia(sample, cia.nf = 2)
然后试着理解你是如何传递的:
至:
它可能是几个矩阵产生的,其中一个与sample_mcoa$mcoa$axis
。我不知道这是否对你有帮助。
https://stackoverflow.com/questions/54555479
复制相似问题