REAL alpha ~ user_gamma(0.25, 0.5)
REAL tau ~ continuous_uniform(0.5, 10)
REAL gamma ~ continuous_uniform(-1.0, 1.0)

ENTITY species, "models/birds/species.data" 
{
    REAL a ~ user_normal(alpha, tau)
    INTEGER waterbird(1)
}

ENTITY cell, "models/birds/gridcell.data", "models/birds/neighbours.data"
{
    REAL waterpercentage(4)
    INTEGER researchgrade(7)
    INTEGER centeredcoordinate(3)
}

ENTITY observations, "models/birds/obs.data", combines (cell, species)
{
    REAL q = SUM(&x)
    REAL p  = EXP(waterbird*a*waterpercentage + gamma*centeredcoordinate + 0.25 * q)/(1 +  EXP(waterbird*a*waterpercentage + gamma*centeredcoordinate + 0.25 * q))
    # 0.25 could be any other constant. defines spatial weighting.

    INTEGER x ~ user_bernoulli(p)
    INTEGER observation(*) ~ user_defined_points2(x, researchgrade)
}