Test of Independence for Baker’s Bivariate Distributions

Baker (2008 Baker, R. (2008). An order-statistics-based method for constructing multivariate distributions with fixed marginals. Journal of Multivariate Analysis 99: 2312–2327.[Crossref], [Web of Science ®] , [Google Scholar]) introduced a new method for constructing multivariate distributions with given marginals based on order statistics. In this paper, we provide a test of independence for a pair of absolutely continuous random variables (X, Y) jointly distributed according to Baker’s bivariate distributions. Our purpose is to test the hypothesis that X and Y are independent versus the alternative that X and Y are positively (negatively) quadrant dependent. The asymptotic distribution of the proposed test statistic is investigated. Also, the powers of the proposed test and the class of distribution-free tests proposed by Kochar and Gupta (1987 Kochar, S. G., Gupta, R. P. (1987). Competitors of Kendall-tau test for testing independence against positive quadrant dependence. Biometrika 74(3): 664–666.[Crossref], [Web of Science ®] , [Google Scholar]) are compared empirically via a simulation study.     

Complete convergence for weighted sums of END random variables and its application to nonparametric regression models

In this article, the complete convergence for weighted sums of extended negatively dependent (END, for short) random variables is investigated. Some sufficient conditions for the complete convergence are provided. In addition, the Marcinkiewicz–Zygmund type strong law of large numbers for weighted sums of END random variables is obtained. The results obtained in the article generalise and improve the corresponding one of Wang et al. [(2014b), ‘On Complete Convergence for an Extended Negatively Dependent Sequence’, Communications in Statistics-Theory and Methods, 43, 2923–2937]. As an application, the complete consistency for the estimator of nonparametric regression model is established.     

A Comparison of Bayesian Models of Heteroscedasticity in Nested Normal Data

We consider the fitting of a Bayesian model to grouped data in which observations are assumed normally distributed around group means that are themselves normally distributed, and consider several alternatives for accommodating the possibility of heteroscedasticity within the data. We consider the case where the underlying distribution of the variances is unknown, and investigate several candidate prior distributions for those variances. In each case, the parameters of the candidate priors (the hyperparameters) are themselves given uninformative priors (hyperpriors). The most mathematically convenient model for the group variances is to assign them inverse gamma distributed priors, the inverse gamma distribution being the conjugate prior distribution for the unknown variance of a normal population. We demonstrate that for a wide class of underlying distributions of the group variances, a model that assigns the variances an inverse gamma-distributed prior displays favorable goodness-of-fit properties relative to other candidate priors, and hence may be used as standard for modeling such data. This allows us to take advantage of the elegant mathematical property of prior conjugacy in a wide variety of contexts without compromising model fitness. We test our findings on nine real world publicly available datasets from different domains, and on a wide range of artificially generated datasets.     

A flexible multi‐domain test with adaptive weights and its application to clinical trials

The design of a clinical trial is often complicated by the multi‐systemic nature of the disease; a single endpoint often cannot capture the spectrum of potential therapeutic benefits. Multi‐domain outcomes which take into account patient heterogeneity of disease presentation through measurements of multiple symptom/functional domains are an attractive alternative to a single endpoint. A multi‐domain test with adaptive weights is proposed to synthesize the evidence of treatment efficacy over numerous disease domains. The test is a weighted sum of domain‐specific test statistics with weights selected adaptively via a data‐driven algorithm. The null distribution of the test statistic is constructed empirically through resampling and does not require estimation of the covariance structure of domain‐specific test statistics. Simulations show that the proposed test controls the type I error rate, and has increased power over other methods such as the O'Brien and Wei‐Lachin tests in scenarios reflective of clinical trial settings. Data from a clinical trial in a rare lysosomal storage disorder were used to illustrate the properties of the proposed test. As a strategy of combining marginal test statistics, the proposed test is flexible and readily applicable to a variety of clinical trial scenarios.     

Exact Inference for a Simple Step-stress Model with Generalized Type-I Hybrid Censored Data from the Exponential Distribution

In this article, the simple step-stress model is considered based on generalized Type-I hybrid censored data from the exponential distribution. The maximum likelihood estimators (MLEs) of the unknown parameters are derived assuming a cumulative exposure model. We then derive the exact distributions of the MLEs of the parameters using conditional moment generating functions. The Bayesian estimators of the parameters are derived and then compared with the MLEs. We also derive confidence intervals for the parameters using these exact distributions, asymptotic distributions of the MLEs, Bayesian, and the parametric bootstrap methods. The problem of determining the optimal stress-changing point is discussed and the MLEs of the pth quantile and reliability functions at the use condition are obtained. Finally, Monte Carlo simulation and some numerical results are presented for illustrating all the inferential methods developed here.     

Numerical computation of bridge probabilities

This paper contains an analysis of the problem of computing the joint probability density of the honor card point count in each of four hands in the game of bridge. Efficient representation of the data is considered. Computational algorithms are given for dealing with a compressed form of the density. From the joint density, the densities of the point count in the best and worst hands are obtained. Also obtained is the conditional distribution that a partnership has m points given that one of the partners has n1 points.     

Cross-country risk quantification of extreme wildfires in Mediterranean Europe

We estimate the country-level risk of extreme wildfires defined by burned area (BA) for Mediterranean Europe and carry out a cross-country comparison. To this end, we avail of the European Forest Fire Information System (EFFIS) geospatial data from 2006 to 2019 to perform an extreme value analysis. More specifically, we apply a point process characterization of wildfire extremes using maximum likelihood estimation. By modeling covariates, we also evaluate potential trends and correlations with commonly known factors that drive or affect wildfire occurrence, such as the Fire Weather Index as a proxy for meteorological conditions, population density, land cover type, and seasonality. We find that the highest risk of extreme wildfires is in Portugal (PT), followed by Greece (GR), Spain (ES), and Italy (IT) with a 10-year BA return level of 50'338 ha, 33'242 ha, 25'165 ha, and 8'966 ha, respectively. Coupling our results with existing estimates of the monetary impact of large wildfires suggests expected losses of 162–439 million € (PT), 81–219 million € (ES), 41–290 million € (GR), and 18–78 million € (IT) for such 10-year return period events.

SUMMARY

We model the risk of extreme wildfires for Italy, Greece, Portugal, and Spain in form of burned area return levels, compare them, and estimate expected losses.     

金融资产的VaR和CVaR风险的优良估计

本文利用统计理论的优良点估计方法来估计金融市场风险的VaR和CVaR,既可避开现有方法中大量的模拟计算和参数估计等工作,又可提高估算精度.在资产-正态模型下,根据不同的风险估计要求,对金融资产的这两种风险分别提供了三种优良估计,即一致最小方差无偏估计,最佳线性次序统计量无偏估计,最佳线性次序统计量同变估计,并提供了实证分析.     

A Robust Approach to Risk Assessment Based on Species Sensitivity Distributions

The guidelines for setting environmental quality standards are increasingly based on probabilistic risk assessment due to a growing general awareness of the need for probabilistic procedures. One of the commonly used tools in probabilistic risk assessment is the species sensitivity distribution (SSD), which represents the proportion of species affected belonging to a biological assemblage as a function of exposure to a specific toxicant. Our focus is on the inverse use of the SSD curve with the aim of estimating the concentration, HCp, of a toxic compound that is hazardous to p% of the biological community under study. Toward this end, we propose the use of robust statistical methods in order to take into account the presence of outliers or apparent skew in the data, which may occur without any ecological basis. A robust approach exploits the full neighborhood of a parametric model, enabling the analyst to account for the typical real‐world deviations from ideal models. We examine two classic HCp estimation approaches and consider robust versions of these estimators. In addition, we also use data transformations in conjunction with robust estimation methods in case of heteroscedasticity. Different scenarios using real data sets as well as simulated data are presented in order to illustrate and compare the proposed approaches. These scenarios illustrate that the use of robust estimation methods enhances HCp estimation.