Asymptotic behavior of RA-estimates in autoregressive 2D processes

In this work we study the asymptotic behavior of a robust class of estimators of the coefficient of a AR-2D process. We establish the precise conditions for the consistency and asymptotic normality of the RA estimator. The AR-2D model has many applications in image modeling and statistical image processing, therefore the relevance of knowing such properties. The adequacy of the AR-2D model is analyzed with real images; we also show the impact of contamination and the capability of the RA estimator to produce useful results even in the presence of spurious data.     

On univariate and bivariate generalized gamma convolutions

This paper has two parts. In the first part some results for generalized gamma convolutions (GGCs) are reviewed. A GGC is a limit distribution for sums of independent gamma variables. In the second part, bivariate gamma distributions and bivariate GGCs are considered. New bivariate gamma distributions are derived from shot-noise models. The remarkable property hyperbolic complete monotonicity (HCM) for a function is considered both in the univariate case and in the bivariate case.     

Fully exponential Laplace approximations for the joint modelling of survival and longitudinal data

Summary.  A common objective in longitudinal studies is the joint modelling of a longitudinal response with a time-to-event outcome. Random effects are typically used in the joint modelling framework to explain the interrelationships between these two processes. However, estimation in the presence of random effects involves intractable integrals requiring numerical integration. We propose a new computational approach for fitting such models that is based on the Laplace method for integrals that makes the consideration of high dimensional random-effects structures feasible. Contrary to the standard Laplace approximation, our method requires much fewer repeated measurements per individual to produce reliable results.     

Skew-normal semiparametric varying coefficient model and score test

In this paper, we consider inference aspects of skew-normal semiparametric varying coefficient models which provide a useful extension of the normal regression models. The maximum likelihood estimation based on B-spline is proposed. Further, we discuss the score test for homogeneity of the variance in skew-normal semiparametric varying coefficient models. Their asymptotical properties are investigated. Some simulated examples are used to examine our proposed methods.     

A simple approach for varying-coefficient model selection

In varying-coefficient models, an important question is to determine whether some of the varying coefficients are actually invariant coefficients. This article proposes a penalized likelihood method in the framework of the smoothing spline ANOVA models, with a penalty designed toward the goal of automatically distinguishing varying coefficients and those which are not varying. Unlike the stepwise procedure, the method simultaneously quantifies and estimates the coefficients. An efficient algorithm is given and ways of choosing the smoothing parameters are discussed. Simulation results and an analysis on the Boston housing data illustrate the usefulness of the method. The proposed approach is further extended to longitudinal data analysis.     

Sampling methods and sensitivity analysis for large parameter sets

Models with large parameter (i.e., hundreds or thousands of parameters) often behave as if they depend upon only a few parameters, with the rest having comparatively little influence. One challenge of sensitivity analysis with such models is screening parameters to identify the influential ones, and then characterizing their influences.

Large models often require significant computing resources to evaluate their output, and so a good screening mechanism should be efficient: it should minimize the number of times a model must be exercised. This paper describes an efficient procedure to perform sensitivity analysis on deterministic models with specified ranges or probability distributions for each parameter.

It is based on repeated exercising of the model, which can be treated as a black box. Statistical checks can ensure that the screening identified parameters that account for the bulk of the model variation. Subsequent sensitivity analysis can use the screening information to reduce the investment required to characterize the influence of influential and other parameters.

The procedure exploits simplifications in the dependence of a model output on model inputs. It works best where a small number of parameters are much more influential than all the rest. The method is much more sensitive to the number of influential parameters than to the total number of parameters. It is most effective when linear or quadratic effects dominate higher order effects and complex interactions.

The paper presents a set of M athematica functions that can be used to create a variety of types of experimental designs useful for sensitivity analysis, including simple random, latin hypercube and fractional factorial sampling. Each sampling method can use discretization, folding, grouping and replication to create composite designs. These techniques have beencombined in a composite approach called Iterated Fractional Factorial Design (IFFD).

The procedure is applied to model of nuclear fuel waste disposal, and to simplified example models to demonstrate the concepts involved.     

Flash Flood Risks and Warning Decisions: A Mental Models Study of Forecasters,Public Officials,and Media Broadcasters in Boulder,Colorado

Timely warning communication and decision making are critical for reducing harm from flash flooding. To help understand and improve extreme weather risk communication and management, this study uses a mental models research approach to investigate the flash flood warning system and its risk decision context. Data were collected in the Boulder, Colorado area from mental models interviews with forecasters, public officials, and media broadcasters, who each make important interacting decisions in the warning system, and from a group modeling session with forecasters. Analysis of the data informed development of a decision‐focused model of the flash flood warning system that integrates the professionals’ perspectives. Comparative analysis of individual and group data with this model characterizes how these professionals conceptualize flash flood risks and associated uncertainty; create and disseminate flash flood warning information; and perceive how warning information is (and should be) used in their own and others’ decisions. The analysis indicates that warning system functioning would benefit from professionals developing a clearer, shared understanding of flash flood risks and the warning system, across their areas of expertise and job roles. Given the challenges in risk communication and decision making for complex, rapidly evolving hazards such as flash floods, another priority is development of improved warning content to help members of the public protect themselves when needed. Also important is professional communication with members of the public about allocation of responsibilities for managing flash flood risks, as well as improved system‐wide management of uncertainty in decisions.     

Analysis of an EOQ inventory model with partial backordering and non-linear unit holding cost

In this paper, an economic order quantity inventory model is analyzed, considering that the unit cumulative holding cost has two significant components: a fixed cost which represents the cost of accommodating the item in the warehouse and a variable cost given by a potential function of the length of time over which the item is held in stock. Shortages are allowed and, during the stockout period, only a fraction of demand is partially backordered. The backordering cost includes a fixed cost and a cost linearly dependent on the length of time for which backorder exists. A solution procedure is developed for determining the optimal inventory policy. Moreover, to illustrate the effects of some parameters on the optimal policy and the minimum total inventory cost, a numerical study is developed.     

时滞交联系统的H_∞滤波设计

基于T-S模糊模型方法,针对一类时滞交联系统,提出了一种集中H∞滤波设计新方案.通过构造新的模糊线积分李雅普诺夫泛函,给出了整体滤波误差系统时滞相关的集中H∞滤波设计LMI条件,降低了性能指标的保守性.最后通过仿真例子说明了所提方法的有效性和优越性.     

Posterior Predictive p-values in Bayesian Hierarchical Models

Abstract.  The present work focuses on extensions of the posterior predictive p -value (ppp-value) for models with hierarchical structure, designed for testing assumptions made on underlying processes. The ppp-values are popular as tools for model criticism, yet their lack of a common interpretation limit their practical use. We discuss different extensions of ppp-values to hierarchical models, allowing for discrepancy measures that can be used for checking properties of the model at all stages. Through analytical derivations and simulation studies on simple models, we show that similar to the standard ppp-values, these extensions are typically far from uniformly distributed under the model assumptions and can give poor power in a hypothesis testing framework. We propose a calibration of the p -values, making the resulting calibrated p -values uniformly distributed under the model conditions. Illustrations are made through a real example of multinomial regression to age distributions of fish.