Overview
Implementation
{% \mu_i = exp[\vec{x}_i^T \vec{\beta}] %}
let beta = pl.irls(data,'num_awards', ['progAcademic','progVocational','math']);
let beta = pl.irls(data,'num_awards', ['progAcademic','progVocational','math']);