A Bayesian hierarchical spatio-temporal rainfall model

A Bayesian hierarchical spatio-temporal rainfall model is presented and analysed. The model has the ability to deal with extensive missing or null values, uses a sophisticated variance stabilising rainfall pre-transformation, incorporates a new elevation model and can provide sub-catchment rainfall estimation and interpolation using a sequential kriging scheme. The model uses a vector autoregressive stochastic process to represent the time dependence of the rainfall field and an exponential covariogram to model the spatial correlation of the rainfall field. The model can be readily generalised to other types of stochastic processes. In this paper, some results of applying the model to a particular rainfall catchment are presented.     

Niveaudurchagangszeiten zur charakterisierung sequentieller schätzverfahren

Let θ be a parameter of a homogenous additive stochastic process. In order to get an unbiased and efficient estimator for a function h(v) one has often to use sequential procedures. In this paper we consider processes of the socalled exponential class. We study level crossing times, which characterize certain sequential estimations. It is shown that the family of level crossing times for an increasing sequence of levels is also a process of the exponential class. The density function of the one-dimensional probability distributions of this new process is given Examples and applications conclude the paper.     

Sequential Designs for Binary Data with the Purpose to Maximize the Probability of Response

Two kinds of sequential designs are proposed for finding the point that maximizes the probability of response assuming a binary response variable and a quadratic logistic regression model. One is a parametric optimal design approach, and the other one is a nonparametric stochastic approximation approach. The suggested sequential designs are evaluated and compared in a simulation study. In summary, the parametric approach performed very well whereas its competitor failed in some cases.     

A UNIFIED FRAMEWORK FOR MODELLING WILDLIFE POPULATION DYNAMICS†

This paper proposes a unified framework for defining and fitting stochastic, discrete‐time, discrete‐stage population dynamics models. The biological system is described by a state‐space model, where the true but unknown state of the population is modelled by a state process, and this is linked to survey data by an observation process. All sources of uncertainty in the inputs, including uncertainty about model specification, are readily incorporated. The paper shows how the state process can be represented as a generalization of the standard Leslie or Lefkovitch matrix. By dividing the state process into subprocesses, complex models can be constructed from manageable building blocks. The paper illustrates the approach with a model of the British grey seal metapopulation, using sequential importance sampling with kernel smoothing to fit the model.     

Bayesian Dynamic Factor Models and Portfolio Allocation

We discuss the development of dynamic factor models for multivariate financial time series, and the incorporation of stochastic volatility components for latent factor processes. Bayesian inference and computation is developed and explored in a study of the dynamic factor structure of daily spot exchange rates for a selection of international currencies. The models are direct generalizations of univariate stochastic volatility models and represent specific varieties of models recently discussed in the growing multivariate stochastic volatility literature. We discuss model fitting based on retrospective data and sequential analysis for forward filtering and short-term forecasting. Analyses are compared with results from the much simpler method of dynamic variance-matrix discounting that, for over a decade, has been a standard approach in applied financial econometrics. We study these models in analysis, forecasting, and sequential portfolio allocation for a selected set of international exchange-rate-return time series. Our goals are to understand a range of modeling questions arising in using these factor models and to explore empirical performance in portfolio construction relative to discount approaches. We report on our experiences and conclude with comments about the practical utility of structured factor models and on future potential model extensions.     

Sequential estimation for a family of counting processes in the nuisance parameter case

The paper deals with the problem of sequential estimation for stochastic processes in the presence of a nuisance parameter. Using the approach to estimation through estimating equations, optimum estimating functions based on a random observation time are investigated in some models for processes appearing in reliability systems theory.     

Risk efficient estimation of fully dependent random coefficient autoregressive models of general order

We consider a stochastic dynamic model with autoregressive progression. The drift coefficients of the autoregressive model are random where the randomness in the coefficients can have any dependence structure. We propose a two-step sequential estimator and study the asymptotic behavior of few important properties. Paradigm of sequential estimation has its own advantage in reducing sample size and plugging estimates of nuisance parameters while inferring about the main parameters. Our proposed estimator is asymptotically optimal as the predictive risk of the proposed estimator attains the risk of the oracle that assumes known nuisance parameters. Extensive simulation confirms our results.     

BERNOULLI SAMPLING

The concepts of the Bernoulli count process of a point process and Bernoulli sampling of a discrete parameter stochastic process are introduced. The Bernoulli count process determines the stochastic structure of the point process, and a process obtained by thinning a discrete parameter stochastic process by Bernoulli sampling satisfies the same property. Stationarity and the Markov property remain invariant under Bernoulli sampling.     

2005年汇改后人民币汇率多次结构突变检验分析——基于循序检验方法

采用最新的多次结构突变循序检验方法,对2005年7月21日汇改后人民币汇率时间序列趋势项是否具有多次结构突变进行研究,并在多次结构突变检验结果的基础上对消除趋势后的人民币汇率数据进行分析,结果发现：人民币汇率时间序列是围绕着4个结构断点的分段趋势平稳的;人民币汇率服从分段趋势平稳的结论对汇率政策有效性、汇率与其他经济总量关系研究及汇率预测具有重要意义。     

Group sequential tests under fractional Brownian motion in monitoring clinical trials

Group sequential tests have been effective tools in monitoring long term clinical trials. There have been several popular discrete sequential boundaries proposed for modeling interim analysis of clinical trials under the assumption of Brownian motion for the stochastic processes generated from test statistics. In this paper, we study the five sequential boundaries in Lan and DeMets (Biometrika 70:659–663, 1983) under the fractional Brownian motion. The fractional Brownian includes the classic Brownian motion as a special case. An example from a real data set is used to illustrate the applications of the boundaries.     

Sequential testing of gradual changes in the drift of a stochastic process

In this paper, some sequential monitoring procedures are constructed and analyzed for detecting a “gradual” change in the drift parameter of a general stochastic process satisfying a certain (weak) invariance principle. It is shown that the tests can be constructed such that the “false alarm rate” attains a prescribed level (say) α and that the tests have “asymptotic power 1”. A more precise analysis of the procedures under the alternative proves that the stopping times, suitably normalized, have a standard normal limiting distribution. A few results from a small simulation study are also presented in order to give an idea of the finite sample behaviour of the suggested procedures.     

Change-points in stochastic ordering

This article studies the problem of testing and locating changepoints in stochas¬tic ordering. We propose a sequential process to detect the changepoints from two multinomial distributions. We also obtain the maximum likelihood estimators of two multinomial probability vectors under the assumption that the cumulative distribu¬tions have a changepoint. Asymptotically unbiased Akaike's information criterion is used to estimate the changepoints of two discrete probability distributions. Finally. we demonstrate our procedure by studying a data set pertaining to average daily insulin dose from the Boston Collaborative Drug Surveillance Program and locate the changepoints in stochastic ordering.     

Interim analysis of clinical trials based on urn models

Clinical trials usually involve efficient and ethical objectives such as maximizing the power and minimizing the total failure number. Interim analysis is now a standard technique in practice to achieve these objectives. Randomized urn models have been extensively studied in the literature. In this paper, we propose to perform interim analysis on clinical trials based on urn models and study its properties. We show that the urn composition, allocation of patients and parameter estimators can be approximated by a joint Gaussian process. Consequently, sequential test statistics of the proposed procedure converge to a Brownian motion in distribution and the sequential test statistics asymptotically satisfy the canonical joint distribution defined in Jennison & Turnbull (Jennison & Turnbull 2000. Group Sequential Methods with Applications to Clinical Trials, Chapman and Hall/CRC). These results provide a solid foundation and open a door to perform the interim analysis on randomized clinical trials with urn models in practice. Furthermore, we demonstrate our proposal through examples and simulations by applying sequential monitoring and stochastic curtailment techniques. The Canadian Journal of Statistics 40: 550–568; 2012 © 2012 Statistical Society of Canada     

On Comparison of Stopping Times in Sequential Procedures for Exponential Families of Stochastic Processes

For curved ( k + 1), k -exponential families of stochastic processes a natural and often studied sequential procedure is to stop observation when a linear combination of the coordinates of the canonical process crosses a prescribed level. For such procedures the model is, approximately or exactly, a non-curved exponential family. Subfamilies of these stopping rules defined by having the same Fisher (expected) information are considered. Within a subfamily the Bartlett correction for a point hypothesis is also constant. Methods for comparing the durations of the sampling periods for the stopping rules in such a subfamily are discussed. It turns out that some stopping times tend to be smaller than others. For exponential families of diffusions and of counting processes the probability that one such stopping time is smaller than another can be given explicity. More generally, an Edgeworth expansion of this probability is given     

Invariance principles for jackknifing U-statistics for finite population sampling and some applications

For simple random sampling (without replacement) from a finite population, suitable stochastic processes are constructed from the entire sequence of jackknife estimators based on smooth functions of U-statistics and these are approximated (in distributions) by some Brownian bridge processes. Strong convergence of the Tukey estimator of the variance of a jackknife U-statistic has been interpreted suitably and established. Some applications of these results in sequential analysis relating to finite population sampling are also considered.     

Testing for a Change of the Innovation Distribution in Nonparametric Autoregression: The Sequential Empirical Process Approach

We consider a nonparametric autoregression model under conditional heteroscedasticity with the aim to test whether the innovation distribution changes in time. To this end, we develop an asymptotic expansion for the sequential empirical process of nonparametrically estimated innovations (residuals). We suggest a Kolmogorov–Smirnov statistic based on the difference of the estimated innovation distributions built from the first ?ns?and the last n ? ?ns? residuals, respectively (0 ≤ s ≤ 1). Weak convergence of the underlying stochastic process to a Gaussian process is proved under the null hypothesis of no change point. The result implies that the test is asymptotically distribution‐free. Consistency against fixed alternatives is shown. The small sample performance of the proposed test is investigated in a simulation study and the test is applied to a data example.     

On drift parameter estimation in models with fractional Brownian motion

We consider a stochastic differential equation involving standard and fractional Brownian motion with unknown drift parameter to be estimated. We investigate the standard maximum likelihood estimate of the drift parameter, two non-standard estimates and three estimates for the sequential estimation. Model strong consistency and some other properties are proved. The linear model and Ornstein–Uhlenbeck model are studied in detail. As an auxiliary result, an asymptotic behaviour of the fractional derivative of the fractional Brownian motion is established.     

Optimal stopping for Brownian motion with applications to sequential analysis and option pricing

Herman Chernoff made fundamental contributions to analytical and computational methods for solving optimal stopping problems for Brownian motion. He also showed how these optimal stopping problems are closely related to some basic problems in sequential analysis and singular stochastic control. This paper gives a survey of these and related developments and describes some recent applications to option valuation in financial economics.     

Repeated Confidence Intervals Under Fractional Brownian Motion in Long-Term Clinical Trials

Repeated confidence interval (RCI) is an important tool for design and monitoring of group sequential trials according to which we do not need to stop the trial with planned statistical stopping rules. In this article, we derive RCIs when data from each stage of the trial are not independent thus it is no longer a Brownian motion (BM) process. Under this assumption, a larger class of stochastic processes fractional Brownian motion (FBM) is considered. Comparisons of RCI width and sample size requirement are made to those under Brownian motion for different analysis times, Type I error rates and number of interim analysis. Power family spending functions including Pocock, O'Brien-Fleming design types are considered for these simulations. Interim data from BHAT and oncology trials is used to illustrate how to derive RCIs under FBM for efficacy and futility monitoring.