Parametric Likelihood Inference for Record Breaking Problems
本文提出一种参数似然推断方法,通过蒙特卡洛积分近似处理记录突破数据中的高维积分,适用于依赖性和均值漂移的观测,并以奥运跳高记录为例展示模型选择与预测。
In this paper we consider the analysis of record breaking data sets, where only observations that exceed, or only those that fall below, the current extreme value are recorded. Example application areas include industrial stress testing, meteorological analysis, sporting and athletic events, and oil and mining surveys. A closely related area is that of threshold modelling, where the observations are those that cross a certain threshold value. The inherent missing data structure present in these problems leads to likelihood functions that contain possibly high-dimensional integrals, rendering traditional maximum likelihood methods difficult or not feasible. Fortunately, we may obtain arbitrarily accurate approximations to the likelihood function by iteratively applying Monte Carlo integration methods (Geyer & Thompson, 1992). Subiteration using the Gibbs sampler may help to evaluate any multivariate integrals encountered during this process. This approach can handle far more sophisticated parametric models than have been used previously in record breaking and threshold data contexts. In particular, the methodology allows for observations that are dependent and subject to mean shifts over time. We present a numerical example involving records in Olympic high jump competition, where besides estimation we also address related issues in model selection and prediction.