Classes of diffusion-type processes with a sufficient reduction

Classes of processes of the diffusion type permitting a sufficient data reduction are derived. None of these classes are exponential families in the usual sense. For one type of such classes the sufficient statistic equals that of a curved exponential family of diffusion-type processes. For a second type the last observation is sufficient. In particular cases both types of classes are defined by means of a RICCATI equation     

On Minimum Distance Estimation in Recurrent Markov Step Processes. I

Consider a Markov step process X=(X_t)_t≥0 whose generator depends on an unknown d -dimensional parameter ϑ. We look at certain empirical measures for recurrent Markov step processes and their a.s. convergence; based on this, we introduce a class of minimum distance estimators. For broad families of sequential observation schemes (at stage n, the trajectory of X is observed up to time S_n, (S_n)_n a sequence of stopping times increasing to ∞), we formulate a stochastic expansion of the suitably rescaled estimation error; for a particular scheme, asymptotic normality is obtained as n →∞. A minimax property under misspecification of the model (in the sense that the true probability law is contiguous to the parametric model but not contained in it) is given.     

Multivariate geometric distributions

Families of multivariate geometric distributions with flexible correlations can be constructed by applying inverse sampling to a sequence of multinomial trials, and counting outcomes in possibly overlapping categories. Further multivariate families can be obtained by considering other stopping rules, with the possibility of different stopping roles for different counts, A simple characterisation is given for stopping rules which produce joint distributions with marginals having the same form as that of the number of trials. The inverse sampling approach provides a unified treatment of diverse results presented by earlier authors, including Goldberg (1934), Bates and Meyman (1952), Edwards and Gurland (1961), Hawkes (1972), Paulson and Uppulori (1972) and Griffiths and Milne (1987). It also provides a basis for investigating the range of possible correlations for a given set of marginal parameters. In the case of more than two joint geometric or negative binomial variables, a convenient matrix formulation is provided.     

A plant-capture approach for achieving complete coverage of a population

Sightings of, any member of a target population of unknown size occur according to a homogeneous Poisson process. The Poisson processes are independent and have a common rate. Individuals are tagged when first sighted. The population is augmented by the insertion of known number of planted individuals who have already been tagged and whose behaviour is identical to that of the other members. Consideration is given to two stopping rules for determining, within a specified error probability, when all members of the target population have been seen. A comparison is made of the average time to termination under these two rules. It is shown that, when the target population is small, the use of plants can provide a useful reduction in the average time taken to achieve complete coverage within the specified error probability.     

Selecting the optimum k in ridge regression

Two ridge rules are proposed for selecting the optimal k in ridge regression . Since the sampling distribution of the proposed rules are mathematically in tractable , a Monte Carlo study is conducted to examine their statisticl properties . Numerical results of the simulations in dicate that the performance of ridge rules depends upon the risk function used. Nevertheless, one of the ridge rules does produce a smaller mean squared error than the least squares estimator with the probability greater than 0.57 for all situations.     

Estimation method of multivariate exponential probabilities based on a simplex coordinates transform

A novel unbiased estimator for estimating the probability mass of a multivariate exponential distribution over a measurable set is introduced and is called the exponential simplex (ES) estimator. For any measurable set and given sample size, the statistical efficiency of the ES estimator is higher than or equal to the statistical efficiency of the well-known Monte Carlo (MC) estimator. For non-radially shaped measurable sets, the ES estimator has a strictly higher statistical efficiency than the MC estimator. For ray-convex sets, such as convex sets, the ES estimator can be expressed in a simple analytical form.     

Some recent results on the distributions of stopping times of compound Poisson processes with linear boundaries

A review of recent results is given, for the distribution of stopping times defined on compound Poisson processes and linear boundaries. Generalization of the results of Zacks (Commun. Statist. Stochastic Models 7 (1991) 233) is given for discrete compound Poisson processes. The main approach in the reviewed papers is that of sample path analysis, which leads to renewal type equations. Explicit solutions are given.     

On a sequential subset selection procedure for the least probable multinomial cell

This paper studies a sequential procedure R for selecting a random size subset that contains the multinomial cell which has the smallest cell probability. The stopping rule of the proposed procedure R is the composite of the stopping rules of curtailed sampling, inverse sampling, and the Ramey-Alam sampling. A reslut on the worst configuration is shown and it is employed in computing the procedure parameters that guarantee certain probability requirements. Tables of these procedure parameters, the corresponding probability of correct selection, the expected sample size, and the expected subset size are given for comparison purpose.     

Empirical Bayes rules and Gaussian processes

In this paper we study empirical Bayes (e.B.) rules from a viewpoint which has not yet got any attention in the literature. Since an e.B. estimator can be seen as an estimate of an unknown function, namely the true Bayes estimator, it is natural to consider e.B. estimators as stochastic processes. In this paper we make a first attempt in the direction of this approach. For a certain class of e.B. estimators for the continuous one-parameter exponential family, we investigate the global behaviour on finite intervals. It is shown that the difference between the e.B. and the true Bayes estimator can be represented as a certain type of Gaussian process plus a remainder which is uniformly of smaller order. Several applications of this result are given.     

Verallgemeinerte bayes-verfahren bei sequentiellen testproblemen

In this paper, problems of sequential decision theory are taken into consideration by extending the definition of the BAYES rule and treating BAYES rules. This generalisation is quite useful for practice. In many cases only BAYES rules can be calculated. The conditions under which such sequential decision procedures exist are demonstrated, as well as how to construct them on a scheme of backward induction resulting in the conclusion that the existence of BAYES rules needs essentially weaker assumptions than the existence of BAYES rules.Futhermore, methods are searched to simplify the construction of optimal stopping rules. Some illustrative examples are given.     

First-passage densities of controlled Gaussian processes

Whittle has proved a theorem that gives the optimal control of Gaussian processes in terms of the mathematical expectation of a function of the time and the place where the uncontrolled processes hit the boundary of the stopping region for the first time. In this paper we obtain formulae for the joint probability density function of the first hitting time and place and, in the time-invariant case, for the moment generating function of the first exit time of the optimally controlled processes. Two particular one-dimensional cases are considered.     

Inversion of the Renewal Density with Dead-time

Stationary renewal point processes are defined by the probability distribution of the distances between successive points (lifetimes) that are independent and identically distributed random variables. For some applications it is also interesting to define the properties of a renewal process by using the renewal density. There are well-known expressions of this density in terms of the probability density of the lifetimes. It is more difficult to solve the inverse problem consisting in the determination of the density of the lifetimes in terms of the renewal density. Theoretical expressions between their Laplace transforms are available but the inversion of these transforms is often very difficult to obtain in closed form. We show that this is possible for renewal processes presenting a dead-time property characterized by the fact that the renewal density is zero in an interval including the origin. We present the principle of a recursive method allowing the solution of this problem and we apply this method to the case of some processes with input dead-time. Computer simulations on Poisson and Erlang (2) processes show quite good agreement between theoretical calculations and experimental measurements on simulated data.     

Pivotal Quantities Based on Sequential Data: A Bootstrap Approach

A bootstrap algorithm is provided for obtaining a confidence interval for the mean of a probability distribution when sequential data are considered. For this kind of data the empirical distribution can be biased but its bias is bounded by the coefficient of variation of the stopping rule associated with the sequential procedure. When using this distribution for resampling the validity of the bootstrap approach is established by means of a series expansion of the corresponding pivotal quantity. A simulation study is carried out using Wang and Tsiatis type tests and considering the normal and exponential distributions to generate the data. This study confirms that for moderate coefficients of variation of the stopping rule, the bootstrap method allows adequate confidence intervals for the parameters to be obtained, whichever is the distribution of data.     

The 'demon' problem of Youden: Exponential case

We consider the problem of finding the probability of a sample mean falling above the (n - k)th-order statistic in a random sample of size n. Explicit expressions are obtained for the exponential distribution. Some applications that pertain to testing for outliers and goodness of fit are given.     

Ordering of Sequentially Sampled Exponential Experiments

Let X ₁, X ₂, ... be a sequence of i.i.d. random variables, X _i∼ F _θ, θ∈Θ. Let N ₁ and N ₂ be two stopping rules. For a class of exponential families { F _θ: θ∈Θ} we show that the experiment Y ₁ = ( X ₁, ..., X _N1) carries more statistical information than Y ₂ = ( X ₁, ..., x _N2) only if N ₁ is stochastically larger then N ₂     

Progressively censored tests for clinical experiments and life testing problems based on weighted empirical distributions

In the context of time-sequential studies, progressively censored tests for a simple regression model based on weighted empirical distributions are considered for ungrouped as well as grouped data situations. Early decision rules based on such tests are formulated. The asymptotic theory of the proposed tests rests on a construction of suitable empirical processes and their convergence (in distribution) to appropriate Gaussian functions. Critical values of the proposed test statistics are obtained by simulation, For a hypothetical example (of practical interest), a comparative study is made for the empirical powers and stopping times for some rival tests.     

Interacting neighbour point processes: Some models for clustering

We introduce a class of spatial point processesinteracting neighbour point (INP) processes, where the density of the process can be written by means of local interactions between a point and subsets of its neighbourhood but where the processes may not be Ripley-Kelly Markov processes with respect to this neighbourhood. We show that the processes are iterated Markov processes defined by Hayat and Gubner (1996). Furthermore, we pay special attention to a subclass of interacting neighbour processes, where the density belongs to the exponential family and all neighbours of a point affect it simultaneously. A simulation study is presented to show that some simple processes of this subclass can produce clustered patterns of great variety. Finally, an empirical example is given.     

Identifiability and estimation of recursive max‐linear models

We address the identifiability and estimation of recursive max‐linear structural equation models represented by an edge‐weighted directed acyclic graph (DAG). Such models are generally unidentifiable and we identify the whole class of DAG s and edge weights corresponding to a given observational distribution. For estimation, standard likelihood theory cannot be applied because the corresponding families of distributions are not dominated. Given the underlying DAG, we present an estimator for the class of edge weights and show that it can be considered a generalized maximum likelihood estimator. In addition, we develop a simple method for identifying the structure of the DAG. With probability tending to one at an exponential rate with the number of observations, this method correctly identifies the class of DAGs and, similarly, exactly identifies the possible edge weights.     

Crossing probabilities for a square root boundary by a bessel process

Probabilities for Wiener and Bessel processes to cross a square root boundary are calculated and the Mellin transform of the distribution of the associated stopping time is given.The transform is inverted and a table of critical values for the crossing probability as a function of the observation time is included. Also an analytic expression for the crossing probability is computed for the case of constant drift.     

Uniform Convergence of Compactly Supported Wavelet Expansions of Gaussian Random Processes

New results on uniform convergence in probability for expansions of Gaussian random processes using compactly supported wavelets are given. The main result is valid for general classes of non stationary processes. An application of the obtained results to stationary processes is also presented. It is shown that the convergence rate of the expansions is exponential.