Bayesian analysis of human movement curves

Summary.  We consider the Bayesian analysis of human movement data, where the subjects perform various reaching tasks. A set of markers is placed on each subject and a system of cameras records the three-dimensional Cartesian co-ordinates of the markers during the reaching movement. It is of interest to describe the mean and variability of the curves that are traced by the markers during one reaching movement, and to identify any differences due to covariates. We propose a methodology based on a hierarchical Bayesian model for the curves. An important part of the method is to obtain identifiable features of the movement so that different curves can be compared after temporal warping. We consider four landmarks and a set of equally spaced pseudolandmarks are located in between. We demonstrate that the algorithm works well in locating the landmarks, and shape analysis techniques are used to describe the posterior distribution of the mean curve. A feature of this type of data is that some parts of the movement data may be missing—the Bayesian methodology is easily adapted to cope with this situation.     

An application of cartographic area interpolation to German administrative data

In many situations the applied researcher wishes to combine different data sources without knowing the exact link and merging rule. This paper considers different cartographic interpolation methods for interpolating attributes from German employment office districts to German counties and vice versa. In particular, we apply dasymetric weighting as an alternative to simple area weighting, both of which are based on estimated intersection areas. We also present conditions under which the choice of interpolation method does not matter and confirm the theoretical results with a simulation study. Our application to German administrative data suggests robustness of estimation results of interpolated attributes with respect to the choice of interpolation method. We provide weighting matrices for regional data sources of the two largest German data producers.     

Easy-to-construct confidence bands for comparing two simple linear regression lines

This paper shows how to construct confidence bands for the difference between two simple linear regression lines. These confidence bands provide directly the information on the magnitude of the difference between the regression lines over an interval of interest and, as a by-product, can be used as a formal test of the difference between the two regression lines. Various different shapes of confidence bands are illustrated, and particular attention is paid towards confidence bands whose construction only involves critical points from standard distributions so that they are consequently easy to construct.     

Generalized model for IGC growth in aluminum alloys

A generalized brick wall model is developed to describe intergranular corrosion in Equi-axed AA7178-T6 and Wingskin AA7178-T6 aluminum alloys. The intergranular corrosion rate is highly related to grain size and shape. High strength aluminum alloys are often elongated and anisotropic, with the fastest nominal IGC growth rate in the longitudinal direction (L) or long transverse direction (T) and the slowest in the short transverse direction (S). We propose a three-way intersection model and use it to simulate the corrosion kinetics for each direction. With a proper combination of model parameters, the generalized IGC model provides a good fit to experimental data developed by the foil penetration technique.     

Wavelet goodness-of-fit test for dependent data

We introduce a new goodness-of-fit test which can be applied to hypothesis testing about the marginal distribution of dependent data. We derive a new test for the equivalent hypothesis in the space of wavelet coefficients. Such properties of the wavelet transform as orthogonality, localisation and sparsity make the hypothesis testing in wavelet domain easier than in the domain of distribution functions. We propose to test the null hypothesis separately at each wavelet decomposition level to overcome the problem of bi-dimensionality of wavelet indices and to be able to find the frequency where the empirical distribution function differs from the null in case the null hypothesis is rejected. We suggest a test statistic and state its asymptotic distribution under the null and under some of the alternative hypotheses.     

Bayesian inference for dynamic social network data

We consider a continuous-time model for the evolution of social networks. A social network is here conceived as a (di-) graph on a set of vertices, representing actors, and the changes of interest are creation and disappearance over time of (arcs) edges in the graph. Hence we model a collection of random edge indicators that are not, in general, independent. We explicitly model the interdependencies between edge indicators that arise from interaction between social entities. A Markov chain is defined in terms of an embedded chain with holding times and transition probabilities. Data are observed at fixed points in time and hence we are not able to observe the embedded chain directly. Introducing a prior distribution for the parameters we may implement an MCMC algorithm for exploring the posterior distribution of the parameters by simulating the evolution of the embedded process between observations.     

Retracted and replaced: Sodium cromoglicate: an ineffective drug or meta‐analysis misused?

The following article from Pharmaceutical Statistics, Sodium cromoglicate: an ineffective drug or meta‐analysis misused? by M. T. Stevens, A. M. Edwards, J. B. L. Howell published online on 15 March 2007 in Wiley InterScience ( www.interscience.wiley.com ), has been retracted by agreement between the author, the journal Editor in Chief, Steven Julious, and John Wiley & Sons, Ltd. The retraction has been agreed because the article is not yet ready for publication and an early version without revisions was published in error. Replacement article pending. Copyright © 2007 John Wiley & Sons, Ltd.     

H-likelihood: problems and solutions

In recent issues of this journal it has been asserted in two papers that the use of h-likelihood is wrong, in the sense of giving unsatisfactory estimates of some parameters for binary data (Kuk and Cheng, 1999; Waddington and Thompson, 2004) or theoretically unsound (Kuk and Cheng, 1999). We wish to refute both these assertions.     

Quantifying the Distribution of Inhalation Exposure in Human Populations: 2. Distributions of Time Spent by Adults, Adolescents, and Children at Home, at Work, and at School

Using distributions of time spent at various ventilation levels, ranges of inhalation exposure in the population can be established. Distributions of exposure time were determined using results of a study by the California Air Resources Board (CARB) which focused on time spent by humans participating in various activities and the locations where the activities occurred. The daily at-home activities from the CARB study were assigned to one of three ventilation levels, generating aggregate time periods. Distinct age and gender populations were identified, and distributions for aggregate time were established for these populations at each of the ventilation levels. In addition to aggregate time spent at home, distributions for various ages and genders were established for aggregate time spent at school and work. By combining distributions of aggregate time with corresponding ventilation rates, the distribution of inhalation rates can be established for at home, at work, and at school exposures.