fusion function

Fit a spatial fusion model