rstanarm icon indicating copy to clipboard operation
rstanarm copied to clipboard

allow using student_t instead of gaussian for likelihood

Open jgabry opened this issue 10 years ago • 5 comments

For stan_glm with a suitable outcome (and maybe stan_lmer), we should provide an option to use the t distribution as the likelihood instead of a Gaussian.

Perhaps we just follow BDA3 17.2 and use the fact that

y_i ~ t_v(m, s^2)

is equivalent to

y_i | V_i ~ N(m, V_i) V_i ~ scaled-inv-chi-sq(v, s^2)

jgabry avatar Nov 09 '15 21:11 jgabry

We're going to also need a prior on the degrees of freedom. Aki suggested gamma with shape 2 and rate 0.1 https://github.com/stan-dev/stan/wiki/Prior-Choice-Recommendations

bgoodri avatar Nov 09 '15 21:11 bgoodri

@bgoodri Do you think the approach in my initial comment is preferable to implementing the student_t likelihood the same way we do the others? Might it be more efficient to not estimate all N of the V_i terms and instead just use y ~ student_t(...)? I'm not sure.

jgabry avatar Jan 18 '16 18:01 jgabry

Yeah, although we would have to create a family object for the t distribution.

On Mon, Jan 18, 2016 at 1:39 PM, Jonah Gabry [email protected] wrote:

@bgoodri https://github.com/bgoodri Do you think the approach in my initial comment is preferable to implementing the student_t likelihood the same way we do the others? Might it be more efficient to not estimate all N of the V_i terms and instead just use y ~ student_t(...)? I'm not sure.

— Reply to this email directly or view it on GitHub https://github.com/stan-dev/rstanarm/issues/35#issuecomment-172619206.

bgoodri avatar Jan 18 '16 18:01 bgoodri

Hi All - just wondering if adding a Student's t distribution for likelihood is still being considered? Or was this effort conclusively discontinued?

maciematt avatar Feb 12 '20 18:02 maciematt

We certainly haven't ruled out including it; we just haven't included it yet.

On Wed, Feb 12, 2020 at 1:09 PM Matt Maciejewski [email protected] wrote:

Hi All - just wondering if adding a Student's t distribution for likelihood is still being considered? Or was this effort conclusively discontinued?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/stan-dev/rstanarm/issues/35?email_source=notifications&email_token=AAZ2XKR2YGCFFDDXLO47O5DRCQ3PNA5CNFSM4BT3X7YKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOELRZLEQ#issuecomment-585340306, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAZ2XKQIUY6NAPDWKT4YSTLRCQ3PNANCNFSM4BT3X7YA .

bgoodri avatar Feb 12 '20 18:02 bgoodri