double logistic(double b, void *data);
double poisson(double b, void *data);

struct glmparm {
   int idx;
   Vector *y, *beta;
   Matrix *X;
};