元件之间我在R.新用户如何检查该矢量A
的元素是矢量B
的例如元件之间:检查哪个一个矢量的元素是另一个的R中
A = c(1.1, 3.2, 5, 8.5, 4.6, 2.2)
B = c(1, 2, 3)
表示A
的哪些元素在1和2之间,那么A
的哪些元素在2和3之间?
元件之间我在R.新用户如何检查该矢量A
的元素是矢量B
的例如元件之间:检查哪个一个矢量的元素是另一个的R中
A = c(1.1, 3.2, 5, 8.5, 4.6, 2.2)
B = c(1, 2, 3)
表示A
的哪些元素在1和2之间,那么A
的哪些元素在2和3之间?
你想值1的findInterval
findInterval(A,B)
[1] 1 3 3 3 3 2
一些变体表示为1和2(B中的最低和下一个最低值) 值2表明它是在2和3之间
因此,要找到哪些是
which(findInterval(A,B) %in% seq_len(length(unique(B))-1))
# [1] 1 6
之间,并从A
A[findInterval(A,B) %in% seq_len(length(unique(B))-1)]
# [1] 1.1 2.2
您也可以使用cut
,这将创建一个因子。
与split
结合这会给
split(A, cut(A,B),drop=FALSE)
$`(1,2]`
[1] 1.1
$`(2,3]`
[1] 2.2
如果我理解正确的话,这是一种可能性:
A = c(1.1, 3.2, 5, 8.5, 4.6, 2.2)
B = c(1, 2, 3,4,10)
B1 <- head(B, -1)
B2 <- tail(B, -1)
outs <- list()
for(i in seq_along(B1)) {
outs[[i]] <- A[B1[i] < A & A < B2[i]]
}
names(outs) <- paste(B1, " & ", B2)
## > outs
## $`1 & 2`
## [1] 1.1
##
## $`2 & 3`
## [1] 2.2
##
## $`3 & 4`
## [1] 3.2
##
## $`4 & 10`
## [1] 5.0 8.5 4.6
你刚才在这里重新创建了'cut'和'split' ...... :) :) – mnel 2013-04-04 04:50:26
我意识到看到你的回应后:)我的可能不是一样快。 – 2013-04-04 04:51:49
我假定这将是连续的数字B中为您将检查
result.pair<-matrix(rep(integer(), 3), ncol = 3)
colnames(result.pair)<-c("B1", "A", "B2")
for(i in 1:(length(B)))
{
for(j in 1:(length(A)))
{
if ((B[i] <= A[j]) & (B[i+1] >= A[j]))
{
result.pair<-rbind(result.pair, c(B[i], A[j], B[i+1]))
}
}
}
result.pair
试试这个:
data.frame(A,cut(A,B))
对于A中的每一个观察,它都会告诉你它之间的哪一对B观测值。
像这样:
> data.frame(A,cut(A,B))
A cut.A..B.
1 1.1 (1,2]
2 3.2 <NA>
3 5.0 <NA>
4 8.5 <NA>
5 4.6 <NA>
6 2.2 (2,3]
NA意味着它不是两个B观测之间。
相关搜索:
data.frame(A,cut(A,c(-Inf,B,Inf)))
的期望是什么输出?一个列表? – Nishanth 2013-04-04 05:17:24