A counterexample to Beder's conjectures about Hadamard matrices

In this note we provide a counterexample which resolves conjectures about Hadamard matrices made in this journal. Beder [1998. Conjectures about Hadamard matrices. Journal of Statistical Planning and Inference 72, 7–14] conjectured that if H$\mathit{H}$ is a maximal m×n$m\times n$ row-Hadamard matrix then m is a multiple of 4; and that if n   is a power of 2 then every row-Hadamard matrix can be extended to a Hadamard matrix. Using binary integer programming we obtain a maximal 13×32$13\times 32$ row-Hadamard matrix, which disproves both conjectures. Additionally for n being a multiple of 4 up to 64, we tabulate values of m   for which we have found a maximal row-Hadamard matrix. Based on the tabulated results we conjecture that a m×n$m\times n$ row-Hadamard matrix with m?n-7$m?n-7$ can be extended to a Hadamard matrix.     

A robust QTL mapping procedure

In quantitative trait linkage studies using experimental crosses, the conventional normal location-shift model or other parameterizations may be unnecessarily restrictive. We generalize the mapping problem to a genuine nonparametric setup and provide a robust estimation procedure for the situation where the underlying phenotype distributions are completely unspecified. Classical Wilcoxon–Mann–Whitney statistics are employed for point and interval estimation of QTL positions and effects.     

Two-dimensional generalisations of dynamic programming for image analysis

Dynamic programming (DP) is a fast, elegant method for solving many one-dimensional optimisation problems but, unfortunately, most problems in image analysis, such as restoration and warping, are two-dimensional. We consider three generalisations of DP. The first is iterated dynamic programming (IDP), where DP is used to recursively solve each of a sequence of one-dimensional problems in turn, to find a local optimum. A second algorithm is an empirical, stochastic optimiser, which is implemented by adding progressively less noise to IDP. The final approach replaces DP by a more computationally intensive Forward-Backward Gibbs Sampler, and uses a simulated annealing cooling schedule. Results are compared with existing pixel-by-pixel methods of iterated conditional modes (ICM) and simulated annealing in two applications: to restore a synthetic aperture radar (SAR) image, and to warp a pulsed-field electrophoresis gel into alignment with a reference image. We find that IDP and its stochastic variant outperform the remaining algorithms.     

Asymptotic distributions for the intervals estimators of the extremal index and the cluster-size probabilities

Under appropriate long range dependence conditions, the point process of exceedances of a stationary sequence weakly converges to a homogeneous compound Poisson point process. This limiting point process can be characterized by the extremal index and the cluster-size probabilities. In this paper we address the problem of estimating these quantities and we consider the intervals estimators introduced in Ferro and Segers [2003. Inference for clusters of extreme values. J. Roy. Statist. Soc. Ser. B 545–556] and in Ferro [2004. Statistical methods for clusters of extreme values. Ph.D. Thesis, Lancaster University]. We establish asymptotic weak convergence to Gaussian random variables and we give their asymptotic variance.     

Nonparametric density deconvolution by weighted kernel estimators

Nonparametric density estimation in the presence of measurement error is considered. The usual kernel deconvolution estimator seeks to account for the contamination in the data by employing a modified kernel. In this paper a new approach based on a weighted kernel density estimator is proposed. Theoretical motivation is provided by the existence of a weight vector that perfectly counteracts the bias in density estimation without generating an excessive increase in variance. In practice a data driven method of weight selection is required. Our strategy is to minimize the discrepancy between a standard kernel estimate from the contaminated data on the one hand, and the convolution of the weighted deconvolution estimate with the measurement error density on the other hand. We consider a direct implementation of this approach, in which the weights are optimized subject to sum and non-negativity constraints, and a regularized version in which the objective function includes a ridge-type penalty. Numerical tests suggest that the weighted kernel estimation can lead to tangible improvements in performance over the usual kernel deconvolution estimator. Furthermore, weighted kernel estimates are free from the problem of negative estimation in the tails that can occur when using modified kernels. The weighted kernel approach generalizes to the case of multivariate deconvolution density estimation in a very straightforward manner.     

Automating and evaluating reversible jump MCMC proposal distributions

The reversible jump Markov chain Monte Carlo (MCMC) sampler (Green in Biometrika 82:711–732, 1995) has become an invaluable device for Bayesian practitioners. However, the primary difficulty with the sampler lies with the efficient construction of transitions between competing models of possibly differing dimensionality and interpretation. We propose the use of a marginal density estimator to construct between-model proposal distributions. This provides both a step towards black-box simulation for reversible jump samplers, and a tool to examine the utility of common between-model mapping strategies. We compare the performance of our approach to well established alternatives in both time series and mixture model examples.     

Cemetery locus as a mechanism of socio-cultural identity

In cultural studies of cemetery locus there is a very important aspect of understanding the cemetery as a tool for the formation of the socio-cultural identity of living people. I adhere to the point of view that the cemetery has always produced, and continues to produce, a variety of identities. While during the pre-Modern period the cemetery was a necessary element of individual self-understanding as a member of a certain community, in the Modern era the cemetery produces more particularistic identities. Modernity generates some universal and abstract schemes of identification, and by means of the repression of death from public consciousness, cemeteries lose their role as a focal point in social communication. I draw attention to the radical utopian ideas of the Russian philosopher Nikolai Fedorov, who not only considered the cemetery as a locus of memory, but also proclaimed the task of the transformation of cemeteries into a base for universal work on resurrection of dead ancestors and the restoration of brotherly relations in all mankind.     

Effects of unmeasured heterogeneity in the linear transformation model for censored data

We investigate the effect of unobserved heterogeneity in the context of the linear transformation model for censored survival data in the clinical trials setting. The unobserved heterogeneity is represented by a frailty term, with unknown distribution, in the linear transformation model. The bias of the estimate under the assumption of no unobserved heterogeneity when it truly is present is obtained. We also derive the asymptotic relative efficiency of the estimate of treatment effect under the incorrect assumption of no unobserved heterogeneity. Additionally we investigate the loss of power for clinical trials that are designed assuming the model without frailty when, in fact, the model with frailty is true. Numerical studies under a proportional odds model show that the loss of efficiency and the loss of power can be substantial when the heterogeneity, as embodied by a frailty, is ignored. An erratum to this article can be found at     

A computationally simple bivariate survival estimator for efficacy and safety

Both treatment efficacy and safety are typically the primary endpoints in Phase II, and even in some Phase III, clinical trials. Efficacy is frequently measured by time to response, death, or some other milestone event and thus is a continuous, possibly censored, outcome. Safety, however, is frequently measured on a discrete scale; in Eastern Cooperative Oncology Group clinical trial E2290, it was measured as the number of weekly rounds of chemotherapy that were tolerable to colorectal cancer patients. For the joint analysis of efficacy and safety, we propose a non-parametric, computationally simple estimator for the bivariate survival function when one time-to-event is continuous, one is discrete, and both are subject to right-censoring. The bivariate censoring times may depend on each other, but they are assumed to be independent of both event times. We derive a closed-form covariance estimator for the survivor function which allows for inference to be based on any of several possible statistics of interest. In addition, we derive its covariance with respect to calendar time of analysis, allowing for its use in sequential studies.     

Is the Study of Happiness a Worthy Scientific Pursuit?

This paper critiques the view that the study of happiness is not a worthy scientific pursuit. The happiness set point and hedonic treadmill theories denote the complexity of increasing happiness levels due to genetic limitations and adaptation, however, there is mounting evidence to suggest that with the use of appropriate measures and specific interventions aimed at fostering strengths and virtues, happiness can be increased. Furthermore, the benefits of investigating methods for increasing happiness include improvements in physical, psychological and social health and well-being. It is concluded that approaching human needs from a top down or holistic standpoint where individuals can use their strengths to overcome life’s challenges, is beneficial to health and well-being. Hence, the study of happiness is a worthy scientific pursuit.