actually i figured it out: regRes[dim(regRes)[1]] should be regRes[dim(regRes)[1], ] -- View this message in context: http://r.789695.n4.nabble.com/pairwise-linear-regression-between-two-large-datasets-tp4527476p4527680.html Sent from the R help mailing list archive at Nabble.com.