On Sat, Apr 2, 2011 at 7:59 PM, David Winsemius <dwinsem...@comcast.net> wrote: > > On Apr 2, 2011, at 10:20 PM, Joshua Wiley wrote: >> VA <- cbind(OS = c(.2, 0, 1, -.1, NA, 3, NA), >> OD = c(.3, -.1, .2, -.1, NA, 0, .1), >> Desired = c(0.225, -0.075, 0.4, -0.1, NA, 0.75, 0.1)) > > VA <- cbind(VA, 0.75*pmin(VA[, "OS"], VA[, "OD"], na.rm=TRUE) + > 0.25*pmax(VA[, "OS"], VA[, "OD"], na.rm=TRUE) ) > >> VA > OS OD Desired > [1,] 0.2 0.3 0.225 0.225 > [2,] 0.0 -0.1 -0.075 -0.075 > [3,] 1.0 0.2 0.400 0.400 > [4,] -0.1 -0.1 -0.100 -0.100 > [5,] NA NA NA NA > [6,] 3.0 0.0 0.750 0.750 > [7,] NA 0.1 0.100 0.100 >
That works wonderfully, and it is so elegant! Thank you, David. Josh > > David Winsemius, MD > West Hartford, CT > ______________________________________________ R-help@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.
