我正在使用一个矩阵 (data1) 中的值来确定如何将第二个矩阵 (data2) 子集化为一个较小的矩阵 (foo)。但是我遇到了问题,因为我无法提前知道 foo 的大小。
这是一个 MWE(实际的数据结构有 100,000 多行)。
data1 <- rbind(c(102,250,'stim1'),c(477,839,'stim2'))
data2 <- rbind(c(99,'xx'),c(105,'yy'),c(230,'zz'),c(312,'aa'),c(587,'bb'),c(846,'cc'))
foo <- NULL
for(i in 1:nrow(data1))
{
foo[i,1] <- subset(data2,as.numeric(data2[,1]) > as.numeric(data1[i,1]) & as.numeric(data2[,1]) < as.numeric(data1[i,2]))
foo[i,2] <- rep.int(data1[i,3],nrow(subset(data2,as.numeric(data2[,1]) > as.numeric(data1[i,1]) & as.numeric(data2[,1]) < as.numeric(data1[i,2]))))
}
我希望 foo 看起来像这样:
105 yy stim1
230 zz stim1
587 bb stim2