From: tintin_et_milou <badgones65_at_hotmail.com>

Date: Mon, 3 Dec 2007 06:19:36 -0800 (PST)

Date: Mon, 3 Dec 2007 06:19:36 -0800 (PST)

Thanks for your help, but there is some more problem. The two vectors have not the same length so there is a problem with cbind. I give you an example. My first vector is

2035.836 652.9494 2035.939 664.5841 2036.043 696.0554 2036.146 719.8969 2036.250 750.7660 2036.767 816.5243 2036.870 806.8539 2036.974 774.2397 2037.491 777.2780 2039.147 589.9075 2042.978 807.7167 2043.082 820.9365 2043.289 849.4942 2043.393 883.8975 2043.495 900.9681 2043.600 922.3238 2043.704 956.5985 2043.911 978.9377 2044.015 969.1999

and the second one is:

> f[1:9,]

M.Z Intensity Echantillon Position 1802.809 1064.1210 1 A1 1865.615 8799.4880 1 A1 1896.426 1667.5908 1 A1 2001.064 515.6214 1 A1 2012.016 837.5599 1 A1 2021.589 4373.6364 1 A1 2028.425 832.6896 1 A1 2036.663 0.0000 1 A1 2043.497 0.0000 1 A1

The two vectors have not the same number of observations and in fact i would like to replace all the intensity of zero of the second vector by the value of the intensity of the nearest m/Z of the first vector. In this case, the value 816.5243 for the m/Z 2036.663 of the second vector and 900.9681 for the intensity corresponding at m/Z=2043.497.

Your method give me the intensity of the farest m/Z.

Ben Bolker wrote:

*>
**>
**>
*

> ## construct sample data

*> m1 =
**> matrix(c(1000.235,1000.356,
**> 125,126.5),ncol=2,
**> dimnames=list(NULL,c("mZ","I")))
**>
**> y = c(995.547,1000.320)
**>
**> ## compute distances, set diagonal to infinity
**> d = as.matrix(dist(cbind(m1[,1],y)))
**> diag(d) <- Inf
**>
**> ## find minimum distances, extract values
**> cbind(y,m1[apply(d,2,which.min),2])
**>
**> Ben Bolker
**>
**>
*

-- View this message in context: http://www.nabble.com/comparison-of-two-vectors-tf4936213.html#a14130713 Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help_at_r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.Received on Mon 03 Dec 2007 - 20:11:31 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Mon 03 Dec 2007 - 23:30:17 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*