The average number of events per time unit and its estimation

Let ζ _t be the number of events which will be observed in the time interval [0;t] and define as the average number of events per time unit if this limit exists. In the case of i.i.d. waiting-times between the events,E[ζ _t ] is the renewal function and it follows from well-known results of renewal theory thatA exists and is equal to 1/τ, if τ>0 is the expectation of the waiting-times. This holds true also when τ = ∞.A may be estimate by ζ _t /t or where is the mean of the firstn waiting-timesX ₁,X ₂, ...,X _n . Both estimators converage with probability 1 to 1/τ if theX _i are i.i.d.; but the expectation of may be infinite for alln and also if it is finite, is in general a positively biased estimator ofA. For a stationary renewal process, ζ _t /t is unbiased for eacht; if theX _i are i.i.d. with densityf(x), then ζ _t /t has this property only iff(x) is of the exponential type and only for this type the numbers of events in consecutive time intervals [0,t], [t, 2t], ... are i.i.d. random variables for arbitraryt > 0.     

Exact bahadur efficiencies of tests for superadditivity of the mean function of a non-homogeneous poisson process

Consider a non-homogeneous Poisson process,N(t), with mean value functionΛ(t) and intensity functionsΛ(t). A conditional test of the hypothesis that the process is homogeneous, versus alternatives for whichΛ(t) is superadditive, was proposed by Hollander and Proschan (1974). Proposing a new test for superadditivity ofΛ(t), Kochar and Ramallingam (1989) have observed the fact that the Pitman asymptotic relative efficiency of their test with respect to the Hollander-Proschan test is unity. In order to distinguish between these competing tests, we shall compute the exact Bahadur slopes of these tests for important alternatives and demonstrate that the new test has high Bahadur efficiencies relative to the test of Hollander and Proschan.     

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     

Largest excess of boundary crossings for martingales

Summary LetX ₁,X ₂, ... form a sequence of martingale differences and denote byZ(a, α) = sup n (S _n −an ^α)⁺ the largest excess forS _n =X ₁ + ... +X _n crossing the boundaryan ^α. We give a sufficient condition for the finiteness ofEZ(a, α)^β which is formulated in terms of bounds forE(X _i ^{+ p} andE(|X _i |γ|X ₁, ...,X _i-1), whereα, β, γ, p are suitably related. This general result is then applied to the case of independent random variables.     

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(α₀).     

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     

An effective step-down algorithm for the construction and the identification of nonisomorphic orthogonal arrays

In this paper we present an effective algorithm for the construction and the identification of two-level nonisomorphic orthogonal arrays. Using this algorithm, we identify and list a full catalogue of nonisomorphic orthogonal arrays with parameters OA(24,7,2,t), OA(28,6,2,t) and OA(32,6,2,t), t ≥ 2. Some statistical properties of these designs are also considered.     

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).     

Moment consistency of estimators in partially linear models under NA samples

Consider the heteroscedastic regression model Y ^(j)(x _in , t _in ) = t _in β + g(x _in ) + σ _in e ^(j)(x _in ), 1 ≤ j ≤ m, 1 ≤ i ≤ n, where s_in²=f(u_in){\sigma_{in}^{2}=f(u_{in})}, (x _in , t _in , u _in ) are fixed design points, β is an unknown parameter, g(·) and f(·) are unknown functions, and the errors {e ^(j)(x _in )} are mean zero NA random variables. The moment consistency for least-squares estimators and weighted least-squares estimators of β is studied. In addition, the moment consistency for estimators of g(·) and f(·) is investigated.     

On Estimators of the Nearest Neighbour Distance Distribution Function for Stationary Point Processes

There are three approaches for the estimation of the distribution function D(r) of distance to the nearest neighbour of a stationary point process: the border method, the Hanisch method and the Kaplan-Meier approach. The corresponding estimators and some modifications are compared with respect to bias and mean squared error (mse). Simulations for Poisson, cluster and hard-core processes show that the classical border estimator has good properties; still better is the Hanisch estimator. Typically, mse depends on r, having small values for small and large r and a maximum in between. The mse is not reduced if the exact intensity λ (if known) or intensity estimators from larger windows are built in the estimators of D(r); in contrast, the intensity estimator should have the same precision as that of λ D(r). In the case of replicated estimation from more than one window the best way of pooling the subwindow estimates is averaging by weights which are proportional to squared point numbers.     

A bound on the expected overshoot for some concave boundaries

LetX ₁,X ₂,… be i.i.d. with finite meanμ>0,S _n =X ₁+…+X _n . Forf(n)=n ^β ,c>0 we consider the stopping timesT _c =inf{n:S _n >c+f(n)} with overshootR _c =S _{T c} −(c+f(T _c )). For 0<β<1 we give a bound for sup _c≥0 ER _c in the spirit of Lorden’s well-known inequality forf=0.     

Limit theorems for the spacings of weak records

Let W(1), W(2), . . . be weak record values obtained from a sample of independent variables with common discrete distribution. In the present paper, we derive weak and strong limit theorems for the spacings W(n + m) − W(n), m ≥ 1, n → ∞.     

Characterization of the multivariate normal distribution by conditional normal distributions

Summary LetX andY be two random vectors with values in ℝ ^k and ℝ∝, respectively. IfZ=(X ^T,Y ^T) ^T is multivariate normal thenX givenY=y andY givenX=x are (multivariate) normal; the converse is wrong. In this paper simple additional conditions are stated such that the converse is true, too. Furthermore, the case is treated that the random vectorZ=(X ₁ ^T , …,X _t ^T ) ^T is splitted intot≥3 partsX ₁, …,X _t.     

On the upper and lower semicontinuity of the Aumann integral

Let (T,τ,μ) be a finite measure space, X be a Banach space, P be a metric space and let L₁(μ,X) denote the space of equivalence classes of X-valued Bochner integrable functions on (T,τ,μ). We show that if φ:T×P→2^X is a set-valued function such that for each fixed pεP, φ(·,p) has a measurable graph and for each fixed tεT, φ(t,·) is either upper or lower semicontinuous then the Aumann integral of φ, i.e.,∫_Tφ(t,p)dμ(t)= {∫_Tx(t)dμ(t):xεS_φ(p)}, where S_φ(p)= {yεL₁(μ,X):y(t)εφ(t,p)μ−a.e.}, is either upper or lower semicontinuous in the variable p as well. Our results generalize those of Aumann (1965, 1976) who has considered the above problem for X=Rⁿ, and they have useful applications in general equilibrium and game theory.     

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 Bayesian-like estimator of the process capability index Cpmk

Pearn et al. (1992) proposed the capability index C_pmk, and investigated the statistical properties of its natural estimator for stable normal processes with constant mean μ. Chen and Hsu (1995) showed that under general conditions the asymptotic distribution of is normal if μ≠m, and is a linear combination of the normal and the folded-normal distributions if μ=m, where m is the mid-point between the upper and the lower specification limits. In this paper, we consider a new estimator for stable processes under a different (more realistic) condition on process mean, namely, P (μ≥m)=p, 0≤p≤1. We obtain the exact distribution, the expected value, and the variance of under normality assumption. We show that for P (μ≥m)=0, or 1, the new estimator is the MLE of C_pmk, which is asymptotically efficient. In addition, we show that under general conditions is consistent and is asymptotically unbiased. We also show that the asymptotic distribution of is a mixture of two normal distributions. RID="*" ID="*"  The research was partially supported by National Science Council of the Republic of China (NSC-89-2213-E-346-003).     

A minimal characterization of the covariance matrix

Summary LetX be ak-dimensional random vector with mean vectorμ and non-singular covariance matrix Σ. We show that among all pairs (a, Δ),a ∈ IR ^k , Δ ∈ IR ^k×k positive definite and symmetric andE(X−a)′ Δ⁻¹(X−a)=k, (μ, Σ) is the unique pair which minimizes det Δ. This motivates certain robust estimators of location and scale. Research supported by the Nuffield Foundation.     

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.     

Change-point analysis through integer-valued autoregressive process with application to some COVID-19 data

In this article, we consider the problem of change-point analysis for the count time series data through an integer-valued autoregressive process of order 1 (INAR(1)) with time-varying covariates. These types of features we observe in many real-life scenarios especially in the COVID-19 data sets, where the number of active cases over time starts falling and then again increases. In order to capture those features, we use Poisson INAR(1) process with a time-varying smoothing covariate. By using such model, we can model both the components in the active cases at time-point t namely, (i) number of nonrecovery cases from the previous time-point and (ii) number of new cases at time-point t. We study some theoretical properties of the proposed model along with forecasting. Some simulation studies are performed to study the effectiveness of the proposed method. Finally, we analyze two COVID-19 data sets and compare our proposed model with another PINAR(1) process which has time-varying covariate but no change-point, to demonstrate the overall performance of our proposed model.     

Multivariate likelihood ratio orderings between spacings of heterogeneous exponential random variables

Let X ₁, X ₂, ..., X _n be independent exponential random variables such that X _i has failure rate λ for i = 1, ..., p and X _j has failure rate λ^* for j = p + 1, ..., n, where p ≥ 1 and q = n − p ≥ 1. Denote by D _i:n (p,q) = X _i:n −X _i-1:n the ith spacing of the order statistics X _1:n ≤ X _2:n ≤ ... ≤ X _n:n , i = 1, ..., n, where X _0:n ≡ 0. The purpose of this paper is to investigate multivariate likelihood ratio orderings between spacings D _i:n (p,q), generalizing univariate comparison results in Wen et al.(J Multivariate Anal 98:743–756, 2007). We also point out that such multivariate likelihood ratio orderings do not hold for order statistics instead of spacings. Supported by National Natural Science Foundation of China, the Program for New Century Excellent Talents in University (No.: NCET-04-0569), and by the Knowledge Innovation Program of the Chinese Academy of Sciences (No.: KJCX3-SYW-S02).