Displaying 3 results from an estimated 3 matches for "dapolono".
2008 Feb 15
0
Poisson-lognormal probability calculations
...ction value
-------------
I tidied up a little (to my eyes only, no offence intended to the estimable
authors of VGAM) and avoided the gamma overflow by using logs - OK so far,
it agrees almost perfectly with VGAM::dpolono for x up to 170, and now
extends the range up to 173 (wow!!).
-------------
dApolono <-
function (x, meanlog = 0, sdlog = 1, ...) {
require(stats)
integrand <- function(t, x, meanlog, sdlog)
exp(-t+x*log(t)-log(sdlog*t*sqrt(2*pi))-0.5*((log(t)-meanlog)/sdlog)^2)
mapply(function(x, meanlog, sdlog, ...){
temp <- try(
integrate(f = integrand, lower...
2008 Feb 15
0
Behaviour of integrate (was 'Poisson-lognormal probability calcul ations')
...; I tidied up a little (to my eyes only, no offence intended to the
> estimable authors of VGAM) and avoided the gamma overflow by using logs
> - OK so far, it agrees almost perfectly with VGAM::dpolono for x up to
> 170, and now extends the range up to 173 (wow!!).
> -------------
> dApolono <-
> function (x, meanlog = 0, sdlog = 1, ...) {
> require(stats)
> integrand <- function(t, x, meanlog, sdlog) exp(-t+x*log(t)-
> log(sdlog*t*sqrt(2*pi))-0.5*((log(t)-meanlog)/sdlog)^2)
> mapply(function(x, meanlog, sdlog, ...){
> temp <- try(
>...
2008 Feb 18
0
Solved (??) Behaviour of integrate (was 'Poisson-lognormal probab ility calculations')
...only, no offence intended to the
> > estimable authors of VGAM) and avoided the gamma overflow by using
> logs
> > - OK so far, it agrees almost perfectly with VGAM::dpolono for x up
> to
> > 170, and now extends the range up to 173 (wow!!).
> > -------------
> > dApolono <-
> > function (x, meanlog = 0, sdlog = 1, ...) {
> > require(stats)
> > integrand <- function(t, x, meanlog, sdlog) exp(-t+x*log(t)-
> > log(sdlog*t*sqrt(2*pi))-0.5*((log(t)-meanlog)/sdlog)^2)
> > mapply(function(x, meanlog, sdlog, ...){
> >...