The size-density hypothesis in Great Britain: Analysis of a deviant case

Recent studies have found size of territorial units to vary inversely with population density, the only exception to this regularity being Great Britain, where size and density were found to be unrelated. The present research accounts for this anomaly by demonstrating Britain's historical conformity to the size-density relation. The size-density hypothesis is further supported by direct test, made possible by a recent reorganization of british counties. Results show that reorganization has restored the expected inverse relation between size and density and has restructured the size and density of counties in a manner precisely specified by underlying theory.     

Tests of fit for normal inverse Gaussian distributions

Goodness-of-fit tests for the family of symmetric normal inverse Gaussian distributions are constructed. The tests are based on a weighted integral incorporating the empirical characteristic function of suitably standardized data. An EM-type algorithm is employed for the estimation of the parameters involved in the test statistic. Monte Carlo results show that the new procedure is competitive with classical goodness-of-fit methods. An application with financial data is also included.     

Renyi's entropy for residual lifetime distribution

In the present paper, we introduce and study Renyi's information measure (entropy) for residual lifetime distributions. It is shown that the proposed measure uniquely determines the distribution. We present characterizations for some lifetime models. Further, we define two new classes of life distributions based on this measure. Various properties of these classes are also given.     

Making predictions from complex longitudinal data, with application to planning monitoring intervals in a national screening programme

When biological or physiological variables change over time, we are often interested in making predictions either of future measurements or of the time taken to reach some threshold value. On the basis of longitudinal data for multiple individuals, we develop Bayesian hierarchical models for making these predictions together with their associated uncertainty. Particular aspects addressed, which include some novel components, are handling curvature in individuals' trends over time, making predictions for both underlying and measured levels, making predictions from a single baseline measurement, making predictions from a series of measurements, allowing flexibility in the error and random-effects distributions, and including covariates. In the context of data on the expansion of abdominal aortic aneurysms over time, where reaching a certain threshold leads to referral for surgery, we discuss the practical application of these models to the planning of monitoring intervals in a national screening programme. Prediction of the time to reach a threshold was too imprecise to be practically useful, and we focus instead on limiting the probability of exceeding the threshold after given time intervals. Although more complex models can be shown to fit the data better, we find that relatively simple models seem to be adequate for planning monitoring intervals.     

Bayes Inference in Constant Partially Accelerated Life Tests for the Generalized Exponential Distribution with Progressive Censoring

In this paper, the problem of constant partially accelerated life tests when the lifetime follows the generalized exponential distribution is considered. Based on progressive type-II censoring scheme, the maximum likelihood and Bayes methods of estimation are used for estimating the distribution parameters and acceleration factor. A Monte Carlo simulation study is carried out to examine the performance of the obtained estimates.     

Bimodality based on the generalized skew-normal distribution

This paper focuses on the development of a new extension of the generalized skew-normal distribution introduced in Gómez et al. [Generalized skew-normal models: properties and inference. Statistics. 2006;40(6):495–505]. To produce the generalization a new parameter is introduced, the signal of which has the flexibility of yielding unimodal as well as bimodal distributions. We study its properties, derive a stochastic representation and state some expressions that facilitate moments derivation. Maximum likelihood is implemented via the EM algorithm which is based on the stochastic representation derived. We show that the Fisher information matrix is singular and discuss ways of getting round this problem. An illustration using real data reveals that the model can capture well special data features such as bimodality and asymmetry.     

A Note on Superiority Comparisons of Homogeneous Linear Estimators

This note extends some results on homogeneous linear estimators to the general, even nonlinear case.A Sufficient condition for the difference of mean square error matrices of minimum conditional mean square error estimator and minimum average risk linear estimator to be postive definite is derived.     

Frailty models for survival data

A frailty model is a random effects model for time variables, where the random effect (the frailty) has a multiplicative effect on the hazard. It can be used for univariate (independent) failure times, i.e. to describe the influence of unobserved covariates in a proportional hazards model. More interesting, however, is to consider multivariate (dependent) failure times generated as conditionally independent times given the frailty. This approach can be used both for survival times for individuals, like twins or family members, and for repeated events for the same individual. The standard assumption is to use a gamma distribution for the frailty, but this is a restriction that implies that the dependence is most important for late events. More generally, the distribution can be stable, inverse Gaussian, or follow a power variance function exponential family. Theoretically, large differences are seen between the choices. In practice, using the largest model makes it possible to allow for more general dependence structures, without making the formulas too complicated.This paper is a revised version of a review, which together with ten papers by the author made up a thesis for a Doctor of Science degree at the University of Copenhagen.