Semiparametric Estimation in Copulas with the Same Marginals

We consider semiparametric multivariate data models based on copula representation of the common distribution function. A copula is characterized by a parameter of association and marginal distribution functions. This parameter and the marginal distributions are unknown. In this article, we study the estimator of the parameter of association in copulas with the marginal distribution functions assumed as nuisance parameters restricted by the assumption that the components are identically distributed. Results of this work could be used to construct special kinds of tests of homogeneity for random vectors having dependent components.     

零膨胀计数数据的联合建模及变量选择

零膨胀计数数据破坏了泊松分布的方差-均值关系，可由取值服从泊松分布的数据和取值为零(退化分布)的数据各占一定比例所构成的混合分布所解释。本文基于自适应弹性网技术， 研究了零膨胀计数数据的联合建模及变量选择问题.对于零膨胀泊松分布，引入潜变量，构造出零膨胀泊松模型的完全似然, 其中由零膨胀部分和泊松部分两项组成.考虑到协变量可能存在共线性和稀疏性，通过对似然函数加自适应弹性网惩罚得到目标函数,然后利用EM算法得到回归系数的稀疏估计量，并用贝叶斯信息准则BIC来确定最优调节参数.本文也给出了估计量的大样本性质的理论证明和模拟研究，最后把所提出的方法应用到实际问题中。     

Objective Bayesian Analysis of an Exponential Regression Model with Constrained Parameters Applied to Animal Digestibility

In animal digestibility the proportion of degraded food along the time has usually been modeled as a normal random variable with mean a function of the time and the following three parameters: the proportion of degraded food almost instantaneously, remaining proportion of food to be degraded, and velocity of degradation. The estimation of these parameters has been carried out mainly from a frequentist viewpoint by using the asymptotic distribution of the maximum likelihood estimator. This may give inadmissible estimates, such as values outside of the range of the parameters. This drawback could not appear if a Bayesian approach were adopted. In this article an objective Bayesian analysis is developed and illustrated on real and simulated data.     

Flexible Modeling in the Koziol-Green Model by a Copula Function

In survival analysis, the classical Koziol-Green random censorship model is commonly used to describe informative censoring. Hereby, it is assumed that the distribution of the censoring time is a power of the distribution of the survival time. In this article, we extend this model by assuming a general function between these distributions. We determine this function from a relationship between the observable random variables which is described by a copula family that depends on an unknown parameter θ. For this setting, we develop a semi-parametric estimator for the distribution of the survival time in which we propose a pseudo-likelihood estimator for the copula parameter θ. As results, we show first the consistency and asymptotic normality of the estimator for θ. Afterwards, we prove the weak convergence of the process associated to the semi-parametric distribution estimator. Furthermore, we investigate the finite sample performance of these estimators through a simulation study and finally apply it to a practical data set on survival with malignant melanoma.     

Kernel density estimation under negative superadditive dependence and its application for real data

In this paper, the kernel density estimator for negatively superadditive dependent random variables is studied. The exponential inequalities and the exponential rate for the kernel estimator of density function with a uniform version, over compact sets are investigated. Also, the optimal bandwidth rate of the estimator is obtained using mean integrated squared error. The results are generalized and used to improve the ones obtained for the case of associated sequences. As an application, FGM sequences that fulfil our assumptions are investigated. Also, the convergence rate of the kernel density estimator is illustrated via a simulation study. Moreover, a real data analysis is presented.     

On the asymptotic normality of the kernel estimators of the density function and its derivatives under censoring

In this paper, we study asymptotic normality of the kernel estimators of the density function and its derivatives as well as the mode in the randomly right censorship model. The mode estimator is defined as the random variable that maximizes the kernel density estimator. Our results are stated under some suitable conditions upon the kernel function, the smoothing parameter and both distributions functions that appear in this model. Here, the Kaplan–Meier estimator of the distribution function is used to build the estimates. We carry out a simulation study which shows how good the normality works.     

Semiparametric Approach to a Random Effects Quantile Regression Model

We consider a random effects quantile regression analysis of clustered data and propose a semiparametric approach using empirical likelihood. The random regression coefficients are assumed independent with a common mean, following parametrically specified distributions. The common mean corresponds to the population-average effects of explanatory variables on the conditional quantile of interest, while the random coefficients represent cluster specific deviations in the covariate effects. We formulate the estimation of the random coefficients as an estimating equations problem and use empirical likelihood to incorporate the parametric likelihood of the random coefficients. A likelihood-like statistical criterion function is yield, which we show is asymptotically concave in a neighborhood of the true parameter value and motivates its maximizer as a natural estimator. We use Markov Chain Monte Carlo (MCMC) samplers in the Bayesian framework, and propose the resulting quasi-posterior mean as an estimator. We show that the proposed estimator of the population-level parameter is asymptotically normal and the estimators of the random coefficients are shrunk toward the population-level parameter in the first order asymptotic sense. These asymptotic results do not require Gaussian random effects, and the empirical likelihood based likelihood-like criterion function is free of parameters related to the error densities. This makes the proposed approach both flexible and computationally simple. We illustrate the methodology with two real data examples.     

A general approach to heteroscedastic linear regression

Our article presents a general treatment of the linear regression model, in which the error distribution is modelled nonparametrically and the error variances may be heteroscedastic, thus eliminating the need to transform the dependent variable in many data sets. The mean and variance components of the model may be either parametric or nonparametric, with parsimony achieved through variable selection and model averaging. A Bayesian approach is used for inference with priors that are data-based so that estimation can be carried out automatically with minimal input by the user. A Dirichlet process mixture prior is used to model the error distribution nonparametrically; when there are no regressors in the model, the method reduces to Bayesian density estimation, and we show that in this case the estimator compares favourably with a well-regarded plug-in density estimator. We also consider a method for checking the fit of the full model. The methodology is applied to a number of simulated and real examples and is shown to work well.     

Bayesian inference using product of spacings function for Progressive hybrid Type-I censoring scheme

This article is devoted to the development of product of spacings estimator for a Progressive hybrid Type-I censoring scheme with binomial removals. The experimental units are assumed to follow inverse Lindley distribution. We propose a Bayes estimator of associated scale parameter based on the product of spacings function and simultaneously compare it with that obtained under a usual Bayesian estimation procedure. The estimators are obtained under the squared error loss function along with corresponding HP intervals evaluated by using the Markov chain Monte-Carlo technique. The classical product of spacings estimator has also been derived and compared with the maximum likelihood estimator in addition to 95% average asymptotic confidence intervals. The applicability of the proposed methods is demonstrated by analysing a real data of guinea pigs affected with tuberculosis for the considered censoring scheme.     

Statistical analysis for Kumaraswamy’s distribution based on record data

In this paper we review some results that have been derived on record values for some well known probability density functions and based on m records from Kumaraswamy’s distribution we obtain estimators for the two parameters and the future sth record value. These estimates are derived using the maximum likelihood and Bayesian approaches. In the Bayesian approach, the two parameters are assumed to be random variables and estimators for the parameters and for the future sth record value are obtained, when we have observed m past record values, using the well known squared error loss (SEL) function and a linear exponential (LINEX) loss function. The findings are illustrated with actual and computer generated data.     

Nonparametric estimation of the survival function from censored data based on relative risk function

We present a nonparametric estimator for distribution function under random censorship from the right. Our approach is based on estimating the relative risk function, and the resulting estimator is closely related with the Kaplan-Meier’s product-limit and Breslow’s exponential hazard estimators. We also consider the general proportional hazards model when the competing risks are variable censored from the left.     

指数族分布中参数的经验贝叶斯估计

指数族分布是一类应用广泛的分布类,包括了泊松分布、Gamma分布、Beta分布、二项分布等常见分布.在非寿险中,索赔额或索赔次数过程常常被假定服从指数族分布,由于风险的非齐次性,指数族分布中的参数θ也为随机变量,假定服从指数族共轭先验分布.此时风险参数的估计落入了Bayes框架,风险参数θ的Bayes估计被表达“信度”形式.然而,在实际运用中,由于先验分布与样本分布中仍然含有结构参数,根据样本的边际分布的似然函数估计结构参数,从而获得风险参数的经验Bayes估计,最后证明了该经验Bayes估计是渐近最优的.     

A semiparametric estimator of the distribution function of a variable measured with error

The estimation of the distribution functon of a random variable X measured with error is studied. Let the i-th observation on X be denoted by Y_iX_i+ε_i where ε_i is the measuremen error. Let {Y_i} (i=1,2,…,n) be a sample of independent observations. It is assumed that {X_i} and {∈_i} are mutually independent and each is identically distributed. As is standard in the literature for this problem, the distribution of e is assumed known in the development of the methodology. In practice, the measurement error distribution is estimated from replicate observations.

The proposed semiparametric estimator is derived by estimating the quantises of X on a set of n transformed V-values and smoothing the estimated quantiles using a spline function. The number of parameters of the spline function is determined by the data with a simple criterion, such as AIC. In a simulation study, the semiparametric estimator dominates an optimal kernel estimator and a normal mixture estimator for a wide class of densities.

The proposed estimator is applied to estimate the distribution function of the mean pH value in a field plot. The density function of the measurement error is estimated from repeated measurements of the pH values in a plot, and is treated as known for the estimation of the distribution function of the mean pH value.     

Estimation and Confidence Intervals for Clock Offset in Networks with Bivariate Exponential Delays

There is a large and increasing literature on statistical modeling-based estimation of the offset between two clocks. Recent work has focused on the construction of confidence intervals for offset. However, in most of this work it has been assumed that the network delays that occur during the synchronization process are independent. The network delays are often modeled as independent exponential random variables. Thus, we introduce the use of a bivariate exponential distribution to capture the anticipated correlation between the network delays and derive a maximum likelihood estimator and a confidence interval procedure for the offset parameter. We then illustrate how use of the independent model for network delays can lead to improper inference about the offset parameter.     

Bayesian inference for the generalized exponential distribution

The two-parameter generalized exponential (GE) distribution was introduced by Gupta and Kundu [Gupta, R.D. and Kundu, D., 1999, Generalized exponential distribution. Australian and New Zealand Journal of Statistics, 41(2), 173–188.]. It was observed that the GE can be used in situations where a skewed distribution for a nonnegative random variable is needed. In this article, the Bayesian estimation and prediction for the GE distribution, using informative priors, have been considered. Importance sampling is used to estimate the parameters, as well as the reliability function, and the Gibbs and Metropolis samplers data sets are used to predict the behavior of further observations from the distribution. Two data sets are used to illustrate the Bayesian procedure.     

A note on Bayesian nonparametric regression function estimation

In this note the problem of nonparametric regression function estimation in a random design regression model with Gaussian errors is considered from the Bayesian perspective. It is assumed that the regression function belongs to a class of functions with a known degree of smoothness. A prior distribution on the given class can be induced by a prior on the coefficients in a series expansion of the regression function through an orthonormal system. The rate of convergence of the resulting posterior distribution is employed to provide a measure of the accuracy of the Bayesian estimation procedure defined by the posterior expected regression function. We show that the Bayes’ estimator achieves the optimal minimax rate of convergence under mean integrated squared error over the involved class of regression functions, thus being comparable to other popular frequentist regression estimators.     

Logspline Deconvolution in Besov Space

ABSTRACT. In this paper we consider logspline density estimation for random variables which are contaminated with random noise. In the logspline density estimation for data without noise, the logarithm of an unknown density function is estimated by a polynomial spline, the unknown parameters of which are given by maximum likelihood. When noise is present, B-splines and the Fourier inversion formula are used to construct the logspline density estimator of the unknown density function. Rates of convergence are established when the log-density function is assumed to be in a Besov space. It is shown that convergence rates depend on the smoothness of the density function and the decay rate of the characteristic function of the noise. Simulated data are used to show the finite-sample performance of inference based on the logspline density estimation.     

On estimation of the PDF and CDF of the Lindley distribution

This article addresses two methods of estimation of the probability density function (PDF) and cumulative distribution function (CDF) for the Lindley distribution. Following estimation methods are considered: uniformly minimum variance unbiased estimator (UMVUE) and maximum likelihood estimator (MLE). Since the Lindley distribution is more flexible than the exponential distribution, the same estimators have been found out for the exponential distribution and compared. Monte Carlo simulations and a real data analysis are performed to compare the performances of the proposed methods of estimation.     

A note on estimation in a minimum-type model with an exponential factor

The problem of estimating the parameter Q appearing in the distribution function of a continuous random variable T=min{T₁,T₂} is considered, when T₁,T₂ are non-negative independent random variables such that T₁ has an exponential distribution with scale parameter θ and T₂ has a possibly defective distribution function G₂(t) that is G₂(∞) < 1. It is shown that the estimator proposed by Gertsbach [(1967) Theory of Probability and Its Applications, 12 is weekly consistent and asymptotically normal. The merit of this estimator is that the above properties do not depend even on the form of G₂(t) except that G₂(t) and its derivative vanish at zero.     

Estimation procedures with delayed observations

The statistical model is considered in which the collection of data from several independent populations is available only at random times determined by order statistics of lifetimes of a given number of objects. Each of the populations is distributed according to a general multiparameter exponential family. The problem is to estimate the mean value vector parameter of the multiparameter exponential family of distributions of the forthcoming observations. Under the loss function involving a weighted squared error loss, the cost proportional to the events appeared and a cost of observing the process, a class of optimal sequential procedures is established. The procedures are derived in two situations: when the lifetime distribution is completely known and in the case when it is unknown but assumed to belong to an exponential subfamily with an unknown failure rate parameter.