Optimal risk rate in a functional inverse problem

In this note, we will consider the problem of recovering an unknown input function when the output function is observed in its entirety, blurred with functional error. An estimator is constructed whose risk converges at an optimal rate. In this functional model, convergence rates of order 1/n (n is the sample size) are possible, provided that the error distribution is sufficiently concentrated so as to compensate for the ill‐posedness of the inverse of the model operator.     

A note on generalized wald’s method

Let {v _n(θ)} be a sequence of statistics such that whenθ =θ ₀,v _n(θ ₀) N _p(0,Σ), whereΣ is of rankp andθ εR ^d. Suppose that underθ =θ ₀, {Σ _n} is a sequence of consistent estimators ofΣ. Wald (1943) shows thatv _n ^T (θ ₀)Σ _n ⁻¹ v _n(θ ₀) x ²(p). It often happens thatv _n(θ ₀) N _p(0,Σ) holds butΣ is singular. Moore (1977) states that under certain assumptionsv _n ^T (θ ₀)Σ _n ⁻ v _n(θ ₀) x ²(k), wherek = rank (Σ) andΣ _n ⁻ is a generalized inverse ofΣ _n. However, Moore’s result as stated is incorrect. It needs the additional assumption that rank (Σ _n) =k forn sufficiently large. In this article, we show that Moore’s result (as corrected) holds under somewhat different, but easier to verify, assumptions. Research partly supported by the U.S. Army Research Office through the Mathematical Sciences Institute at Cornell University.     

The problem of optimum stratification and allocation withq=n/N>0

Summary Dalenius/Gurney [1951] published necessary conditions for the stratum boundaries, so that with Neyman's optimal allocation of the sample sizen the variance of the sample mean will become a minimum. They introduced in the variance of the sample mean for the sample sizesn _h the opti mal values according to Neyman and differentiated this variance with respect to the stratum boundaries. Because Neyman's allocation formula yields only feasible solutions forn _h N _h , the conditions ofDalenius result in wrong, i.e. nonfeasible solutions, if one of the restrictionsn _h N _h (h=1 (1) L) is violated.By the example of a logarithmic normal distribution with =0, =1,5 forL=2 the behaviour of the Dalenius-Neyman-minimum and that of the feasible minimum will be shown in dependence on the sampling fractionq=n/N and a critical valueq _c will be given. For valuesq>q _c the Dalenius-Neyman-minimum is no longer feasible.For the same logarithmic normal distribution andL=2 (1) 10 this critical sampling fractionq _c will be given (section 5).For different values of andq the optimal stratum boundaries and sampling fractions are listed in section 6 forL=2;3;4.     

A method to obtain new copulas from a given one

Given a strictly increasing continuous function φ from [0, 1] to [0, 1] and its pseudo-inverse φ^[−1], conditions that φ must satisfy for C_φ(x₁, . . . ,x_n)=φ^[−1](C(φ(x₁), . . . ,φ(x_n))) to be a copula for any copula C are studied. Some basic properties of the copulas obtained in this way are analyzed and several examples of generator functions φ that can be used to construct copulas C_φ are presented. In this manner, a method to obtain from a given copula C a variety of new copulas is provided. This method generalizes that used to construct Archimedean copulas in which the original copula C is the product copula, and it is related with mixtures     

Het verdelen van steekproeven over subpopulaties bij accountantscontroles

Summary “Stratificationprocedures for a typical auditing problem”. During the past ten years, much experience was gained in The Netherlands in using random sampling methods for typical auditing problems. Especially, a method suggested by VAN. HEERDEN [2] turned out to be very fruitful. In this method a register of entries is considered to be a population of T guilders, if all entries total up to T guilders. The sample size n₀ is determined in such a way that the probability β not to find any mistake in the sample, if a fraction p₀ or more of T is incorrect, is smaller than a preassigned value β₀. So n₀ should satisfy (l-p)ⁿ⁰≤β₀ for p≥p₀. A complication arises if it is not possible to postpone sampling until the whole population T is available. One then wants to take samples from a population which is growing up to T. Suppose one is going to take samples n_i from e.g. r subpopulations Using the minimax procedure, it is shown, that in this case one should choose the sizes n_i equal to (T_i/T)n₀. The minimax-value of the probability not to find any incorrect guilder in the r samples, taken together is equal to β₀.     

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 .     

ON THE DISTRIBUTIONS OF ORDER STATISTICS FOR A RANDOM SAMPLE SIZE*

The probability distribution of the i –th and j–th order statistics and of the range R of a sample of size n, taken from a population with probability density function f (x) have been obtained when the sample size n is a random variable N and has: (i) a generalized Poisson distribution; and (ii) a generalized negative bonimial distribution. Specific results are then obtained; (a) when f (x) is uniform over (0,1); and (b) when f(x) is exponential. All the results for N, being a Poisson, binomial and negative binomial rv follow as special cases.     

A consistent test for multivariate normality based on the empirical characteristic function

LetX ₁,X ₂, …,X _n be independent identically distributed random vectors in IR ^d ,d ⩾ 1, with sample mean and sample covariance matrixS _n. We present a practicable and consistent test for the composite hypothesisH _d: the law ofX ₁ is a non-degenerate normal distribution, based on a weighted integral of the squared modulus of the difference between the empirical characteristic function of the residualsS _n ^−1/2 (X _j − ) and its pointwise limit exp (−1/2|t|²) underH _d. The limiting null distribution of the test statistic is obtained, and a table with critical values for various choices ofn andd based on extensive simulations is supplied.     

Slow convergence of the number of near-maxima for Burr XII distributions

A near-maximum is an observation which falls within a distance a of the maximum observation in an independent and identically distributed sample of size n. Subject to some conditions on the tail thickness of the population distribution, the number K _n (a) of near-maxima is known to converge in probability to one or infinity, or in distribution to a shifted geometric law. In this paper we show that for all Burr XII distributions K _n (a) converges almost surely to unity, but this convergence property may not become clear under certain cases even for very large n. We explore the reason of such slow convergence by studying a distributional continuity between Burr XII and Weibull distributions. We have also given a theoretical explanation of slow convergence of K _n (a) for the Burr XII distributions by showing that the rate of convergence in terms of P{K _n (a) > 1} tending to zero changes very little with the sample size n. Illustrations of the limiting behaviour K _n (a) for the Burr XII and the Weibull distributions are given by simulations and real data. The study also raises an important issue that although the Burr XII provides overall better fit to a given data set than the Weibull distribution, cautions should be taken for the extrapolation of the upper tail behaviour in the case of slow convergence.      

Learning curves—an axiomatic approach

A systematic mathematical analysis of learning curves is presented. It is shown that while a learning curve with learning factor k does necessarily satisfy the functional equation f(2x) = kf(x), this equation admits numerous other analytic, convex solutions as well so that it cannot be used to uniquely characterize learning curves. Rather, it is shown that a pair of functional equations f(rx) = kf(x) and f(sx) = jf(x) where rⁿ ≠ s^m for any positive integers n and m suffice to uniquely determine the learning curve; compatibility of the two equations requires that log_rk = log_sj or there will be no learning curve satisfying the pair of equations. Two classes of almost learning curves are generated and studied by means of suitable perturbation terms in the differential equation y = by/x of the true or standard learning curve and these curves are applied to describe data not exhibiting exact learning-curve behavior. Finally the concept of average marginal hours or cost is introduced and its behavior is found to also exhibit the learning-curve phenomenon except for an initial deviation.     

Asymptotic results and tests for the choice of approximative models in nonlinear two-phases regression models, heteroscedatic case

This paper is devoted to the study of the least squares estimator of f for the classical, fixed design, nonlinear model X (t _i)=f(t _i)+ε(t _i), i=1,2,…,n, where the (ε(t _i))_i=1,…,n are independent second order r.v.. The estimation of f is based upon a given parametric form. In Brodeau (1993) this subject has been studied in the homoscedastic case. This time we assume that the ε(t _i) have non constant and unknown variances σ²(t _i). Our main goal is to develop two statistical tests, one for testing that f belongs to a given class of functions possibly discontinuous in their first derivative, and another for comparing two such classes. The fundamental tool is an approximation of the elements of these classes by more regular functions, which leads to asymptotic properties of estimators based on the least squares estimator of the unknown parameters. We point out that Neubauer and Zwanzig (1995) have obtained interesting results for connected subjects by using the same technique of approximation. Received: February 1996     

Multivariate multistage methodologies for simultaneous all pairwise comparisons

We consider the problem of constructing simultaneous fixed-width confidence intervals for all pairwise treatment differences μ₁−μ _J , in the presence ofk(≥2) independent populationsN _p (μ₁,Σ), 1≤i≠j≤k. Appropriate purely sequential, accelerated sequential and three-stage sampling strategies have been developed and variousfirst-order asymptotic properties are then derived when Σ _pxp is completely unknown, but positive definite (p.d.). In the two special cases when the largest component variance in Σ is a known multiple of one of the variances or Σ=σ² H where σ(>0) is unknown, butH _pxp is known and p.d., the original multistage sampling strategies are specialized. Under such special circumstances, associatedsecond-order characteristics are then developed. It is to be noted that our present formulation and the methodologies fill important voids in the context of multivariate multiple comparisons which is a challenging area that has not yet been fully explored. Moderate sample performances of the proposed techniques were very encouraging and detailed remarks on these were included in Mukhopadhyay and Aoshima (1997).     

Testing a parametric quantile-regression model with an endogenous explanatory variable against a nonparametric alternative

This paper is concerned with inference about a function g that is identified by a conditional quantile restriction involving instrumental variables. The paper presents a test of the hypothesis that g belongs to a finite-dimensional parametric family against a nonparametric alternative. The test is not subject to the ill-posed inverse problem of nonparametric instrumental variable estimation. Under mild conditions, the test is consistent against any alternative model. In large samples, its power is arbitrarily close to 1 uniformly over a class of alternatives whose distance from the null hypothesis is proportional to n^−1/2, where n is the sample size. Monte Carlo simulations illustrate the finite-sample performance of the test.     

Minimax invariant estimation of a continuous distribution function under entropy loss

For the invariant decision problem of estimating a continuous distribution function F with two entropy loss functions, it is proved that the best invariant estimators d ₀ exist and are the same as the best invariant estimator of a continuous distribution function under the squared error loss function L (F, d)=∫|F (t) −d (t) |² dF (t). They are minimax for any sample size n≥1.     

Berry-Esséen rate in asymptotic normality for perturbed sample quantiles

In estimating quantiles with a sample of sizeN obtained from a distributionF, the perturbed sample quantiles based on a kernel functionk have been investigated by many authors. It is well known that their behaviour depends on the choices of “window-width”, sayw _N. Under suitable and reasonably mild assumptions onF andk, Ralescu and Sun (1993) have recently proven that lim _N→∞ N ^1/4w_N=0 is the necessary and sufficient condition for the asymptotic normality of the perturbed sample quantiles. In this paper, their rate of convergence is investigated. It turns out that the optimal Berry-Esséen rate ofO(N^?1/2) can be achieved by choosing the window-width suitably, sayw _N=O(N^?1/2). The obtained results, in addition to being explicit enough to verify the sufficient condition for the asymptotic normality, improve Ralescu's (1992) result of which the rate is of order (logN)N ^?1/2.     

Price asymptotics

This paper examines the pricing decisions of a seller facing an unknown demand function. It is assumed that partial information, in the form of an independent random sample of values, is available. The optimal price for the inferred demand satisfies a consistency property—as the size of the sample increases, the maximum profit and price approach the values for the case where demand is known. The main results deduced here are asymptotics for prices. Prices converge at a rate of O _p (n ^−1/3) with a limit that can be expressed as a functional of a Gaussian process. Implications for the comparison of mechanisms are discussed.      

The discrepancy distribution of stationary multiplier rules for rounding probabilities

The problem of rounding finitely many (continuous) probabilities to (discrete) proportions N _i/n is considered, for some fixed rounding accuracy n. It is well known that the rounded proportions need not sum to unity, and instead may leave a nonzero discrepancy D=(∑N _i) −n. We determine the distribution of D, assuming that the rounding function used is stationary and that the original probabilities follow a uniform distribution. Received: April 1999     

Estimating a restricted normal mean

Let X ₁, X ₂, ..., X _n be a random sample from a normal distribution with unknown mean μ and known variance σ ². In many practical situations, μ is known a priori to be restricted to a bounded interval, say [−m, m] for some m > 0. The sample mean , then, becomes an inadmissible estimator for μ. It is also not minimax with respect to the squared error loss function. Minimax and other estimators for this problem have been studied by Casella and Strawderman (Ann Stat 9:870–878, 1981), Bickel (Ann Stat 9:1301–1309, 1981) and Gatsonis et al. (Stat Prob Lett 6:21–30, 1987) etc. In this paper, we obtain some new estimators for μ. The case when the variance σ ² is unknown is also studied and various estimators for μ are proposed. Risk performance of all estimators is numerically compared for both the cases when σ ² may be known and unknown.     

TypeA operating characteristic functions of defects sampling plans

Investigations on acceptance sampling have attached rather few attention to defects inspection. As to sampling models and the corresponding operating characteristic (OC) function of single defects sampling plans, generally only typeB (process model) has been considered: sampling from a production process where the OC is conceived as a function of sample sizen, acceptance numberc, and process average number of defects per itemλ. For modern quality control, both the steadily increasing complexity of products and the need for differentiated cost calculation involve a clear demand for defects sampling in its practically most relevant form: lot-by-lot single sampling plans, where the OC (typeA, lot OC) is considered as a function of lot sizeN, sample sizen, acceptance numberc, number of defects in the lotK. The present paper develops two lot OC functions from suitable assumptions on the arrangements of the total numberK of defects over theN elements of the lot. Limiting theorems on these OC functions are used to justify to some extent the customary assumption that the Poisson process OC can serve as an approximation for typeA. The dependence of the OC functions on sample sizen, acceptance numberc, total number of defects in the lotK is described by simple formulae.