Sorry
should read:
ip - function(x) 1240 * x^(-0.88)
Rainer M Krug wrote:
Hi
System:
Linux, SuSE 10
R 2.3.0
I try to do the following
ip - function(x) 1240*dip(x, 0.88)
iip - function(x) integrate(ip, x - 0.045, x + 0.045)$value
f - integrate(iip, 10, 100)
Error in
to do it.
Med venlig hilsen
Frede Aakmann Tøgersen
-Oprindelig meddelelse-
Fra: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] På vegne af Rainer M Krug
Sendt: 22. juni 2006 14:03
Til: [EMAIL PROTECTED]
Cc: R help list
Emne: [R] problem with integrate() - correction
Sorry
:[EMAIL PROTECTED] På vegne af Rainer M Krug
Sendt: 22. juni 2006 14:03
Til: [EMAIL PROTECTED]
Cc: R help list
Emne: [R] problem with integrate() - correction
Sorry
should read:
ip - function(x) 1240 * x^(-0.88)
Rainer M Krug wrote:
Hi
System:
Linux, SuSE 10
R 2.3.0
I try to do
The background: I'm trying to fit a Poisson-lognormal distrbutuion to
some data. This is a way of modelling species abundances:
N ~ Pois(lam)
log(lam) ~ N(mu, sigma2)
The number of individuals are Poisson distributed with an abundance
drawn from a log-normal distrbution.
To fit this to data, I
Have you done a search of www.r-project.org - search - R site
search for hermite quadrature? I just got 11 hits on this, the
second of which referred to gauss.quad {statmod}. This said, See
also ... gauss.quad.prob {statmod}. I haven't tried it, but it looks
like what you want.
On Tue, 2 Mar 2004, Anon. wrote:
The background: I'm trying to fit a Poisson-lognormal distrbutuion to
some data. This is a way of modelling species abundances:
N ~ Pois(lam)
log(lam) ~ N(mu, sigma2)
The number of individuals are Poisson distributed with an abundance
drawn from a
Thomas Lumley wrote:
On Tue, 2 Mar 2004, Anon. wrote:
The background: I'm trying to fit a Poisson-lognormal distrbutuion to
some data. This is a way of modelling species abundances:
N ~ Pois(lam)
log(lam) ~ N(mu, sigma2)
The number of individuals are Poisson distributed with an abundance
drawn
On Tue, 2 Mar 2004, Anon. wrote:
Thomas Lumley wrote:
The help page for integrate() says
When integrating over infinite intervals do so explicitly, rather
than just using a large number as the endpoint. This increases
the chance of a correct answer - any function whose