The Extended Skew-Exponential Power Distribution and Its Derivation

We consider an extended family of asymmetric univariate distributions generated using a symmetric density, f, and the cumulative distribution function, G, of a symmetric distribution, which depends on two real-valued parameters λ and β and is such that when β = 0 it includes the entire class of distributions with densities of the form g(z | λ) = 2 G(λ z) f(z). A key element in the construction of random variables distributed according to the family is that they can be represented stochastically as the product of two random variables. From this representation we can readily derive theoretical properties, easy-to-implement simulation schemes, as well as extensions to the multivariate case and an explicit procedure for obtaining the moments. We give special attention to the extended skew-exponential power distribution. We derive its information matrix in order to obtain the asymptotic covariance matrix of the maximum likelihood estimators. Finally, an application to a real data set is reported, which shows that the extended skew-exponential power model can provide a better fit than the skew-exponential power distribution.     

Divergence versus decision P-values: A distinction worth making in theory and keeping in practice: Or,how divergence P-values measure evidence even when decision P-values do not

There are two distinct definitions of “P-value” for evaluating a proposed hypothesis or model for the process generating an observed dataset. The original definition starts with a measure of the divergence of the dataset from what was expected under the model, such as a sum of squares or a deviance statistic. A P-value is then the ordinal location of the measure in a reference distribution computed from the model and the data, and is treated as a unit-scaled index of compatibility between the data and the model. In the other definition, a P-value is a random variable on the unit interval whose realizations can be compared to a cutoff α to generate a decision rule with known error rates under the model and specific alternatives. It is commonly assumed that realizations of such decision P-values always correspond to divergence P-values. But this need not be so: Decision P-values can violate intuitive single-sample coherence criteria where divergence P-values do not. It is thus argued that divergence and decision P-values should be carefully distinguished in teaching, and that divergence P-values are the relevant choice when the analysis goal is to summarize evidence rather than implement a decision rule.     

A Novel Bayesian Parameter Mapping Method for Estimating the Parameters of an Underlying Scientific Model

Population-parameter mapping (PPM) is a method for estimating the parameters of latent scientific models that describe the statistical likelihood function. The PPM method involves a Bayesian inference in terms of the statistical parameters and the mapping from the statistical parameter space to the parameter space of the latent scientific parameters, and obtains a model coherence estimate, P(coh). The P(coh) statistic can be valuable for designing experiments, comparing competing models, and can be helpful in redesigning flawed models. Examples are provided where greater estimation precision was found for small sample sizes for the PPM point estimates relative to the maximum likelihood estimator (MLE).     

On the Sum of Symmetric Random Variables

Let U and V be two symmetric (about zero) random variables with U + V symmetric about C; here C is a constant. It is easy to see that if U and V are mutually independent, or if both U and V satisfy the weak law of large numbers, then C = 0. So, intuitively, we would suspect that C = 0 in general. However, we show that there exist two random variables U and V symmetric about 0 with U + V symmetric about C ≠ 0 The example given is closely related to one given by Alejandro D. De Acosta in another context.     

The Beta Product Distribution with Complex Parameters

The family consisting of the distributions of products of two independent beta variables is extended to include cases where some of the parameters are not positive but negative or complex. This “beta product” distribution is expressible as a Meijer G function. An example (from risk theory) where such a distribution arises is given: an infinite sum of products of independent random variables is shown to have a distribution that is the product convolution of a complex-parameter beta product and an independent exponential. The distribution of the infinite sum is a new explicit solution of the stochastic equation X = (in law) B(X + C). Characterizations of some G distributions are also proved.     

Nonparametric estimation of the multivariate survivor function: the multivariate Kaplan–Meier estimator

The Dabrowska (Ann Stat 16:1475–1489, 1988) product integral representation of the multivariate survivor function is extended, leading to a nonparametric survivor function estimator for an arbitrary number of failure time variates that has a simple recursive formula for its calculation. Empirical process methods are used to sketch proofs for this estimator’s strong consistency and weak convergence properties. Summary measures of pairwise and higher-order dependencies are also defined and nonparametrically estimated. Simulation evaluation is given for the special case of three failure time variates.     

Log-gamma-generated families of distributions

We introduce two new general families of continuous distributions, generated by a distribution F and two positive real parameters α and β which control the skewness and tail weight of the distribution. The construction is motivated by the distribution of k-record statistics and can be derived by applying the inverse probability integral transformation to the log-gamma distribution. The introduced families are suitable for modelling the data with a significantly skewed and heavy-tailed distribution. Various properties of the introduced families are studied and a number of estimations and data fitness on real data are given to illustrate the results.     

The Rate of Convergence of some Asymptotic Expansions for Distribution Approximations via an Esseen Type Estimate

Some asymptotic expansions not necessarily related to the central limit theorem are studied. We first observe that the smoothing inequality of Esseen implies the proximity, in the Kolmogorov distance sense, of the distributions of the random variables of two random sequences satisfying a sort of general asymptotic relation. We then present several instances of this observation. A first example, partially motivated by the the statistical theory of high precision measurements, is given by a uniform asymptotic approximation to (g(X + μ _n )) _n∈?, where g is some smooth function, X is a random variable and (μ _n ) _n∈? is a sequence going to infinity; a multivariate version is also stated and proved. We finally present a second class of examples given by a randomization of the interesting parameter in some classical asymptotic formulas; namely, a generic Laplace's type integral, randomized by the sequence (μ _n X) _n∈?, X being a Gamma distributed random variable.     

A robust process capability index

The existing process capability indices (PCI's) assume that the distribution of the process being investigated is normal. For non-normal distributions, PCI's become unreliable in that PCI's may indicate the process is capable when in fact it is not. In this paper, we propose a new index which can be applied to any distribution. The proposed indexC_f:, is directly related to the probability of non-conformance of the process. For a given random sample, the estimation of C_f boils down to estimating non-parametrically the tail probabilities of an unknown distribution. The approach discussed in this paper is based on the works by Pickands (1975) and Smith (1987). We also discuss the construction of bootstrap confidence intervals of C_f: based on the so-called accelerated bias correction method (BC _a:). Several simulations are carried out to demonstrate the flexibility and applicability of C_f:. Two real life data sets are analyzed using the proposed index.     

Distribution of the product of a pair of independent two-sided power variates

ABSTRACT

The distributions of algebraic functions of random variables are important in theory of probability and statistics and other areas such as engineering, reliability, and actuarial applications, and many results based on various distributions are available in the literature. The two-sided power distribution is defined on a bounded range, and it is a generalization of the uniform, triangular, and power-function probability distributions. This paper gives the exact distribution of the product of two independent two-sided power-distributed random variables in a computable representation. The percentiles of the product are then computed, and a real data application is given.     

A characterization of probability distributions similar to the exponential

A concept of the lack-of-memory property at a given time point c > 0 is introduced. It is equivalent to the concept of the almost-lack-of-memory (ALM) property of the random variables. A representation theorem is given for the cumulative distribution function of such random variables as well as for corresponding decompositions in terms of independent random variables. It is shown that a periodic failure rate for a random variable is equivalent to the ALM property. In addition some properties of the service time of an unreliable server are observed.     

Marshall–Olkin Extended Lomax Distribution and Its Application to Censored Data

The problem of making statistical inference about θ =P(X > Y) has been under great investigation in the literature using simple random sampling (SRS) data. This problem arises naturally in the area of reliability for a system with strength X and stress Y. In this study, we will consider making statistical inference about θ using ranked set sampling (RSS) data. Several estimators are proposed to estimate θ using RSS. The properties of these estimators are investigated and compared with known estimators based on simple random sample (SRS) data. The proposed estimators based on RSS dominate those based on SRS. A motivated example using real data set is given to illustrate the computation of the newly suggested estimators.     

Inference of stress-strength for the Type-II generalized logistic distribution under progressively Type-II censored samples

This article studies the estimation of the reliability R = P[Y < X] when X and Y come from two independent generalized logistic distributions of Type-II with different parameters, based on progressively Type-II censored samples. When the common scale parameter is unknown, the maximum likelihood estimator and its asymptotic distribution are proposed. The asymptotic distribution is used to construct an asymptotic confidence interval of R. Bayes estimator of R and the corresponding credible interval using the Gibbs sampling technique have been proposed too. Assuming that the common scale parameter is known, the maximum likelihood estimator, uniformly minimum variance unbiased estimator, Bayes estimation, and confidence interval of R are extracted. Monte Carlo simulations are performed to compare the different proposed methods. Analysis of a real dataset is given for illustrative purposes. Finally, methods are extended for proportional hazard rate models.     

A Class of Multivariate Bilateral Selection t Distributions and Its Properties

This article proposes a class of multivariate bilateral selection t distributions useful for analyzing non-normal (skewed and/or bimodal) multivariate data. The class is associated with a bilateral selection mechanism, and it is obtained from a marginal distribution of the centrally truncated multivariate t. It is flexible enough to include the multivariate t and multivariate skew-t distributions and mathematically tractable enough to account for central truncation of a hidden t variable. The class, closed under linear transformation, marginal, and conditional operations, is studied from several aspects such as shape of the probability density function, conditioning of a distribution, scale mixtures of multivariate normal, and a probabilistic representation. The relationships among these aspects are given, and various properties of the class are also discussed. Necessary theories and two applications are provided.     

Weak disintegrability as a form of preservation of coherence

Summary Weak disintegrations are investigated from various points of view. Kolmogorov's definition of conditional probability is critically analysed, and it is noted how the notion of disintegrability plays some role in connecting Kolmogorov's definition with the one given in line with de Finetti's coherence principle. Conditions are given, on the domain of a prevision, implying the equivalence between weak disintegrability and conglomerability. Moreover, weak sintegrations are characterized in terms of coherence, in de Finetti's sense, of, a suitable function. This fact enables us to give, an interpretation of weak disintegrability as a form of “preservation of coherence”. The previous results are also applied to a hypothetical inferential problem. In particular, an inference is shown to be coherent, in the sense of Heath and Sudderth, if and only if a suitable function is coherent, in de Finetti's sense. Research partially supported by: M.U.R.S.T. 40% “Problemi di inferenza pura”.     

On the mean L1-error in the heteroscedastic deconvolution problem with compactly supported noises

We study the heteroscedastic deconvolution problem when random noises have compactly supported densities. In this context, the Fourier transforms of the densities can vanish on the real line. We propose a truncated type of estimator for target density and derive the convergence rate of the mean L¹-error uniformly over a class of target densities. A lower bound for the mean L¹-error is also established. Some simulations will be given to illustrate the performance of the proposed estimator.     

Bounds for moments through a general orthogonal expansion in a pre-hilbert space-II

In an earlier article [Canad. J. Statist., Vol, 3, No. 1, 1975, 13–34] bounds are obtained for the product moments of an arbitrary finite number of ordered random variables. These bounds are obtained with the help of a representation of an arbitrary function in terms of a complete orthonormal system in a pre-Hilbert space of square integrable functions defined in a k-dimensional unit cube. These results are extended to symmetric cases in this article.     

An RKHS framework for functional data analysis

Linear combinations of random variables play a crucial role in multivariate analysis. Two extension of this concept are considered for functional data and shown to coincide using the Loève–Parzen reproducing kernel Hilbert space representation of a stochastic process. This theory is then used to provide an extension of the multivariate concept of canonical correlation. A solution to the regression problem of best linear unbiased prediction is obtained from this abstract canonical correlation formulation. The classical identities of Lawley and Rao that lead to canonical factor analysis are also generalized to the functional data setting. Finally, the relationship between Fisher's linear discriminant analysis and canonical correlation analysis for random vectors is extended to include situations with function-valued random elements. This allows for classification using the canonical Y scores and related distance measures.     

Characterizations based on conditional expectations

For given real functionsg andh, first we give necessary and sufficient conditions such that there exists a random variableX satisfying thatE(g(X)|X≥y)=h(y)r _x (y),∀y ∈ C _x , whereC _x andT _X are the support and the failure rate function ofX, respectively. These extend the results of Ruiz and Navarro (1994) and Ghitany et al. (1995). Next we investigate necessary and sufficient conditions such thath(y)=E(g(X)|X≥y), for a given functionh. Support for this research was provided in part by the National Science Council of the Republic of China, Grant No. NSC 86-2115-M-110-014 and NSC 88-2118-M-110-001     

Bivariate Contours of Copula

Quantile functions associated with bivariate copulas are considered. Some of their structural properties are studied. Quantile functions allow one to express the cdf of the random variable C(X, Y), where (X, Y) is distributed as C(x, y) and where C is a copula. Quantile functions provide also a simple algorithm for simulating random observations.