One of the larger problems in modern statistics is that that telemetry (i.e., the data source) is chosen before analysis conducted.

[Yes it can be updated, yes a lot of good work goes into choices where events are recorded and how that data is recorded, and yes I really don’t care about the point you are trying to make right now… just read my article!]

The result of incomplete telemetry is that even if there are adequate model selection procedures (e.g., rjMCMC) there is still no guarantee that the data is available to accurately answer the parameter estimate question. This problem is has been discussed on this blog before; but that was in the context of GLM modelling.

This post explores the same problem within the context of survival analysis where for some bizarre accident of history dispersion is called frailty.

To be more precise, frailty is when two phone calls form the same population (i.e., have the same predictive statistics) have different true hang up time distributions. There could ether be unobservable sub-groups or a continuous spectrum.

The frailty model is where is a parameter of the survival function for the ith observation. There is an obvious parallel to Bayes theorem here. If is thought of as the prior distribution and the likelihood then is the evidence.

Note that I don’t really know why the custom is to express frailty models in terms of survival function. For example note that is the cdf of the distribution of the ith observation. Thus . Which implies that the process could have been started with the specification of the frail pdf of the ith observation .

An informative (simple-ish) example is the exponential-gamma frailty model. Consider a a population of phone calls with exponentially distributed hang-up times; i.e., . The complication is that the are distributed gamma(; i.e,

Thus to get the frailty model . this integral can be evaluated to be

That is it for now. Just a word of warning you still want a prior for and as they will not be constrained by the data. So don’t just MCMC sample from this thing yet!

### Like this:

Like Loading...

*Related*