Just a curiosity question:
In the documentation for the nlm procedure
a find this example of defining a function
with a gradient attribute:
-----------
f <- function(x, a)
{
res <- sum((x-a)^2)
attr(res, "gradient") <- 2*(x-a)
res
}
-----------
I get the gradient with
attr(f(3,2),"gradient")
but how do I get the function value it self?
Gerrit