how to pass a function to a function inside a function

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120925/8e3d3a49/attachment.pl>
Hello,

Your function creates the vector 'x', then creates the function fs(), 
then returns the function just created. It _never_ reaches the line

     intgrt <- function(lower) {

and therefore never creates this last function. Now supposed it did. In 
the very end you make the same mistake, you simply return the function 
without ever calling it.

Hope this helps,

Rui Barradas

Em 25-09-2012 13:28, Al Ehan escreveu:
Hi,

I'm trying to compile two functions into one function. the first funtion is
called 'fs' which is self-made function, another function is from the
built-in 'integration' function that is copy-paste-edited. If built
separatey, these functions work well. However that is not the case if
combines together, where certainly I made mistake somewhere when
constructing the code.

can someone help point out the mistakes for me please?

Thanks!
Aehan

____________________

intgfun<- function(a,b,c,mu,alpha,xi,upper, ...,subdivisions = 100, rel.tol
= .Machine$double.eps^0.25,
     abs.tol = rel.tol, stop.on.error = TRUE, keep.xy = FALSE,
     aux = NULL)

{
x<-vector()
fs<-function(x){
temp<-(-xi^(-1))*log(1-xi*(x-mu)/alpha)

(((alpha^(-1))*(exp((-(1-xi)*temp)-exp(-temp))))*((a*(x-c)^0.5)+(b*(x-c))))
}
return(fs)

intgrt<-function(lower)
{
call(fs)
     fs <- match.fun(f)
     ff <- function(x) f(x, ...)
     limit <- as.integer(subdivisions)
     if (limit < 1 || (abs.tol <= 0 && rel.tol < max(50 *
.Machine$double.eps,
         5e-29)))
         stop("invalid parameter values")
     if (is.finite(lower) && is.finite(upper)) {
         wk <- .External("call_dqags", ff, rho = environment(),
             as.double(lower), as.double(upper), as.double(abs.tol),
             as.double(rel.tol), limit = limit, PACKAGE = "base")
     }
     else {
         if (is.na(lower) || is.na(upper))
             stop("a limit is missing")
         if (is.finite(lower)) {
             inf <- 1
             bound <- lower
         }
         else if (is.finite(upper)) {
             inf <- -1
             bound <- upper
         }
         else {
             inf <- 2
             bound <- 0
         }
         wk <- .External("call_dqagi", ff, rho = environment(),
             as.double(bound), as.integer(inf), as.double(abs.tol),
             as.double(rel.tol), limit = limit, PACKAGE = "base")
     }
     res <- wk[c("value", "abs.error", "subdivisions")]
     res$message <- switch(wk$ierr + 1, "OK", "maximum number of
subdivisions reached",
         "roundoff error was detected", "extremely bad integrand behaviour",
         "roundoff error is detected in the extrapolation table",
         "the integral is probably divergent", "the input is invalid")
     if (wk$ierr == 6 || (wk$ierr > 0 && stop.on.error))
         stop(res$message)
     res$call <- match.call()
     class(res) <- "integrate"
     res
}
return (intgrt)}

	[[alternative HTML version deleted]]

______________________________________________
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.