• predict(..., newdata) can take a matrix now