Optimal estimation in random coefficient regression models

In linear regression models with random coefficients, the score function usually involves unknown nuisance parameters in the form of weights. Conditioning with respect to the sufficient statistics for the nuisance parameter, when the parameter of interest is held fixed, eliminates the nuisance parameters and is expected to give reasonably good estimating functions. The present paper adopts this approach to the problem of estimation of average slope in random coefficient regression models. Four sampling situations are discussed. Some asymptotic results are also obtained for a model where neither the regressors nor the random regression coefficients replicate. Simulation studies for normal as well as non-normal models show that the performance of the suggested estimating functions is quite satisfactory.     

Semiparametric Bayesian inference for mean-covariance regression models

In this paper, we propose a Bayesian semiparametric mean-covariance regression model with known covariance structures. A mixture model is used to describe the potential non-normal distribution of the regression errors. Moreover, an empirical likelihood adjusted mixture of Dirichlet process model is constructed to produce distributions with given mean and variance constraints. We illustrate through simulation studies that the proposed method provides better estimations in some non-normal cases. We also demonstrate the implementation of our method by analyzing the data set from a sleep deprivation study.     

带随机效应非线性模型的影响分析

Abstract. In this paper,a unified diagnostic method for the nonlinear models with random ef-fects based upon the joint likelihood given by Robinson in 1991 is presented. It is shown that thecase deletion model is equivalent to the mean shift outlier model. From this point of view ,sever-al diagnostic measures, such as Cook distance, score statistics are derived. The local influencemeasure of Cook is also presented. A numerical example illustrates that the method is avail-able     

Additive hazards regression with random effects for clustered failure times

Additive hazards model with random effects is proposed for modelling the correlated failure time data when focus is on comparing the failure times within clusters and on estimating the correlation between failure times from the same cluster, as well as the marginal regression parameters. Our model features that, when marginalized over the random effect variable, it still enjoys the structure of the additive hazards model. We develop the estimating equations for inferring the regression parameters. The proposed estimators are shown to be consistent and asymptotically normal under appropriate regularity conditions. Furthermore, the estimator of the baseline hazards function is proposed and its asymptotic properties are also established. We propose a class of diagnostic methods to assess the overall fitting adequacy of the additive hazards model with random effects. We conduct simulation studies to evaluate the finite sample behaviors of the proposed estimators in various scenarios. Analysis of the Diabetic Retinopathy Study is provided as an illustration for the proposed method.     

The estimation of multivariate random coefficient autoregressive models

Using a two stage regression procedure estimates of the unknown parameters of a class of multivariate random coefficient autoregressive models are obtained. The estimates are shown, under fairly general conditions, to be strongly consistent and to have a distribution which converges to that of a normally distributed random vector.     

Hierarchical linear regression models for conditional quantiles

The quantile regression has several useful features and therefore is gradually developing into a comprehensive approach to the statistical analysis of linear and nonlinear response models, but it cannot deal effectively with the data with a hierarchical structure. In practice, the existence of such data hierarchies is neither accidental nor ignorable, it is a common phenomenon. To ignore this hierarchical data structure risks overlooking the importance of group effects, and may also render many of the traditional statistical analysis techniques used for studying data relationships invalid. On the other hand, the hierarchical models take a hierarchical data structure into account and have also many applications in statistics, ranging from overdispersion to constructing min-max estimators. However, the hierarchical models are virtually the mean regression, therefore, they cannot be used to characterize the entire conditional distribution of a dependent variable given high-dimensional covariates. Furthermore, the estimated coefficient vector (marginal effects) is sensitive to an outlier observation on the dependent variable. In this article, a new approach, which is based on the Gauss-Seidel iteration and taking a full advantage of the quantile regression and hierarchical models, is developed. On the theoretical front, we also consider the asymptotic properties of the new method, obtaining the simple conditions for an n1/2-convergence and an asymptotic normality. We also illustrate the use of the technique with the real educational data which is hierarchical and how the results can be explained.     

Empirical likelihood estimators for the error distribution in nonparametric regression models

The aim of this paper is to show that existing estimators for the error distribution in non-parametric regression models can be improved when additional information about the distribution is included by the empirical likelihood method. The weak convergence of the resulting new estimator to a Gaussian process is shown and the performance is investigated by comparison of asymptotic mean squared errors and by means of a simulation study.      

基于纵向数据的Poisson-Gamma模型系数的随机性检验

应用随机系数化方法对纵向数据的Poisson-Gamma回归模型进行了研究,采用Laplace展开方法得到了关于响应变量的近似似然函数并得到模型系数的随机性检验的Score统计量.通过Monte Carlo模拟分析了检验的渐近功效.最后把得到的检验统计量应用到具体的数值实例分析中说明其有效性.     

随机效应模型中方差分量的经验Bayes检验问题

给出了双向分类随机效应模型中方差分量的Bayes检验的判决函数,利用核估计的方法,构造了相应的经验Bayes(EB)检验的判决函数．在适当的条件下证明了EB判决函数是渐近最优的且有收敛速度．给出了模型的特例和推广．最后,举出一个满足定理条件的例子．     

Shrinkage estimation and variable selection in multiple regression models with random coefficient autoregressive errors

In this paper, we consider improved estimation strategies for the parameter vector in multiple regression models with first-order random coefficient autoregressive errors (RCAR(1)). We propose a shrinkage estimation strategy and implement variable selection methods such as lasso and adaptive lasso strategies. The simulation results reveal that the shrinkage estimators perform better than both lasso and adaptive lasso when and only when there are many nuisance variables in the model.     

Empilrical likelihood for non-parametric regression models with missing responses: Multiple design case

Empirical likelihood(EL) ratio statistic on θ = g(x) is constructed based on the inverse probability weighted imputation approach in a nonparametric regression model Y = g(x) + ε(x ∈ [0,1]p) with fixed designs and missing responses,which asymptotically has χ12 distribution.This result is used to obtain a EL based confidence interval on θ.     

Two-component Poisson mixture regression modelling of count data with bivariate random effects

Two-component Poisson mixture regression is typically used to model heterogeneous count outcomes that arise from two underlying sub-populations. Furthermore, a random component can be incorporated into the linear predictor to account for the clustering data structure. However, when including random effects in both components of the mixture model, the two random effects are often assumed to be independent for simplicity. A two-component Poisson mixture regression model with bivariate random effects is proposed to deal with the correlated situation. A restricted maximum quasi-likelihood estimation procedure is provided to obtain the parameter estimates of the model. A simulation study shows both fixed effects and variance component estimates perform well under different conditions. An application to childhood gastroenteritis data demonstrates the usefulness of the proposed methodology, and suggests that neglecting the inherent correlation between random effects may lead to incorrect inferences concerning the count outcomes.     

On familial longitudinal Poisson mixed models with gamma random effects

Poisson mixed models are used to analyze a wide variety of cluster count data. These models are commonly developed based on the assumption that the random effects have either the log-normal or the gamma distribution. Obtaining consistent as well as efficient estimates for the parameters involved in such Poisson mixed models has, however, proven to be difficult. Further problem gets mounted when the data are collected repeatedly from the individuals of the same cluster or family. In this paper, we introduce a generalized quasilikelihood approach to analyze the repeated familial data based on the familial structure caused by gamma random effects. This approach provides estimates of the regression parameters and the variance component of the random effects after taking the longitudinal correlations of the data into account. The estimators are consistent as well as highly efficient.     

Testing for random effects in linear mixed models for longitudinal data under moment conditions

In this paper, we consider whether the random effect exists in linear mixed models （LMMs） when only moment conditions are assumed. Based on the estimators of parameters and their asymptotic properties, a Wald-type test is constructed. It is consistent against global alternatives and is sensitive to the local alternatives converging to the null hypothesis at parametric rates, a fastest possibly rate for goodness-of-fit testing. Moreover, a simulation study shows the performance of the test is good. The procedure also applies to a real data.     

SOME ASYMPTOTIC INFERENCE ON NONLINEAR MODELS WITH RANDOM EFFECTS （A GEOMETRIC APPROACH）

ＳＯＭＥＡＳＹＭＰＴＯＴＩＣＩＮＦＥＲＥＮＣＥＯＮＮＯＮＬＩＮＥＡＲＭＯＤＥＬＳＷＩＴＨＲＡＮＤＯＭＥＦＦＥＣＴＳ（ＡＧＥＯＭＥＴＲＩＣＡＰＰＲＯＡＣＨ）ＺＨＯＮＧＸＵＰＩＮＧＡＮＤＷＥＩＢＯＣＨＥＮＧＡｂｓｔｒａｃｔ．Ａｇｅｏｍｅｔｒｉｃｆｒａｍｅ...     

Symmetric sum and symmetric product of two independent random variables

The symmetry of the sum and product of two independent random variablesX andY, whenX orY is not symmetric, is studied.     

Complete moment convergence for weighted sums of widely orthant-dependent random variables and its application in nonparametric regression models

We establish some results on the complete moment convergence for weighted sums of widely orthant-dependent (WOD) random variables, which improve and extend the corresponding results of Y. F. Wu, M. G. Zhai, and J. Y. Peng [J. Math. Inequal., 2019, 13(1): 251–260]. As an application of the main results, we investigate the complete consistency for the estimator in a nonparametric regression model based on WOD errors and provide some simulations to verify our theoretical results.     

Methods for generating random variates with Polya characteristic functions

Polya has shown that real even continuous functions that are convex on (0,∞), for 1 t = 0, and decreasing to 0 as t → ∞ are characteristic functions. Dugué and Girault (1955) have shown that the corresponding random variables are distributed as $\text{Y}\text{Z}$ where Y is a random variable with density (2π)^?1$\text{(}\text{sin}\text{(}\text{x}\text{2}\text{)/(}\text{x}\text{2}\text{))}{}^2$, and Z is independent of Y and has distribution function 1 ? φ + tφ′, t > 0. This property allows us to develop fast algorithms for this class of distributions. This is illustrated for the symmetric stable distribution, Linnik's distribution and a few other distributions. We pay special attention to the generation of Y.