R help - Jul 2004 - Integration with "adapt"

Hi all,

I need to calculate  a multidimensional integration on R. I am using the
command "adapt" (from library adapt), although sometimes I get the
following error message:

Ifail=2, lenwrk was too small. -- fix adapt() !
 Check the returned relerr! in: adapt(3, linf, lsup, functn = Integrando1)

I guess it happens because the domain of integration is too small,
although I tried a change of variables to avoid this problem and it didn?t
help. The command adapt calls a fortran routine, but I don?t know fortran
enough to fix the problem.
Can someone help me?
Thanks a lot
Rodrigo Drummond


____________________________________________
Rodrigo D. Drummond
Laboratorio Genoma Funcional
Centro de Biologia Molecular e Eng. Genetica
Universidade Estadual de Campinas
Caixa Postal 6010
13083-875 - Campinas - SP - Brasil
Tel: xx-19-3788-1119 Fax: xx-19-3788-1089

>>>>> "Rodrigo" == Rodrigo Drummond <rduarte at
unicamp.br>
>>>>>     on Tue, 27 Jul 2004 15:31:07 -0300 (BRT) writes:
    Rodrigo> Hi all, I need to calculate a multidimensional
    Rodrigo> integration on R. I am using the command
    Rodrigo> "adapt" (from library adapt), although

it's a "package", not a library.

    Rodrigo> sometimes I get the following error message:

    Rodrigo> Ifail=2, lenwrk was too small. -- fix adapt() !
    Rodrigo> Check the returned relerr! in: adapt(3, linf, lsup,
    Rodrigo> functn = Integrando1)

If you could give as a *reproducible* example,

we (the adapt authors) would have chance to do what the above
message says, namely "fix adapt()" ..


    Rodrigo> I guess it happens because the domain of
    Rodrigo> integration is too small,

maybe, maybe not. We need an example we can reproduce, see above.

    Rodrigo> although I tried a change of variables to avoid
    Rodrigo> this problem and it didn?t help. The command adapt
    Rodrigo> calls a fortran routine, but I don?t know fortran
    Rodrigo> enough to fix the problem.

Martin Maechler, ETH Zurich