Hi all,
 I am running a deSolve based simulation (epidemiological & demographic model).
Some demographic parameters like birth follow a normal distribution to better reflect
the natural dynamic. Basically implemented as follows:

## B =  birth
## h =  total numbers of juveniles
## my = mortality rate

model <- function (t, y, p) {
  with(as.list(c(y, p)),{

      B <- birth(t)

      dN = B*h - my*y[1]
    return(list(c(dN)))
  })
}

birth <- approxfun(x = 0:100, y = dnorm(0:100, 50, 25))

out <- ode(y = 100, times = 0:100, parms = c(my = 4e-2, h = 5), func = model)
plot(out, ylab = "Population size", xlab = "Time”)

My question is rather post-modeling related. I am struggling to explain this in the
method section of a manuscript. Basically I would like to receive a recommendation
how to describe this process formally and mathematically right i.e. birth follows
a normal distribution derived by multiplying the number of hatchlings with a density
function of a normal distribution …. ?!

Thanks,
Simeon




	[[alternative HTML version deleted]]

