merge (chondro [1:10,, 600], chondro [5:15,, 600], by = c("x", "y"))$.
tmp <- merge (chondro [1:10,, 610], chondro [5:15,, 610],
by = c("x", "y"), all = TRUE)
tmp$.
wl (tmp)
stopifnot (nrow (merge (chondro [1:10], chondro [5:15], all = TRUE)) == 15)
stopifnot (nrow (merge (chondro [1:10], chondro [5:15])) == 6)
## remove duplicated wavelengths:
approxfun <- function (y, wl, new.wl){
approx (wl, y, new.wl, method = "constant",
ties = function (x) mean (x, na.rm = TRUE)
)$y
}
merged <- merge (chondro [1:7,, 610 ~ 620], chondro [5:10,, 615 ~ 625], all = TRUE)
merged$.
merged <- apply (merged, 1, approxfun,
wl = wl (merged), new.wl = unique (wl (merged)),
new.wavelength = "new.wl")
merged$.
Run the code above in your browser using DataLab