/code>,
where the same intercept <code>alpha</code> is used for all observations. The
number of rows of the independent variable matrix <code>x</code> needs to match
the size of the dependent variable vector <code>y</code> and the number of
columns of <code>x</code> needs to match the size of the coefficient vector <code>beta</code>.
If <code>x</code> and <code>y</code> are data (not parameters) this function can be executed on a GPU.</p>
<!-- real; poisson_log_glm_lpmf; (int[] y | matrix x, vector alpha, v