//@{ const ArrayXd probitLink::linkFun(const ArrayXd& mu) const {return mu.unaryExpr(probit<double>());}
const ArrayXd logitLink::muEta( const ArrayXd& eta) const {return eta.unaryExpr(logitmueta<double>());}
const ArrayXd logitLink::linkInv(const ArrayXd& eta) const {return eta.unaryExpr(logitinv<double>());}
const ArrayXd logLink::muEta( const ArrayXd& eta) const {return eta.unaryExpr(boundexp<double>());}
//@{ const ArrayXd cauchitLink::linkFun(const ArrayXd& mu) const {return mu.unaryExpr(cauchit<double>());}
const ArrayXd binomialDist::variance(const ArrayXd& mu) const {return mu.unaryExpr(x1mx<double>());}