On the modeling of size distributions when technologies are complex

The study considers a stochastic R&D process where the invented production technologies consist of a large number $n$ of complementary components. The degree of complementarity is captured by the elasticity of substitution of the CES aggregator function. Drawing from the Central Limit Theorem and the Extreme Value Theory we find, under very general assumptions, that the cross-sectional distributions of technological productivity are well-approximated either by the lognormal, Weibull, or a novel “CES/Normal” distribution, depending on the underlying elasticity of substitution between technology components. We find the tail of the “CES/Normal” distribution to be fatter than the Weibull tail but qualitatively thinner than the Pareto (power law) one. We also numerically assess the rate of convergence of the true technological productivity distribution to the theoretical limit with $n$ as fast in the body but slow in the tail.     

Multivariate gamma distributions with one-factorial accompanying correlation matrices and applications to the distribution of the multivariate range

Summary A new representation for the characteristic function of the joint distribution of the Mahalanobis distances betweenk independentN(μ, Σ)-distributed points is given. Especially fork=3 the corresponding distribution function is obtained as a special case of multivariate gamma distributions whose accompanying normal distribution has a positive semidefinite correlation matrix with correlationsϱ _ij=−a _i a _j. These gamma distribution functions are given here by one-dimensional parameter integrals. With some further trivariate gamma distributions third order Bonferroni inequalities are derived for the upper tails of the distribution function of the multivariate range ofk independentN(μ, I)-distributed points. From these inequalities very accurate (conservative) approximations to upperα-level bounds can also be computed for studentized multivariate ranges.     

Linearity of regression for non-adjacent order statistics

Let X ₁,X ₂,…,X _n be a random sample from a continuous distribution with the corresponding order statistics X _1:n≤X _2:n≤…≤X _n:n. All the distributions for which E(X _{k+r: n}|X _k:n)=a X _k:n+b are identified, which solves the problem stated in Ferguson (1967). Received February 1998     

One-Sample Bayesian Predictive Interval of Future Ordered Observations for the Pareto Distribution

Nigm et al. (2003, statistics 37: 527–536) proposed Bayesian method to obtain predictive interval of future ordered observation Y _(j) (r < j≤ n ) based on the right type II censored samples Y ₍₁₎ < Y ₍₂₎ < ... < Y _(r) from the Pareto distribution. If some of Y ₍₁₎ < ... < Y _(r-1) are missing or false due to artificial negligence of typist or recorder, then Nigm et al.’s method may not be an appropriate choice. Moreover, the conditional probability density function (p.d.f.) of the ordered observation Y _(j) (r < j ≤ n ) given Y ₍₁₎ <Y ₍₂₎ < ... < Y _(r) is equivalent to the conditional p.d.f. of Y _(j) (r < j ≤ n ) given Y _(r). Therefore, we propose another Bayesian method to obtain predictive interval of future ordered observations based on the only ordered observation Y _(r), then compares the length of the predictive intervals when using the method of Nigm et al. (2003, statistics 37: 527–536) and our proposed method. Numerical examples are provided to illustrate these results.     

Two families of kurtosis measures

Two families of kurtosis measures are defined as K ₁(b)=E[ab ^−|z|] and K ₂(b)=E[a(1−|z|^b)] where z denotes the standardized variable and a is a normalizing constant chosen such that the kurtosis is equal to 3 for normal distributions. K ₂(b) is an extension of Stavig's robust kurtosis. As with Pearson's measure of kurtosis β₂=E[z ⁴], both measures are expected values of continuous functions of z that are even, convex or linear and strictly monotonic in ℜ⁻ and in ℜ⁺. In contrast to β₂, our proposed kurtosis measures give more importance to the central part of the distribution instead of the tails. Tests of normality based on these new measures are more sensitive with respect to the peak of the distribution. K ₁(b) and K ₂(b) satisfy Van Zwet's ordering and correlate highly with other kurtosis measures such as L-kurtosis and quantile kurtosis. RID="*" ID="*"  The authors thank the referees for their insightful comments that significantly improved the clarity of the article.     

A note on reliability properties of k-record statistics

LetY _k,n denote the nth (upper) k-record value of an infinite sequence of independent, identically distributed random variables with common continuous distribution function F. We show that if the nth k-record valueY _k,n has an increasing failure rate (IFR), thenY _l,n (l<k) andY _k+1,n+1(n≤k+1) also have IFR distributions. On the other hand, ifY _k,n has a decreasing failure rate (DFR), thenY _l,n (1>k) has also a DFR distribution. We also present some results concerning log convexity and log concavity ofY _k,n .     

Two-sample rank tests with truncated populations

LetX ₁,…,X _m andY ₁,…,Y _n be two independent samples from continuous distributionsF andG respectively. Using a Hoeffding (1951) type theorem, we obtain the distributions of the vector S=(S ₍₁₎,…,S _(n)), whereS _(j)=# (X _i ’s≤Y _(j)) andY _(j) is thej-th order statistic ofY sample, under three truncation models: (a)G is a left truncation ofF orG is a right truncation ofF, (b)F is a right truncation ofH andG is a left truncation ofH, whereH is some continuous distribution function, (c)G is a two tail truncation ofF. Exploiting the relation between S and the vectorR of the ranks of the order statistics of theY-sample in the pooled sample, we can obtain exact distributions of many rank tests. We use these to compare powers of the Hajek test (Hajek 1967), the Sidak Vondracek test (1957) and the Mann-Whitney-Wilcoxon test. We derive some order relations between the values of the probagility-functions under each model. Hence find that the tests based onS ₍₁₎ andS _(n) are the UMP rank tests for the alternative (a). We also find LMP rank tests under the alternatives (b) and (c).     

Maximum likelihood methods for fitting the burr type XII distribution to multiply (progressively) censored life test data

Summary This paper develops mathematical and computational methodology for fitting, by the method of maximum likelihood (ML), the Burr Type XII distribution to multiply (or progressively) censored life test data. Mathematical expressions are given for approximating the asymptotic variances and covariances of the ML estimates (MLEs) of the parameters of the Burr Type XII distribution. A rigorous mathematical analysis is undertaken to investigate the existence and uniqueness of the MLEs for arbitrary sample data. The methodology of this paper is applied to progressively censored sample data arising in a life test experiment.     

Asymptotic normality of M-estimators in a semiparametric model with longitudinal data

Suppose that the longitudinal observations (Y _ij , X _ij , t _ij ) for i = 1, . . . ,n; j = 1, . . . ,m _i are modeled by the semiparamtric model where β ₀ is a k × 1 vector of unknown parameters, g(·) is an unknown estimated function and e _ij are unobserved disturbances. This article consider M-type regressions which include mean, median and quantile regressions. The M-estimator of the slope parameter β ₀ is obtained through piecewise local polynomial approximation of the nonparametric component. The local M-estimator of g(·) is also obtained by replacing β ₀ in model with its M-estimator and using local linear approximation. The asymptotic distribution of the estimator of β ₀ is derived. The asymptotic distributions of the local M-estimators of g(·) at both interior and boundary points are also established. Various applications of our main results are given. The research is supported in part by National Natural Science Foundation of China (Grant No. 10671089).     

Prediction of <Emphasis Type="Italic">k</Emphasis>-records from a general class of distributions under balanced type loss functions

We study the problem of predicting future k-records based on k-record data for a large class of distributions, which includes several well-known distributions such as: Exponential, Weibull (one parameter), Pareto, Burr type XII, among others. With both Bayesian and non-Bayesian approaches being investigated here, we pay more attention to Bayesian predictors under balanced type loss functions as introduced by Jafari Jozani et al. (Stat Probab Lett 76:773–780, 2006a). The results are presented under the balanced versions of some well-known loss functions, namely squared error loss, Varian’s linear-exponential loss and absolute error loss or L ₁ loss functions. Some of the previous results in the literatures such as Ahmadi et al. (Commun Stat Theory Methods 34:795–805, 2005), and Raqab et al. (Statistics 41:105–108, 2007) can be achieved as special cases of our results. Partial support from Ordered and Spatial Data Center of Excellence of Ferdowsi University of Mashhad is acknowledged by J. Ahmadi. M. J. Jozani’s research supported partially by a grant of Statistical Research and Training Center. é. Marchand’s research supported by NSERC of Canada. A. Parsian’s research supported by a grant of the Research Council of the University of Tehran.     

Characterization of some types of completeness resp. Total completeness and their conservation under direct products

Some notions ofL _p (μ)-completeness resp. totally L _p (μ)-completeness (1≦p≦∞) are characterized for families of probability distributions dominated by aσ-finite measureμ and their conservation with respect to direct products is proved. Furthermore, it is shown that totallyL _∞(μ)-completeness does not implyL ₁(μ)-completeness and that there are families of probability distributions in the i.i.d. case induced by the order statistic, which are L₁(μ)-complete but not totallyL _∞(μ)-complete.     

Minimax estimation of a cumulative distribution function by converting to a parametric problem

Let X = (X ₁,...,X _n ) be a sample from an unknown cumulative distribution function F defined on the real line . The problem of estimating the cumulative distribution function F is considered using a decision theoretic approach. No assumptions are imposed on the unknown function F. A general method of finding a minimax estimator d(t;X) of F under the loss function of a general form is presented. The method of solution is based on converting the nonparametric problem of searching for minimax estimators of a distribution function to the parametric problem of searching for minimax estimators of the probability of success for a binomial distribution. The solution uses also the completeness property of the class of monotone decision procedures in a monotone decision problem. Some special cases of the underlying problem are considered in the situation when the loss function in the nonparametric problem is defined by a weighted squared, LINEX or a weighted absolute error.     

On characterizing discrete distributions via conditional expectations of weak record values

Let (W _n ,n ≥ 0) denote the sequence of weak records from a distribution with support S = { α₀,α₁,...,α _N }. In this paper, we consider regression functions of the form ψ _n (x) = E(h(W _n ) |W _n+1 = x), where h(·) is some strictly increasing function. We show that a single function ψ _n (·) determines F uniquely up to F(α₀). Then we derive an inversion formula which enables us to obtain F from knowledge of ψ _n (·), ψ _n-1(·), h(·) and F(α₀).     

An identity for the product moments of order statistics

A general identity for the product moments of successive order statistics is given, which is valid in a class of probability distributions including Weibull, Pareto, exponential and Burr distributions.     

Convergence of a branching type recursion with non-stationary immigration

The asymptotic distribution of a branching type recursion with non-stationary immigration is investigated. The recursion is given by , where (X _l ) is a random sequence, (L _n ⁻¹⁽¹⁾ ) are iid copies ofL _n−1,K is a random number andK, (L _n ⁻¹⁽¹⁾ ), {(X _l ),Y _n } are independent. This recursion has been studied intensively in the literature in the case thatX _l ≥0,K is nonrandom andY _n =0 ∀n. Cramer, Rüschendorf (1996b) treat the above recursion without immigration with starting conditionL ₀=1, and easy to handle cases of the recursion with stationary immigration (i.e. the distribution ofY _n does not depend on the timen). In this paper a general limit theorem will be deduced under natural conditions including square-integrability of the involved random variables. The treatment of the recursion is based on a contraction method. The conditions of the limit theorem are built upon the knowledge of the first two moments ofL _n . In case of stationary immigration a detailed analysis of the first two moments ofL _n leads one to consider 15 different cases. These cases are illustrated graphically and provide a straight forward means to check the conditions and to determine the operator whose unique fixed point is the limit distribution of the normalizedL _n .     

Minimax estimation under convex loss when the parameter interval is bounded

In the present paper families of truncated distributions with a Lebesgue density forx=(x ₁,...,x _n ) ε ℝ ⁿ are considered, wheref ₀:ℝ → (0, ∞) is a known continuous function andC _n (ϑ) denotes a normalization constant. The unknown truncation parameterϑ which is assumed to belong to a bounded parameter intervalΘ=[0,d] is to be estimated under a convex loss function. It is studied whether a two point prior and a corresponding Bayes estimator form a saddle point when the parameter interval is sufficiently small.     

Goodness-of-fit tests for discrete models based on the integrated distribution function

This paper presents a new widely applicable omnibus test for discrete distributions which is based on the difference between the integrated distribution function Ψ(t)=∫_t ^∞ (1−F(x))dx and its empirical counterpart. A bootstrap version of the test for common lattice models has accurate error rates even for small samples and exhibits high power with respect to competitive procedures over a large range of alternatives. Received: July 1998