New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
rstantools::loo_pit sometimes returns values larger than 1 #11
Comments
I'll look into this. Thanks for letting me know.
…On Thu, Aug 31, 2017 at 4:05 PM Aki Vehtari ***@***.***> wrote:
@jgabry <https://github.com/jgabry> has access to the data and then
should be able to reproduce
mod3 = stan_glmer(log_pm25 ~ log_sat_2014 + (log_sat_2014 |
cluster_region),
***@***.***)
psis3 <- psislw(-log_lik(mod3))
pit3 <- rstantools::loo_pit(object = posterior_predict(mod3), y =
***@***.***$pm25), lw = psis3$lw_smooth)
max(pit3)-1
[1] 8.881784e-16
pit values should be between 0 and 1
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#11>, or mute the thread
<https://github.com/notifications/unsubscribe-auth/AHb4Qy0dwh1GAvF665MA6zyi_s6Seqxmks5sdxH_gaJpZM4PJVt4>
.
|
Probably not safe to close, but @TeemuSailynoja is making a new pit function anyway |
Closed by #121 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
@jgabry has access to the data and then should be able to reproduce
mod3 = stan_glmer(log_pm25 ~ log_sat_2014 + (log_sat_2014 | cluster_region),
data=GMs@data)
psis3 <- psislw(-log_lik(mod3))
pit3 <- rstantools::loo_pit(object = posterior_predict(mod3), y = log(GMs@data$pm25), lw = psis3$lw_smooth)
pit values should be between 0 and 1
The text was updated successfully, but these errors were encountered: