Backfitting Random Varying-Coefficient Models with Time-Dependent Smoothing Covariates

Abstract.  In this paper, we propose a random varying-coefficient model for longitudinal data. This model is different from the standard varying-coefficient model in the sense that the time-varying coefficients are assumed to be subject-specific, and can be considered as realizations of stochastic processes. This modelling strategy allows us to employ powerful mixed-effects modelling techniques to efficiently incorporate the within-subject and between-subject variations in the estimators of time-varying coefficients. Thus, the subject-specific feature of longitudinal data is effectively considered in the proposed model. A backfitting algorithm is proposed to estimate the coefficient functions. Simulation studies show that the proposed estimation methods are more efficient in finite-sample performance compared with the standard local least squares method. An application to an AIDS clinical study is presented to illustrate the proposed methodologies.     

对基础研究投入问题的思考

文章首先从公共财政的角度对基础研究的投入进行分析,其次从布什和斯托克斯的科研模式角度对基础研究的投入进行分析。在此基础上对加大我国基础研究的投入提出了相关建议。     

Local influence in multilevel models

The authors study the local influence of observations in multilevel regression models. To this end, they perturb simultaneously the variances, responses and design matrix. To measure the local change caused by these perturbations, they use generalized Cook statistics for the fixed and random parameter estimates. Closed form local influence measures also allow them to assess the joint influence of various observations. They suggest a simple computation method and illustrate their results using two examples.     

Quantifying uncertainty in the biospheric carbon flux for England and Wales

Summary.  A crucial issue in the current global warming debate is the effect of vegetation and soils on carbon dioxide (CO₂) concentrations in the atmosphere. Vegetation can extract CO₂ through photosynthesis, but respiration, decay of soil organic matter and disturbance effects such as fire return it to the atmosphere. The balance of these processes is the net carbon flux. To estimate the biospheric carbon flux for England and Wales, we address the statistical problem of inference for the sum of multiple outputs from a complex deterministic computer code whose input parameters are uncertain. The code is a process model which simulates the carbon dynamics of vegetation and soils, including the amount of carbon that is stored as a result of photosynthesis and the amount that is returned to the atmosphere through respiration. The aggregation of outputs corresponding to multiple sites and types of vegetation in a region gives an estimate of the total carbon flux for that region over a period of time. Expert prior opinions are elicited for marginal uncertainty about the relevant input parameters and for correlations of inputs between sites. A Gaussian process model is used to build emulators of the multiple code outputs and Bayesian uncertainty analysis is then used to propagate uncertainty in the input parameters through to uncertainty on the aggregated output. Numerical results are presented for England and Wales in the year 2000. It is estimated that vegetation and soils in England and Wales constituted a net sink of 7.55 Mt C (1 Mt C = 10¹² g of carbon) in 2000, with standard deviation 0.56 Mt C resulting from the sources of uncertainty that are considered.     

Modelling heterogeneous space–time occurrences of earthquakes and its residual analysis

Summary. Earthquake intensities are modelled as a function of previous activity whose specific form is based on established empirical laws in seismology, but whose parameter values can vary from place to place. This model is used for characterizing regional features of seismic activities in and around Japan, and also for exploring regions where the actual seismicity rate systematically deviates from that of the modelled rate.     

Choice of the primary analysis in longitudinal clinical trials

Missing data, and the bias they can cause, are an almost ever‐present concern in clinical trials. The last observation carried forward (LOCF) approach has been frequently utilized to handle missing data in clinical trials, and is often specified in conjunction with analysis of variance (LOCF ANOVA) for the primary analysis. Considerable advances in statistical methodology, and in our ability to implement these methods, have been made in recent years. Likelihood‐based, mixed‐effects model approaches implemented under the missing at random (MAR) framework are now easy to implement, and are commonly used to analyse clinical trial data. Furthermore, such approaches are more robust to the biases from missing data, and provide better control of Type I and Type II errors than LOCF ANOVA. Empirical research and analytic proof have demonstrated that the behaviour of LOCF is uncertain, and in many situations it has not been conservative. Using LOCF as a composite measure of safety, tolerability and efficacy can lead to erroneous conclusions regarding the effectiveness of a drug. This approach also violates the fundamental basis of statistics as it involves testing an outcome that is not a physical parameter of the population, but rather a quantity that can be influenced by investigator behaviour, trial design, etc. Practice should shift away from using LOCF ANOVA as the primary analysis and focus on likelihood‐based, mixed‐effects model approaches developed under the MAR framework, with missing not at random methods used to assess robustness of the primary analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

BURR CRITICAL VALUE APPROXIMATIONS FOR TESTS OF AUTOCORRELATION AND HETEROSCEDASTICITY

The paper evaluates the accuracy of Burr approximations of critical values and p-values for test a of autocorrelation and heteroscedasticity in the linear regression model.     

用声学测量预报等温淬火球墨铸铁强度

等温淬火球墨铸铁的机械性能与其声学性能相关。本文将两种碳当量球铁在８４５，８９９和９４５℃温度进行奥氏体处理，在２６０，３１６和３７１℃温度进行等温淬火得到证实。声学性能测量包括共振频率，声速和内耗。由共振频率转换成动态模量表示出这种关系可以用来预报球墨铸铁不同热处理条件下的机械性能。     

关于线性互补问题的迭代算法

建立了一类求解线性互补问题的迭代算法。在一定条件下,研究了保证原问题的解存在唯一的充分条件,并且证明了新算法的收敛性。     

Some remarks on the physical models concerning two different approaches to inference in statistical process control

Summary For technological applications it can be useful to identify some simple physical mechanisms, which, on the basis of the available knowledge of the production process, may suggest the most appropriate approach to statistical control of the random quantities of interest. For this purpose the notion of rupture point is introduced firstly. A rupture point is characterized bym randomly arising out of control states, assumed to be mutually exclusive and stochastically independent. Shewhart's control charts seem to represent the natural statistical tool for controlling a rupture point; however it is shown that they are fully justified only when the hazard rates attached to the causes of failure are constant. Otherwise, typically in the presence of time increasing hazard rates, Shewhart's control charts should be completed by a preventive intervention rule (preventive maintenance). In the second place, the notion of dynamic instability point is introduced, which is specifically characterized by assuming that the random quantity of interest is ruled by a stochastic differential equation with constant coefficients. By discretization, developed according to a possibly new approach, it is shown that the former model reduces to an equation error model, which is among the simplest used in adaptive control, and thus particularly easy to deal with in regard to parameter estimation and the definition of the optimum control rule.