To incorporate bias correction, extra parameters are used to model smooth intensity variations. A field modelling the variation at element i is denoted i(), where is a vector of unknown parameters. Intensities from the ith cluster are assumed normally distributed, with mean k/i() and variance (k/i())2: