add basic GLM prediction method for new X matrix

Addresses #76.
Can this be `linkinv!(mm.rr.l, eta, eta)` to avoid allocating a second array?