@staticmethod
def _model(nplayer: int, nmatch: int, minutes: jnp.array, y: jnp.array, alpha: jnp.array):
    theta = dist.Dirichlet(concentration=alpha)
    with numpyro.plate('nplayer', nplayer):
        dprobs = numpyro.sample('probs', theta)
    prob_score = numpyro.deterministic('prob_score', (dprobs[:, 0, None] * (minutes / 90.0)))
    prob_assist = numpyro.deterministic('prob_assist', (dprobs[:, 1, None] * (minutes / 90.0)))
    prob_neither = numpyro.deterministic('prob_neither', ((dprobs[:, 2, None] * (minutes / 90.0)) + (90.0 - minutes)))
    theta_mins = dist.Multinomial(probs=jnp.moveaxis(jnp.array([prob_score, prob_assist, prob_neither]), 0, 2))
    return numpyro.sample('obs', theta_mins, obs=y)
@staticmethod
def _model(nplayer: int, nmatch: int, minutes: jnp.array, y: jnp.array, alpha: jnp.array):
    theta = dist.Dirichlet(concentration=alpha)
    with numpyro.plate('nplayer', nplayer):
        dprobs = numpyro.sample('probs', theta)
    prob_score = numpyro.deterministic('prob_score', (dprobs[:, 0, None] * (minutes / 90.0)))
    prob_assist = numpyro.deterministic('prob_assist', (dprobs[:, 1, None] * (minutes / 90.0)))
    prob_neither = numpyro.deterministic('prob_neither', ((dprobs[:, 2, None] * (minutes / 90.0)) + (90.0 - minutes)))
    theta_mins = dist.Multinomial(probs=jnp.moveaxis(jnp.array([prob_score, prob_assist, prob_neither]), 0, 2))
    return numpyro.sample('obs', theta_mins, obs=y)