Communication in Statistics-Theory and methods improved GQL estimation method for the generalised BINMA(1) model

In a recent research, the quasi-likelihood estimation methodology was developed to estimate the regression effects in the Generalized BINMA(1) (GBINMA(1)) process. The method provides consistent parameter estimates but, in the intermediate computations, moment estimating equations were used to estimate the serial- and cross-correlation parameters. This procedure may not result optimal parameter estimates, in particular, for the regression effects. This paper provides an alternative simpler GBINMA(1) process based on multivariate thinning properties where the main effects are estimated via a robust generalized quasi-likelihood (GQL) estimation approach. The two techniques are compared through some simulation experiments. A real-life data application is studied.     

Cytometry inference through adaptive atomic deconvolution

In this paper, we consider a statistical estimation problem known as atomic deconvolution. Introduced in reliability, this model has a direct application when considering biological data produced by flow cytometers. From a statistical point of view, we aim at inferring the percentage of cells expressing the selected molecule and the probability distribution function associated with its fluorescence emission. We propose here an adaptive estimation procedure based on a previous deconvolution procedure introduced by Es, Gugushvili, and Spreij [(2008), ‘Deconvolution for an atomic distribution’, Electronic Journal of Statistics, 2, 265–297] and Gugushvili, Es, and Spreij [(2011), ‘Deconvolution for an atomic distribution: rates of convergence’, Journal of Nonparametric Statistics, 23, 1003–1029]. For both estimating the mixing parameter and the mixing density automatically, we use the Lepskii method based on the optimal choice of a bandwidth using a bias-variance decomposition. We then derive some convergence rates that are shown to be minimax optimal (up to some log terms) in Sobolev classes. Finally, we apply our algorithm on the simulated and real biological data.     

A new regression model for bimodal data and applications in agriculture

We define the odd log-logistic exponential Gaussian regression with two systematic components, which extends the heteroscedastic Gaussian regression and it is suitable for bimodal data quite common in the agriculture area. We estimate the parameters by the method of maximum likelihood. Some simulations indicate that the maximum-likelihood estimators are accurate. The model assumptions are checked through case deletion and quantile residuals. The usefulness of the new regression model is illustrated by means of three real data sets in different areas of agriculture, where the data present bimodality.     

A comparison study on modeling of clustered and overdispersed count data for multiple comparisons

Data collected in various scientific fields are count data. One way to analyze such data is to compare the individual levels of the factor treatment using multiple comparisons. However, the measured individuals are often clustered – e.g. according to litter or rearing. This must be considered when estimating the parameters by a repeated measurement model. In addition, ignoring the overdispersion to which count data is prone leads to an increase of the type one error rate. We carry out simulation studies using several different data settings and compare different multiple contrast tests with parameter estimates from generalized estimation equations and generalized linear mixed models in order to observe coverage and rejection probabilities. We generate overdispersed, clustered count data in small samples as can be observed in many biological settings. We have found that the generalized estimation equations outperform generalized linear mixed models if the variance-sandwich estimator is correctly specified. Furthermore, generalized linear mixed models show problems with the convergence rate under certain data settings, but there are model implementations with lower implications exists. Finally, we use an example of genetic data to demonstrate the application of the multiple contrast test and the problems of ignoring strong overdispersion.     

非线性面板数据聚类方法研究

对于一类变量非线性相关的面板数据,现有的基于线性算法的面板数据聚类方法并不能准确地度量样本间的相似性,且聚类结果的可解释性低。综合考虑变量非线性相关问题及聚类结果可解释性问题,提出一种非线性面板数据的聚类方法,通过非线性核主成分算法实现对样本相似性的测度,并基于混合高斯模型进行样本概率聚类,实证表明该方法的有效性及其对聚类结果的可解释性有所提高。     

Testing for bioequivalence of highly variable drugs from TR‐RT crossover designs with heterogeneous residual variances

Traditional bioavailability studies assess average bioequivalence (ABE) between the test (T) and reference (R) products under the crossover design with TR and RT sequences. With highly variable (HV) drugs whose intrasubject coefficient of variation in pharmacokinetic measures is 30% or greater, assertion of ABE becomes difficult due to the large sample sizes needed to achieve adequate power. In 2011, the FDA adopted a more relaxed, yet complex, ABE criterion and supplied a procedure to assess this criterion exclusively under TRR‐RTR‐RRT and TRTR‐RTRT designs. However, designs with more than 2 periods are not always feasible. This present work investigates how to evaluate HV drugs under TR‐RT designs. A mixed model with heterogeneous residual variances is used to fit data from TR‐RT designs. Under the assumption of zero subject‐by‐formulation interaction, this basic model is comparable to the FDA‐recommended model for TRR‐RTR‐RRT and TRTR‐RTRT designs, suggesting the conceptual plausibility of our approach. To overcome the distributional dependency among summary statistics of model parameters, we develop statistical tests via the generalized pivotal quantity (GPQ). A real‐world data example is given to illustrate the utility of the resulting procedures. Our simulation study identifies a GPQ‐based testing procedure that evaluates HV drugs under practical TR‐RT designs with desirable type I error rate and reasonable power. In comparison to the FDA's approach, this GPQ‐based procedure gives similar performance when the product's intersubject standard deviation is low (≤0.4) and is most useful when practical considerations restrict the crossover design to 2 periods.     

我国经济欠发达地区农村社会养老保险制度的构建

与发达地区相比,我国经济欠发达地区农村养老保险发展十分缓慢,覆盖面极其有限。实际上,经济欠发达地区的农民最需要社会养老保险,因为他们的基本生活需求没有保证,这些地区未富先老特征更为明显,传统的保障功能在弱化。然而,源于农民收入低,增长慢,支付保险费的能力弱;集体经济薄弱,无法为农民养老保险提供资金支持;地方财力单薄,难以全面提供以政府为主的养老基金融资;经济欠发达地区农民获得社会养老保险能力又最弱。为此,应将经济欠发达地区农村社会养老保险作为建设覆盖城乡居民社会保障体系的推进重点,变政府扶持为财政拨款,并向经济欠发达地区倾斜;大力进行宣传教育,让农民充分了解参加社会养老保险的必要性。     

Bayesian robustness of the compound Poisson distribution under bidimensional prior: an application to the collective risk model

The distribution of the aggregate claims in one year plays an important role in Actuarial Statistics for computing, for example, insurance premiums when both the number and size of the claims must be implemented into the model. When the number of claims follows a Poisson distribution the aggregated distribution is called the compound Poisson distribution. In this article we assume that the claim size follows an exponential distribution and later we make an extensive study of this model by assuming a bidimensional prior distribution for the parameters of the Poisson and exponential distribution with marginal gamma. This study carries us to obtain expressions for net premiums, marginal and posterior distributions in terms of some well-known special functions used in statistics. Later, a Bayesian robustness study of this model is made. Bayesian robustness on bidimensional models was deeply treated in the 1990s, producing numerous results, but few applications dealing with this problem can be found in the literature.     

道德的经济学——亲缘选择与广义效用理论

早期功利主义者认为效用不仅包括人们获得的物质利益,还包括精神上的追求和道德上的满足感。然而经济学数学化和公理化的趋势使经济学家更愿意把效用解释为物质利益上的满足,这极大地限制了经济学的解释力度。广义效用理论是在现代经济学的框架下对效用含义的扩展。通过对当前广义效用理论的回顾,结合生物学中的亲缘选择理论,提出了一个广义效用模型,并对模型的结果进行了分析。     

Copula-based regression models: A survey

In this review paper we collect several results about copula-based models, especially concerning regression models, by focusing on some insurance applications.