Maximum penalized likelihood estimation of mixed proportional hazard models

Unobservable individual effects in models of duration will cause estimation bias that include the structural parameters as well as the duration dependence. The maximum penalized likelihood estimator is examined as an estimator for the survivor model with heterogeneity. Proofs of the existence and uniqueness of the maximum penalized likelihood estimator in duration model with general forms of unobserved heterogeneity are provided. Some small sample evidence on the behavior of the maximum penalized likelihood estimator is given. The maximum penalized likelihood estimator is shown to be computationally feasible and to provide reasonable estimates in most cases.     

Combined estimators for the correlation coefficient

Three combined estimators for the bivariate normal correlation parameter are considered. The data consist of k independent sample correlation coefficients and it is assumed that the underlying correlation parameters are all equal to ρ. Based upon the joint density function of the sample correlations a combined estimator of ρ is obtained as an approximation to the maximum likelihood solution. Two linearly combined estimators are also considered. One of them is based on Fisher's z-transformation of the sample correlations and the other on an unbiased estimator of ρ. The comparison of these three estimators indicates that the combined (approximate) MLE has a slightly smaller estimated mean squared error relative to the other two combined methods of estimation, but it does so at the expense of a relatively larger bias.     

Estimating spatial means with an application to remotely sensed data

The paper examines alternative estimators for the mean of a spatial process where observations are not independent. Properties of the sample mean and its standard error are contrasted with those of maximum likelihood estimators derived for three spatial models. The information loss caused by spatial dependency in the data is examined. The distribution theory for the estimators is reviewed and the paper concludes with an empirical example illustrating the properties of the estimators and the practical benefits of the maximum likelihood procedure.     

On the characterization of maximum likelihood estimators for location-scale families

Given maximum likelihood equations for location and scale parameters, one determines conditions under which there exists a uniquely defined parametric statistical model, whose location and scale maximum likelihood estimators are the given ones. The constructive approach is exemplified at several kinds of mean estimators including the mean, mean square, mean mean and stretched power mean. The possible extension of the method to more general situations is discussed and illustrated at the sample median maximum likelihood estimator.     

Analysis of spatial frailty models by a weighted estimating equation

Spatially correlated survival data are frequently observed in ecological and epidemiological studies. An assumption in the clustered survival models is inter-cluster independence, which may not be adequate to model the dependence in spatial settings. For survival data, the likelihood function based on a spatial frailty may be complicated. In this paper, we develop a weighted estimating equation for spatially right-censored data. Some large sample properties for the estimate are developed. We also conduct simulations to compare estimation performance with other methods. A data set from a study of forest decline in Wisconsin is used to illustrate the proposed method.     

Classification Rules for Multivariate Repeated Measures Data with Equicorrelated Correlation Structure on Both Time and Spatial Repeated Measurements

We study the problem of classification for multivariate repeated measures data with structured correlations on both time and spatial repeated measurements. This is a very important problem in many biomedical as well as in engineering field. Classification rules as well as the algorithm to compute the maximum likelihood estimates of the required parameters are given.     

Likelihood-based estimation for Gaussian MRFs

Markov random fields (MRFs) express spatial dependence through conditional distributions, although their stochastic behavior is defined by their joint distribution. These joint distributions are typically difficult to obtain in closed form, the problem being a normalizing constant that is a function of unknown parameters. The Gaussian MRF (or conditional autoregressive model) is one case where the normalizing constant is available in closed form; however, when sample sizes are moderate to large (thousands to tens of thousands), and beyond, its computation can be problematic. Because the conditional autoregressive (CAR) model is often used for spatial-data modeling, we develop likelihood-inference methodology for this model in situations where the sample size is too large for its normalizing constant to be computed directly. In particular, we use simulation methodology to obtain maximum likelihood estimators of mean, variance, and spatial-depencence parameters (including their asymptotic variances and covariances) of CAR models.     

Maximum likelihood estimation foe two parameter exponentials under type I censoring

The paper considers maximum likelihood estimation of the location parameters and the failure rates of two parameter exponentials under type I censoring. Both the one and two sample cases are considered. It turns out that by debiasing the maximum likelihood estimators of the location parameters, one can achieve asymptotically 50% mean squared error reduction.     

Analysis of clustered spatially correlated binary data using autologistic model and Bayesian method with an application to dental caries of 3–5-year-old children

The autologistic model, first introduced by Besag, is a popular tool for analyzing binary data in spatial lattices. However, no investigation was found to consider modeling of binary data clustered in uncorrelated lattices. Owing to spatial dependency of responses, the exact likelihood estimation of parameters is not possible. For circumventing this difficulty, many studies have been designed to approximate the likelihood and the related partition function of the model. So, the traditional and Bayesian estimation methods based on the likelihood function are often time-consuming and require heavy computations and recursive techniques. Some investigators have introduced and implemented data augmentation and latent variable model to reduce computational complications in parameter estimation. In this work, the spatially correlated binary data distributed in uncorrelated lattices were modeled using autologistic regression, a Bayesian inference was developed with contribution of data augmentation and the proposed models were applied to caries experiences of deciduous dents.     

Maximum likelihood estimation and prediction mean square error in the spatial linear model

This paper is concerned with prediction in the spatial linear model using the maximum likelihood estimation of parameters in this model. In particular, we give some properties of predictors obtained on substituting the maximum likelihood estimators of model parameters into the form of the best-in the sense of minimum mean square prediction error-predictor. Such predictors are not optimal but we show them to be asymptotically equivalent to the optimum. We discuss practical aspects of this work and conclude by considering the connection with other areas.     

LEVERAGE ADJUSTMENTS FOR DISPERSION MODELLING IN GENERALIZED NONLINEAR MODELS

For normal linear models, it is generally accepted that residual maximum likelihood estimation is appropriate when covariance components require estimation. This paper considers generalized linear models in which both the mean and the dispersion are allowed to depend on unknown parameters and on covariates. For these models there is no closed form equivalent to residual maximum likelihood except in very special cases. Using a modified profile likelihood for the dispersion parameters, an adjusted score vector and adjusted information matrix are found under an asymptotic development that holds as the leverages in the mean model become small. Subsequently, the expectation of the fitted deviances is obtained directly to show that the adjusted score vector is unbiased at least to O(1/n) . Exact results are obtained in the single‐sample case. The results reduce to residual maximum likelihood estimation in the normal linear case.     

Reml estimation for repeated measures analysis

Models for repeated measures or growth curves consist of a mean response plus error and the errors are usually correlated. Both maximum likelihood and residual maximum likelihood (REML) estimators of a regression model with dependent errors are derived for cases in which the variance matrix of the error model admits a convenient Cholesky factorisation. This factorisation may be linked to methods for producing recursive estimates of the regression parameters and recursive residuals to provide a convenient computational method. The method is used to develop a general approach to repeated measures analysis.     

Classification of Higher-order Data with Separable Covariance and Structured Multiplicative or Additive Mean Models

Although devised in 1936 by Fisher, discriminant analysis is still rapidly evolving, as the complexity of contemporary data sets grows exponentially. Our classification rules explore these complexities by modeling various correlations in higher-order data. Moreover, our classification rules are suitable to data sets where the number of response variables is comparable or larger than the number of observations. We assume that the higher-order observations have a separable variance-covariance matrix and two different Kronecker product structures on the mean vector. In this article, we develop quadratic classification rules among g different populations where each individual has κth order (κ ≥2) measurements. We also provide the computational algorithms to compute the maximum likelihood estimates for the model parameters and eventually the sample classification rules.     

Using a correlated probit model approximation to estimate the variance for binary matched pairs

A correlated probit model approximation for conditional probabilities (Mendell and Elston 1974) is used to estimate the variance for binary matched pairs data by maximum likelihood. Using asymptotic data, the bias of the estimates is shown to be small for a wide range of intra-class correlations and incidences. This approximation is also compared with other recently published, or implemented, improved approximations. For the small sample examples presented, it shows a substantial advantage over other approximations. The method is extended to allow covariates for each observation, and fitting by iteratively reweighted least squares.     

A Note on Distinguishing Random Trees Populations

ABSTRACT

Several new results are presented for a class of univariate distributions for which the maximum likelihood estimate of the population mean is the sample mean. It is shown that the convolution of any two such distributions also belongs to this class of functions. It is also shown that the marginal distribution for the sample mean captures all of the Fisher information for the population mean contained in the full distribution. Parameters orthogonal to the mean are found for special cases of these distributions. If the distribution is conditioned on the sample mean, the conditional distribution depends on the parameters only through parameters orthogonal to the mean.     

双重滞后随机前沿模型技术效率的估计

首次在随机前沿模型中同时引入因变量间(或双边误差间)和技术效率间的空间相关性并构造了双重滞后随机前沿模型,使用极大似然估计方法和JLMS方法得出参数和技术效率的估计。蒙特卡罗模拟表明:忽略技术效率的空间相关性,参数估计和技术效率的估计均表现欠佳。本研究能以较高的精度估计参数和技术效率。随着样本容量的增加,估计效果更优。     

Estimating the correlation coefficient in the presence of correlated observations from a bivariate normal population

Estimation of the correlation coefficient between two variates (p) in the presence of correlated observations from a bivar iate normal population is considered The estimated maximum likelihood estimator (EMLE), an estimate based on the maximum likelihood estimator (MLE), is proposed and studied for the estimation of p For the large sample case , approximate expressions foi the variance and the bias of the Pearson estimate of the correlation coefficient are derived. These expressions suggests that the Pearson’s estimator possesses high mean square error (MSE) in estimating ρ in comparison to the MLE The MSE is particularly high when the observations within clusters aie highly correlated. The Pearson’s estimate, the MLE, and the EMLE aie evaluated in a simulation study This study shows that the proposed EMLE pefoims bettei than the Pearson’s correlation coefficient except when the number of clusters is small.     

A bootstrap approximation to the joint distribution of sum and maximum of a stationary sequence

This paper establishes the asymptotic validity for the moving block bootstrap as an approximation to the joint distribution of the sum and the maximum of a stationary sequence. An application is made to statistical inference for a positive time series where an extreme value statistic and sample mean provide the maximum likelihood estimates for the model parameters. A simulation study illustrates small sample size behavior of the bootstrap approximation.     

短面板约束下固定效应空间动态面板模型的估计与模拟

空间动态面板数据(SDPD)模型中被解释变量初值极易带有内生性,采用一般拟极大似然(QML)方法容易造成参数估计偏误,特别是当样本结构为n大T小的时候。鉴于此,本文在一般QML基础上,通过重塑误差项的方差-协方差矩阵,修正拟似然函数表达式,得到修正QML,进而估计短面板下含空间、时间、误差三类关联项的固定效应SDPD模型,基于数值模拟和实例应用检验一般QML与修正QML的估计效果。数值模拟结果表明：修正QML比一般QML更精确、更稳健,均方误差修正率随样本短面板结构的增大而增大。实例应用不仅重新评估环境规制与技术创新之间的空间效应,回归结果也再次证实从数值模拟中得出的结论。     

Asymptotics for REML estimation of spatial covariance parameters

In agricultural field trials, restricted maximum likelihood estimation (REML) of the spatial covariance parameters is often preferred to maximum likelihood. Although it has either been conjectured or assumed that REML estimators are asymptotically Gaussian, conditions under which such asymptotic results hold are clearly needed. This article gives checkable conditions for spatial regression when sampling locations are either on a rectangular grid or are irregularly spaced but satisfy certain growth conditions.