search

lnccbrown / HSSM

Development of HSSM package
Other
71 stars 10 forks source link

link_settings="log_logit" works only with LANs #354

Closed igrahek closed 1 month ago

igrahek commented 4 months ago

HSSM version: v0.2.0

Code from the current tutorial on using logits runs with LANs:

cav_data = hssm.load_data("cavanagh_theta")

model_safe = hssm.HSSM(
    data=cav_data,
    hierarchical=True,
    link_settings="log_logit",
    loglik_kind="approx_differentiable",
)

model_safe.sample()

However, it throws: "pymc.logprob.utils.ParameterValueError: a >= 0" when using the analytical likelihood:

cav_data = hssm.load_data("cavanagh_theta")

model_safe = hssm.HSSM(
    data=cav_data,
    hierarchical=True,
    link_settings="log_logit",
    loglik_kind="analytical",
)

model_safe.sample()