具有一般协方差结构线性模型的局部影响评价

该文讨论了具有一般协方差结构线性模型的局部影响分析问题. 通过对广义Cook统计量中M/c的适当选取, 文章给出了一种对扰动的参数变换具有不变性质的局部影响度量. 在具协方差扰动模式下, 该文给出了回归系数和方差系数估计、最佳线性预测的局部影响诊断统计量.该结果与数据删除法进行了比较, 并通过实例进行了分析和说明.     

Application of multilevel analysis in animal sciences

Multilevel modeling is considerably useful way to analyze hierarchical data sets. The main purpose of this paper is to apply multilevel analysis in animal science and also show that this modeling technique is appropriate to analyze this kind of data. Thus multilevel modeling technique is used to analyze the milk yield data which has hierarchical structures, sires nested within cows. As a result of the analysis done in this paper, it is obvious that multilevel modeling is needed to use for analyzing this data. This illustrates that it is a convenient way to use multilevel analysis for the data which obtained from animals when the data have hierarchies.     

离散型广义非线性纵向数据模型中偏离名义离差的检验及其功效模拟

离散型广义非线性模型包括Poisson,二项,负二项模型.本文讨论离散型广义非线性纵向数据模型中偏离名义离差的检验问题,得到了检验的score统计量,并利用MonteCarlo方法研究了检验统计量的性质.最后,利用杀虫剂数据说明了检验方法的应用.     

Bayesian networks with a logistic regression model for the conditional probabilities

Logistic regression techniques can be used to restrict the conditional probabilities of a Bayesian network for discrete variables. More specifically, each variable of the network can be modeled through a logistic regression model, in which the parents of the variable define the covariates. When all main effects and interactions between the parent variables are incorporated as covariates, the conditional probabilities are estimated without restrictions, as in a traditional Bayesian network. By incorporating interaction terms up to a specific order only, the number of parameters can be drastically reduced. Furthermore, ordered logistic regression can be used when the categories of a variable are ordered, resulting in even more parsimonious models. Parameters are estimated by a modified junction tree algorithm. The approach is illustrated with the Alarm network.     

Empirical model selection in generalized linear mixed effects models

This paper focuses on model selection in generalized linear mixed models using an information criterion approach. In these models in general, the response marginal distribution cannot be analytically derived. Thus, for parameter estimation, two approximations are revisited both leading to iterative model linearizations. We propose simple model selection criteria adapted from information criteria and based on the linearized model obtained at convergence of the algorithm. The quality of derived criteria are evaluated through simulations.     

Asymptotic properties in semiparametric partially linear regression models for functional data

We consider the semiparametric partially linear regression models with mean function XTβ + g(z), where X and z are functional data. The new estimators of β and g(z) are presented and some asymptotic results are given. The strong convergence rates of the proposed estimators are obtained. In our estimation, the observation number of each subject will be completely flexible. Some simulation study is conducted to investigate the finite sample performance of the proposed estimators.     

数据矩阵条件指数的局部影响分析

以广义影响函数和广义Cook统计量为工具,考察对数据矩阵同时进行微小扰动在条件指数上产生的局部影响,导出了寻找条件指数最大影响方向的计算公式,并用一个实例说明该方法的应用．     

Local asymptotic behavior of regression splines for marginal semiparametric models with longitudinal data

In this paper, we study the local asymptotic behavior of the regression spline estimator in the framework of marginal semiparametric model. Similarly to Zhu, Fung and He (2008), we give explicit expression for the asymptotic bias of regression spline estimator for nonparametric function f. Our results also show that the asymptotic bias of the regression spline estimator does not depend on the working covariance matrix, which distinguishes the regression splines from the smoothing splines and the seemingly u...     

非线性纵向数据模型中自相关性和随机效应的存在性检验

刻画纵向数据协方差结构有三种可能因素 ,即序列相关 (特别是一阶自相关 )、随机效应和常规的随机误差 (Diggleetal,2 0 0 2 ) .本文研究非线性纵向数据模型的自相关性和随机效应存在性的单个和联合检验 ,得到了检验的score统计量 ,并利用血浆药物渗透数据 (Davidian&Gilinan ,1 995)说明检验方法的应用 .     

Generalized partially linear mixed-effects models incorporating mismeasured covariates

In this article we consider a semiparametric generalized mixed-effects model, and propose combining local linear regression, and penalized quasilikelihood and local quasilikelihood techniques to estimate both population and individual parameters and nonparametric curves. The proposed estimators take into account the local correlation structure of the longitudinal data. We establish normality for the estimators of the parameter and asymptotic expansion for the estimators of the nonparametric part. For practical implementation, we propose an appropriate algorithm. We also consider the measurement error problem in covariates in our model, and suggest a strategy for adjusting the effects of measurement errors. We apply the proposed models and methods to study the relation between virologic and immunologic responses in AIDS clinical trials, in which virologic response is classified into binary variables. A dataset from an AIDS clinical study is analyzed.     

Shrinkage estimation for identification of linear components in additive models

In this short paper, we demonstrate that the popular penalized estimation method typically used for variable selection in parametric or semiparametric models can actually provide a way to identify linear components in additive models. Unlike most studies in the literature, we are NOT performing variable selection. Due to the difficulty in a priori deciding which predictors should enter the partially linear additive model as the linear components, such a method will prove useful in practice.     

Bayesian‐type count data models with varying coefficients: estimation and testing in the presence of overdispersion

In this paper we study varying‐coefficient models for count data. A Bayesian approach is taken to model the variability of the regression parameters. Based on a Kalman filter procedure the varying coefficients are estimated as the mode of the posterior distribution. All hyperparameters, including an overdispersion parameter in the negative binomial varying‐coefficient model (NBVC), are estimated as ML‐estimators using an EM‐type algorithm. A bootstrapping test of the fixed‐coefficient hypothesis against a varying‐coefficient alternative is proposed, which is evaluated running a simulation study. The study shows that the choice of a suitable count data model is of special importance in the framework of varying‐coefficient models. The methodology is illustrated analysing the determinants of the number of individual doctor visits. Copyright © 2001 John Wiley & Sons, Ltd.     

An imputation method for dealing with missing data in regression

We develop two methods for imputing missing values in regression situations. We examine the standard fixed-effects linear-regression model Y = X β + ?, where the regressors X are fixed and ? is the error term. This research focuses on the problem of missing X values. A particular component of market-share analysis has motivated this research where the price and other promotional instruments of each brand are allowed to have their own impact on the total sales volume in a consumer-products category. When a brand is not distributed in a particular week, only a few of the many measures occurring in that observation are missing. ‘What values should be imputed for the missing measures?’ is the central question this paper addresses. This context creates a unique problem in the missing-data literature, i.e. there is no true value for the missing measure. Using influence functions, from robust statistics we develop two loss functions, each of which is a function of the missing and existing X values. These loss functions turn out to be sums of ratios of low-order polynomials. The minimization of either loss function is an unconstrained non-linear-optimization problem. The solution to this non-linear optimization leads to imputed values that have minimal influence on the estimates of the parameters of the regression model. Estimates using the method for replacing missing values are compared with estimates obtained via some conventional methods.     

纵向数据下非参数带测量误差的部分线性变系数模型的估计

本文研究纵向数据下非参数部分带有测量误差的部分线性变系数模型的估计.利用B样条函数近似模型中的变系数函数,构造偏差修正的二次推断函数,得到模型中未知参数和变系数函数的估计.证明变系数函数估计量的相合性和参数估计量的渐近正态性.数值模拟和实例分析结果表明所提估计方法在有限样本下的有效性.     

指数族广义非线性随机系数模型的变离差检验

指数族广义非线性随机系数模型是Smith &; Heitjan［１０］和 Wei et al［１１］所研究模型的推广。该文分别在模型离差 (dispersion) 的权不变和变异时，讨论了指数族 广义非线性随机系数模型的变离差的检验问题，得到了score检验统计量。并利用欧洲野兔数据，分别对正态分布模型、Γ 分布模型和 逆高斯分布模型说明检验方法的有效性。     

A hidden Markov model for informative dropout in longitudinal response data with crisis states

We adopt a hidden state approach for the analysis of longitudinal data subject to dropout. Motivated by two applied studies, we assume that subjects can move between three states: stable, crisis, dropout. Dropout is observed but the other two states are not. During a possibly transient crisis state both the longitudinal response distribution and the probability of dropout can differ from those for the stable state. We adopt a linear mixed effects model with subject-specific trajectories during stable periods and additional random jumps during crises. We place the model in the context of Rubin’s taxonomy and develop the associated likelihood. The methods are illustrated using the two motivating examples.     

Robust Estimation in Partial Linear Mixed Model for Longitudinal Data

In this article, robust generalized estimating equation for the analysis of partial linear mixed model for longitudinal data is used. The authors approximate the nonparametric function by a regression spline. Under some regular conditions, the asymptotic properties of the estimators are obtained. To avoid the computation of high-dimensional integral, a robust Monte Carlo Newton-Raphson algorithm is used. Some simulations are carried out to study the performance of the proposed robust estimators. In addition, the authors also study the robustness and the efficiency of the proposed estimators by simulation. Finally, two real longitudinal data sets are analyzed.     

Residual Diagnostics for Covariate Effects in Spatial Point Process Models

For a spatial point process model in which the intensity depends on spatial covariates, we develop graphical diagnostics for validating the covariate effect term in the model, and for assessing whether another covariate should be added to the model. The diagnostics are point-process counterparts of the well-known partial residual plots (component-plus-residual plots) and added variable plots for generalized linear models. The new diagnostics can be derived as limits of these classical techniques under increasingly fine discretization, which leads to efficient numerical approximations. The diagnostics can also be recognized as integrals of the point process residuals, enabling us to prove asymptotic results. The diagnostics perform correctly in a simulation experiment. We demonstrate their utility in an application to geological exploration, in which a point pattern of gold deposits is modeled as a point process with intensity depending on the distance to the nearest geological fault. Online supplementary materials include technical proofs, computer code, and results of a simulation study.     

Quasi-Monte Carlo for Highly Structured Generalised Response Models

Highly structured generalised response models, such as generalised linear mixed models and generalised linear models for time series regression, have become an indispensable vehicle for data analysis and inference in many areas of application. However, their use in practice is hindered by high-dimensional intractable integrals. Quasi-Monte Carlo (QMC) is a dynamic research area in the general problem of high-dimensional numerical integration, although its potential for statistical applications is yet to be fully explored. We survey recent research in QMC, particularly lattice rules, and report on its application to highly structured generalised response models. New challenges for QMC are identified and new methodologies are developed. QMC methods are seen to provide significant improvements compared with ordinary Monte Carlo methods.      

我国宏观经济非参数联立模型的局部线性工具变量变窗宽估计

联立方程模型在经济政策制定、经济结构分析和预测方面起重要作用，目前关于非参数计量经济模型的研究主要停留在单方程模型上，而联立方程模型的研究在国际上刚刚起步，本将非参数回归模型的局部线性估计方法与传统联立方程模型估计方法相结合，首次提出了非参数计量经济联立模型的局部线性工具变量变窗宽估计并应用于我国宏观经济非参数联立模型，结果表明：我国宏观经济非参数联立模型优于线性联立模型且线性模型将造成不必要的人为设定误差；对于非参数联立模型，局部线性工具变量变窗宽估计优于局部线性估计。