How to use pymc3 for survival analysis?

Questions
6

In that example, the survival function for the gumbel distribution is implemented like this:

def gumbel_sf(y, μ, σ):
    return 1. - tt.exp(-tt.exp(-(y - μ) / σ))

Am I missing something, or doesn’t this need to be log-transformed to get the log-likelihood?

1 Like