Our aim was to examine the territorial dependence of risk for household insurances. Besides the classical risk factors such as type of wall, type of building, etc., we consider the location associated to each contract. A Markov random field model seems to be appropriate to describe the spatial effect. Basically there are two ways of fitting the model; we fit a GLM to the counts of claims with the classical risk factors and regarding their effects as fixed we fit the spatial model. Alternatively we can estimate the effects of all covariates (including location) jointly. Although this latter approach may seem to be more accurate, its high complexity and computational demands makes it unfeasible in our case. To overcome the disadvantages of the distinct estimation of the classical and the spatial risk factors proceed as follows: use first a GLM for the non-spatial covariates, and then fit the spatial model by MCMC. Refit next the GLM with keeping the obtained spatial effect fixed and afterwards refit the spatial model, too. Iterate this procedure several times. We achieve much better fit by performing eight iterations.