Michael originally suggested ?outer. I think that was enough in this case (no need for mv or sapply):
wlpk3 <- outer(k3,wl,'+') ln.phiDIC <- log(k1)+k2/wlpk3 phiDIC<- t(exp(ln.phiDIC)) colnames(phiDIC)<- stations str(phiDIC) Cheers On Thu, Apr 12, 2012 at 8:58 PM, R. Michael Weylandt <michael.weyla...@gmail.com> wrote: > Strange, that isn't the error I get: > >> mouter(wl, k1, k2, k3, FUN = function(w, k1, k2, k3) k1 *exp(k2 / (w + k3))) > Error in FUN(X, Y, ...) : argument "k2" is missing, with no default > > Still, it's a problem with my mouter() function which was only tested > on binary operators (and then only really to teach myself to use > Recall()). This should work better: > > `mouter` <- function(..., FUN){ > dotArgs <-as.list(do.call("expand.grid", list(..., KEEP.OUT.ATTRS= F))) > names(dotArgs) <- names(formals(FUN)) > ans <- do.call(FUN, dotArgs) > dim(ans) <- sapply(list(...), length) > ans > } > > This only works when the function can take all the elements at a time > though; I'm sure some fooling around could combine them nicely...it > seems to test right on your data, but I haven't checked it more > generally. > > Michael > > On Thu, Apr 12, 2012 at 3:29 PM, Filoche <pmassico...@hotmail.com> wrote: >> Hi and thank you for your time. >> >> I got this error when trying your function. >> >> mouter(wl, k1, k2, k3, FUN = function(w, k1, k2, k3) k1 *exp(k2 / (w + k3))) >> "Error in k3/(w + k3) : 'k3' is missing" >> >> Regards, >> Phil >> >> -- >> View this message in context: >> http://r.789695.n4.nabble.com/Help-with-vectorization-tp4552638p4552833.html >> Sent from the R help mailing list archive at Nabble.com. >> >> ______________________________________________ >> 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. > > ______________________________________________ > 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. ______________________________________________ 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.