Bayes shrinkage estimator for consistency assessment of treatment effects in multi-regional clinical trials

The primary objective of a multi-regional clinical trial is to investigate the overall efficacy of the drug across regions and evaluate the possibility of applying the overall trial result to some specific region. A challenge arises when there is not enough regional sample size. We focus on the problem of evaluating applicability of a drug to a specific region of interest under the criterion of preserving a certain proportion of the overall treatment effect in the region. We propose a variant of James-Stein shrinkage estimator in the empirical Bayes context for the region-specific treatment effect. The estimator has the features of accommodating the between-region variation and finiteness correction of bias. We also propose a truncated version of the proposed shrinkage estimator to further protect risk in the presence of extreme value of regional treatment effect. Based on the proposed estimator, we provide the consistency assessment criterion and sample size calculation for the region of interest. Simulations are conducted to demonstrate the performance of the proposed estimators in comparison with some existing methods. A hypothetical example is presented to illustrate the application of the proposed method.     

Infill asymptotics for adaptive kernel estimators of spatial intensity

We apply the Abramson principle to define adaptive kernel estimators for the intensity function of a spatial point process. We derive asymptotic expansions for the bias and variance under the regime that n independent copies of a simple point process in Euclidean space are superposed. The method is illustrated by means of a simple example and applied to tornado data.     

Smooth tail-index estimation

The two parametric distribution functions appearing in the extreme-value theory – the generalized extreme-value distribution and the generalized Pareto distribution – have log-concave densities if the extreme-value index γ∈[?1, 0]. Replacing the order statistics in tail-index estimators by their corresponding quantiles from the distribution function that is based on the estimated log-concave density ? f _n leads to novel smooth quantile and tail-index estimators. These new estimators aim at estimating the tail index especially in small samples. Acting as a smoother of the empirical distribution function, the log-concave distribution function estimator reduces estimation variability to a much greater extent than it introduces bias. As a consequence, Monte Carlo simulations demonstrate that the smoothed version of the estimators are well superior to their non-smoothed counterparts, in terms of mean-squared error.     

Rank invariant tests with left truncated and interval censored data

Peto and Peto (1972) have studied rank invariant tests to compare two survival curves for right censored data. We apply their tests, including the logrank test and the generalized Wilcoxon test, to left truncated and interval censored data. The significance levels of the tests are approximated by Monte Carlo permutation tests. Simulation studies are conducted to show their size and power under different distributional differences. In particular, the logrank test works well under the Cox proportional hazards alternatives, as for the usual right censored data. The methods are illustrated by the analysis of the Massachusetts Health Care Panel Study dataset.     

Estimating the variance of the sample median

The small-sample bias and root mean squared error of several distribution-free estimators of the variance of the sample median are examined. A new estimator is proposed that is easy to compute and tends to have the smallest bias and root mean squared error.     

On Pitman's measure of closeness of k-records

Pitman closeness of both the upper and lower k-record statistics to the population quantiles of a location–scale family of distributions is studied. For the population median, the Pitman-closest k-record is also determined. In the case of symmetric distributions, the Pitman closeness probabilities of k-record statistics are shown to be distribution-free, and explicit expressions are also derived for these probabilities. Exact expressions are derived for the required probabilities for uniform and exponential distributions. Numerical results are given for these families and also the Pitman-closest k-record is determined.     

Daily Inhalation Rate and Time‐Activity/Location Pattern in Japanese Preschool Children

Lack of data on daily inhalation rate and activity of children has been an issue in health risk assessment of air pollutants. This study aimed to obtain the daily inhalation rate and intensity and frequency of physical activity in relation to the environment in Japanese preschool children. Children aged four–six years (n= 138) in the suburbs of Tokyo participated in this study, which involved three days' continuous monitoring of physical activity using a tri‐axial accelerometer and parent's completion of a time/location diary during daily life. The estimated three‐day mean daily inhalation rate (body temperature, pressure, saturated with water vapor) was 9.9 ± 1.6 m³/day (0.52 ± 0.09 m³/kg/day). The current daily inhalation rate value of 0.580 m³/kg/day proposed for use in health risk assessment in Japan is confirmed to be valid to calculate central value of inhaled dose of air pollutants in five‐ to six‐year‐old children. However, the 95th percentile daily inhalation rate of 0.83 m³/kg/day based on measurement for five‐year‐old children is recommended to be used to provide an upper bound estimate of exposure that ensure the protection of all five‐ to six‐year‐old children from the health risk of air pollutants. Children spent the majority of their time in sedentary and light level of physical activity (LPA) when indoors, while 85% of their time when outdoors was spent in LPA and moderate‐to‐vigorous physical activity. The results suggest the need to consider variability of minute respiratory ventilation rate according to the environment for more refined short‐term health risk assessment.     

On tests of symmetry,marginal homogeneity and quasi-symmetry in two-way contingency tables based on minimum φ-divergence estimator with constraints

The restricted minimum φ-divergence estimator, [Pardo, J.A., Pardo, L. and Zografos, K., 2002, Minimum φ-divergence estimators with constraints in multinomial populations. Journal of Statistical Planning and Inference, 104, 221–237], is employed to obtain estimates of the cell frequencies of an I×I contingency table under hypotheses of symmetry, marginal homogeneity or quasi-symmetry. The associated φ-divergence statistics are distributed asymptotically as chi-squared distributions under the null hypothesis. The new estimators and test statistics contain, as particular cases, the classical estimators and test statistics previously presented in the literature for the cited problems. A simulation study is presented, for the symmetry problem, to choose the best function φ₂ for estimation and the best function φ₁ for testing.     

Optimal determination of the parameters of some biased estimators using genetic algorithm

In this paper, some new algorithms for estimating the biasing parameters of the ridge, Liu and two-parameter estimators are introduced with the help of genetic algorithm (GA). The proposed algorithms are based on minimizing some statistical measures such as mean square error (MSE), mean absolute error (MAE) and mean absolute prediction error (MAPE). At the same time, the new algorithms allow one to keep the condition number and variance inflation factors to be less than or equal to ten by means of the GA. A numerical example is presented to show the utility of the new algorithms. In addition, an extensive Monte Carlo experiment is conducted. The numerical findings prove that the proposed algorithms enable to eliminate the problem of multicollinearity and minimize the MSE, MAE and MAPE.     

Editorial collaborators

A simulation study of the binomial-logit model with correlated random effects is carried out based on the generalized linear mixed model (GLMM) methodology. Simulated data with various numbers of regression parameters and different values of the variance component are considered. The performance of approximate maximum likelihood (ML) and residual maximum likelihood (REML) estimators is evaluated. For a range of true parameter values, we report the average biases of estimators, the standard error of the average bias and the standard error of estimates over the simulations. In general, in terms of bias, the two methods do not show significant differences in estimating regression parameters. The REML estimation method is slightly better in reducing the bias of variance component estimates.