The SAR model for very large datasets: A reduced-rank approach