Semiparametric estimation for count data through weighted distributions

This paper is concerned with semiparametric discrete kernel estimators when the unknown count distribution can be considered to have a general weighted Poisson form. The estimator is constructed by multiplying the Poisson estimate with a nonparametric discrete kernel-type estimate of the Poisson weight function. Comparisons are then carried out with the ordinary discrete kernel probability mass function estimators. The Poisson weight function is thus a local multiplicative correction factor, and is considered as the uniform measure to detect departures from the equidispersed Poisson distribution. In this way, the effects of dispersion and zero-proportion with respect to the standard Poisson distribution are also minimized. This method of estimation is also applied to the weighted binomial form for the count distribution having a finite support. The proposed estimators, in addition to being simple, easy-to-implement and effective, also outperform the competing nonparametric and parametric estimators in finite-sample situations. Two examples illustrate this new semiparametric estimation.     

Analysis of robust design experiments with time-dependent ordinal response characteristics: a quality improvement study from the horticulture industry

An approach to the analysis of time-dependent ordinal quality score data from robust design experiments is developed and applied to an experiment from commercial horticultural research, using concepts of product robustness and longevity that are familiar to analysts in engineering research. A two-stage analysis is used to develop models describing the effects of a number of experimental treatments on the rate of post-sales product quality decline. The first stage uses a polynomial function on a transformed scale to approximate the quality decline for an individual experimental unit using derived coefficients and the second stage uses a joint mean and dispersion model to investigate the effects of the experimental treatments on these derived coefficients. The approach, developed specifically for an application in horticulture, is exemplified with data from a trial testing ornamental plants that are subjected to a range of treatments during production and home-life. The results of the analysis show how a number of control and noise factors affect the rate of post-production quality decline. Although the model is used to analyse quality data from a trial on ornamental plants, the approach developed is expected to be more generally applicable to a wide range of other complex production systems.     

Designing and integrating composite networks for monitoring multivariate gaussian pollution fields

Networks of ambient monitoring stations are used to monitor environmental pollution fields such as those for acid rain and air pollution. Such stations provide regular measurements of pollutant concentrations. The networks are established for a variety of purposes at various times so often several stations measuring different subsets of pollutant concentrations can be found in compact geographical regions. The problem of statistically combining these disparate information sources into a single 'network' then arises. Capitalizing on the efficiencies so achieved can then lead to the secondary problem of extending this network. The subject of this paper is a set of 31 air pollution monitoring stations in southern Ontario. Each of these regularly measures a particular subset of ionic sulphate, sulphite, nitrite and ozone. However, this subset varies from station to station. For example only two stations measure all four. Some measure just one. We describe a Bayesian framework for integrating the measurements of these stations to yield a spatial predictive distribution for unmonitored sites and unmeasured concentrations at existing stations. Furthermore we show how this network can be extended by using an entropy maximization criterion. The methods assume that the multivariate response field being measured has a joint Gaussian distribution conditional on its mean and covariance function. A conjugate prior is used for these parameters, some of its hyperparameters being fitted empirically.     

Factors influencing Soay sheep survival

We present a survival analysis of Soay sheep mark recapture and recovery data. Unlike previous conditional analyses, it is not necessary to assume equality of recovery and recapture probabilities; instead these are estimated by maximum likelihood. Male and female sheep are treated separately, with the higher numbers and survival probabilities of the females resulting in a more complex model than that used for the males. In both cases, however, age and time aspects need to be included and there is a strong indication of a reduction in survival for sheep aged 7 years or more. Time variation in survival is related to the size of the population and selected weather variables, by using logistic regression. The size of the population significantly affects the survival probabilities of male and female lambs, and of female sheep aged 7 or more years. March rainfall and a measure of the North Atlantic oscillation are found to influence survival significantly for all age groups considered, for both males and females. Either of these weather variables can be used in a model. Several phenotypic and genotypic individual covariates are also fitted. The only covariate which is found to influence survival significantly is the type of horn of first-year female sheep. There is a substantial variation in the recovery probabilities over time, reflecting in part the increased effort when a population crash was expected. The goodness of fit of the model is checked by using graphical procedures.     

Binary probability maps using a hidden conditional autoregressive Gaussian process with an application to Finnish common toad data

The Finnish common toad data of Heikkinen and Hogmander are reanalysed using an alternative fully Bayesian model that does not require a pseudolikelihood approximation and an alternative prior distribution for the true presence or absence status of toads in each 10 km×10 km square. Markov chain Monte Carlo methods are used to obtain posterior probability estimates of the square-specific presences of the common toad and these are presented as a map. The results are different from those of Heikkinen and Hogmander and we offer an explanation in terms of the prior used for square-specific presence of the toads. We suggest that our approach is more faithful to the data and avoids unnecessary confounding of effects. We demonstrate how to extend our model efficiently with square-specific covariates and illustrate this by introducing deterministic spatial changes.     

Procedures and tables for the construction and selection of chain sampling plans ChSP-4A(c1, c2,)r-Part 4

In this part of the series, the notion of average total inspection (ATI) is introduced to ChSP-4A(c1, c2)r plans, minimizing ATI at a given process average, while protection to the consumer is given in terms of average outgoing quality limit and limiting quality level. The procedure is similar to that of Dodge and Roming.     

农户生计策略识别及其动态转型

利用中国健康与养老追踪调查(CHARLS)2013和2015年两期调查数据,构建收入—资产指数识别了中国农户的生计策略及其跨期动态转型,并从生计阶梯的视角阐释了生计转型对贫困治理的意义。研究发现:(1)农户兼业化和生计多样化程度在不断提高;(2)低收益生计策略退出率和高收益生计策略进入率比较低,生计转型主要发生在中等收益的生计策略中,最终整个生计阶梯演化为“葫芦型”,是对“金字塔”型生计阶梯的一种改善;(3)农户家庭当期的生计策略选择对上一期生计策略具有显著的依赖性。自然资本、人力资本、社会资本和金融资本的积累显著促进了农户生计策略的向高收益流动,而物质资本中的农用固定资产所形成的沉没成本,其捆绑效应阻碍了农户的生计多样化。因此,现阶段扶贫政策应该对农户生计资本提供保障,防止各类外生冲击(包括政策、自然)对农户资产的侵蚀,降低农户生计策略向低收益流动或者单一化的风险。     

Dose‐Response Models Incorporating Aerosol Size Dependency for Francisella tularensis

The effect of bioaerosol size was incorporated into predictive dose‐response models for the effects of inhaled aerosols of Francisella tularensis (the causative agent of tularemia) on rhesus monkeys and guinea pigs with bioaerosol diameters ranging between 1.0 and 24 μm. Aerosol‐size‐dependent models were formulated as modification of the exponential and β‐Poisson dose‐response models and model parameters were estimated using maximum likelihood methods and multiple data sets of quantal dose‐response data for which aerosol sizes of inhaled doses were known. Analysis of F. tularensis dose‐response data was best fit by an exponential dose‐response model with a power function including the particle diameter size substituting for the rate parameter k scaling the applied dose. There were differences in the pathogen's aerosol‐size‐dependence equation and models that better represent the observed dose‐response results than the estimate derived from applying the model developed by the International Commission on Radiological Protection (ICRP, 1994) that relies on differential regional lung deposition for human particle exposure.     

An improved alternative estimation procedure for current population mean in presence of positively and negatively correlated auxiliary variables in two-occasion rotation patterns

Abstract

The present study confirms the influential role of a positively and a negatively correlated auxiliary variables in enhancing the precision of estimates of current population mean in two occasion rotation (successive) sampling. Exponential-type estimators of current population mean have been proposed for three different situations: (i) the information on a positively correlated auxiliary variable is readily available on both occasions (ii) the information on a negatively correlated auxiliary variable is readily available on both occasions and (iii) the information on both positively and negatively correlated auxiliary variables are readily available on both the occasions. The characteristics of the proposed estimators have been explored and their efficacious performances are compared with the natural and recent contemporary estimators. Optimum replacement strategies of the proposed estimation procedures have been formulated. Simulation and empirical studies are carried out to justify the proposition of the proposed estimators and appropriate recommendations have been put forward to the survey practitioners.