Fitting spatial regressions to large datasets using unilateral approximations