# [R] Index out SNP position

David L Carlson dcarlson at tamu.edu
Thu Jan 3 23:22:57 CET 2013

```Something like this?

> indx <- sapply(1:nrow(B), function(i) any(B[i]>A[,1] & B[i]<A[,2]))
> SNP <- B[indx]
> SNP
[1] 36003918 35838399 35838589

----------------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77843-4352

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of JiangZhengyu
> Sent: Thursday, January 03, 2013 3:55 PM
> To: r-help at r-project.org
> Subject: [R] Index out SNP position
>
>
>
>
>
>
>
> 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]}}
> }
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help