fitting spatial regression to large datasets using unilateral approximations