[R] Index out SNP position
David Winsemius
dwinsemius at comcast.net
Fri Jan 4 01:59:02 CET 2013
On Jan 3, 2013, at 1:54 PM, JiangZhengyu wrote:
>
>
>
>
>
>
> Dear R experts,
>
> I have 2 matix: A& B. I am trying to index B against A - (1) find out B rows that fall between the col 1 and 2 of A& put them into a new vector SNP.I made code as below, but I cannot think of a right way to do it. Could anyone help me with the code? Thanks,Jiang----
>
> A <- matrix(c(35838396,35838674,36003908,36004090,36150188,36151202,35838584,35838674,36003908,36003992), ncol = 2)
> B <- matrix(c(36003918,35838399,35838589,36262559),ncol = 1) nr=nrow(A)
> rn=nrow(B) for (i in 1:nr)
> {
> for (j in 1:rn){if (B[i,1]<=A[j,1] && B[i,1]>=A[j,2]){SNP[i]=B[i,1]}}
> }
> sapply(B, function(x) apply(A, 1, function(two) x %in% two[1]:two[2]))
[,1] [,2] [,3] [,4]
[1,] TRUE TRUE TRUE FALSE
[2,] FALSE FALSE TRUE FALSE
[3,] FALSE FALSE FALSE FALSE
[4,] TRUE FALSE FALSE FALSE
[5,] FALSE FALSE FALSE FALSE
So the first and third B-locations are in the range of two of the rows, the second-B in one range and the third is in none of them. There is also a bioconductor package called `IRanges` that will undoubtedly be more efficient. (This works because the problem is of necessity dealing with integers.)
--
David Winsemius
Alameda, CA, USA
More information about the R-help
mailing list