fusion.dstan function

Fit a spatial fusion model using Stan