def growth(self, x, y, t):
a, b, c, sd = self.growth_params
m = b * c * x**(c-1.0) / ( b + x**c/a )**2
return dnorm(y, mu=x+m, sd=sd)
def growth(self, x, y, t):
a0, a1, b0, b1 = self.growth_params
m = a0 + a1 * x
v = b0 + b1 * x
return dnorm(y, mu=m, sd=sqrt(v))
def growth(self, x, y, t):
pptWin2 = self.covariat('pptWin2')
TmeanSum1 = self.covariat('TmeanSum1')
a0, a1, a2, a3 = self.growth_params
m = a0 + a1 * x + a2 * pptWin2 + a3 * TmeanSum1
v = 1.79
return dnorm(y, mu=m, sd=sqrt(v))
def n0(self, x):
return dnorm(x, mu=50.0, sd=sqrt(10.0))
def kernel(self, x, y, t, **kwargs):
self.increment_count(x, y)
return exp(-self._theta * x**2) * dnorm(y, mu=x+self._mu, sd=sqrt(self._ss))
def n0(self, x):
return dnorm(x, mu=self.mu0, sd=sqrt(self.ss0))