REAL alpha ~ user_normal(0, 0.00001) REAL sigma ~ user_gamma(1.00001, 1.0) ENTITY county, "models/cancer/lips.data", "models/cancer/lips.neighbours" { REAL x(3) REAL expected(2) REAL n = COUNT(&b) REAL q = SUM(&b) REAL t = q/n REAL b ~ user_normal(t, n) REAL mu = EXP(LOG(expected) + alpha * x / 10 + sigma * b) INTEGER obs(1) ~ poisson(mu) } # x = the persentage of farmland # expected = expected count of lip cancer patients # b = something spatial # mu = the real count of lip cancer patients # obs = observed count of lip cancer patients