在R

时间:2018-04-07 02:11:59

标签: r list iteration montecarlo

我正在处理RNA数据中的一些计算特征提取问题,我发现自己无法处理这个问题:

我有n个序列(例如两个),我从中获得了迭代统计量i次(进行蒙特卡罗迭代以分析获得的统计数据与原始数据的比较)。

示例:

假设我们迭代10次

n <- 10

我得到了一个包含所有迭代的20个值的向量,但是这个向量对应于两个不同的序列,所以我必须将这个向量分成两个相等的部分(每个序列的迭代顺序为1:10 - 1:10)

MFEit <- c(10, 12, 34, 32, 12 .....)  ## vector of length 20
MFEit.split <- split(MFEit, ceiling(MFEit.along/n5))

这会生成一个包含两个项目的列表,每个项目包含10个值,名为$ 1和$ 2

另一方面,我有两个值的向量,这两个值是原始统计数据,每个值对应于每个原始序列

MFE <- c(25, 15)

我想要做的是知道列表MFEit.split中第一项的多少值等于或小于MFE的第一个值,并且迭代地知道列表MFEit中第二项的多少值。 split,等于或小于MFE的第二个值,依此类推,只要我有两个以上的值或项。

我知道怎么一个接一个地说,

R <- length(subset(MFEit.split$`1`, MFEit.split$`1`<=MFE[1]))
R <- length(subset(MFEit.split$`2`, MFEit.split$`1`<=MFE[2]))

但是......如何将它包含在一个循环中,以便我可以迭代地进行每次比较,无论我有多少MFE值或列表中的项目?

所需的输出将是一个名为R的向量,其中n个值对应于每个比较。

任何帮助?...

0 个答案:

没有答案