The bootstrap and Bayesian bootstrap method in assessing bioequivalence

Parametric method for assessing individual bioequivalence (IBE) may concentrate on the hypothesis that the PK responses are normal. Nonparametric method for evaluating IBE would be bootstrap method. In 2001, the United States Food and Drug Administration (FDA) proposed a draft guidance. The purpose of this article is to evaluate the IBE between test drug and reference drug by bootstrap and Bayesian bootstrap method. We study the power of bootstrap test procedures and the parametric test procedures in FDA (2001). We find that the Bayesian bootstrap method is the most excellent.     

The semi-Markov beta-Stacy process: a Bayesian non-parametric prior for semi-Markov processes.

The literature on Bayesian methods for the analysis of discrete-time semi-Markov processes is sparse. In this paper, we introduce the semi-Markov beta-Stacy process, a stochastic process useful for the Bayesian non-parametric analysis of semi-Markov processes. The semi-Markov beta-Stacy process is conjugate with respect to data generated by a semi-Markov process, a property which makes it easy to obtain probabilistic forecasts. Its predictive distributions are characterized by a reinforced random walk on a system of urns.

     

Bayesian Monte Carlo estimation for profile hidden Markov models

Hidden Markov models are used as tools for pattern recognition in a number of areas, ranging from speech processing to biological sequence analysis. Profile hidden Markov models represent a class of so-called “left–right” models that have an architecture that is specifically relevant to classification of proteins into structural families based on their amino acid sequences. Standard learning methods for such models employ a variety of heuristics applied to the expectation-maximization implementation of the maximum likelihood estimation procedure in order to find the global maximum of the likelihood function. Here, we compare maximum likelihood estimation to fully Bayesian estimation of parameters for profile hidden Markov models with a small number of parameters. We find that, relative to maximum likelihood methods, Bayesian methods assign higher scores to data sequences that are distantly related to the pattern consensus, show better performance in classifying these sequences correctly, and continue to perform robustly with regard to misspecification of the number of model parameters. Though our study is limited in scope, we expect our results to remain relevant for models with a large number of parameters and other types of left–right hidden Markov models.     

Frequentist and Bayesian measures of confidence via multiscale bootstrap for testing three regions

A new computation method of frequentist p values and Bayesian posterior probabilities based on the bootstrap probability is discussed for the multivariate normal model with unknown expectation parameter vector. The null hypothesis is represented as an arbitrary-shaped region of the parameter vector. We introduce new functional forms for the scaling-law of bootstrap probability so that the multiscale bootstrap method, which was designed for a one-sided test, can also compute confidence measures of a two-sided test, extending applicability to a wider class of hypotheses. Parameter estimation for the scaling-law is improved by the two-step multiscale bootstrap and also by including higher order terms. Model selection is important not only as a motivating application of our method, but also as an essential ingredient in the method. A compromise between frequentist and Bayesian is attempted by showing that the Bayesian posterior probability with a noninformative prior is interpreted as a frequentist p value of “zero-sided” test.     

A Bayesian approach for assessing process precision based on multiple samples

Using process capability indices to quantify manufacturing process precision (consistency) and performance, is an essential part of implementing any quality improvement program. Most research works for testing the capability indices have focused on using the traditional distribution frequency approaches. Cheng and Spiring [IIE Trans. 21 (1) 97] proposed a Bayesian procedure for assessing process capability index C_p based on one single sample. In practice, manufacturing information regarding product quality characteristic is often derived from multiple samples, particularly, when a routine-based quality control plan is implemented for monitoring process stability. In this paper, we consider estimating and testing C_p with multiple samples using Bayesian approach, and propose accordingly a Bayesian procedure for capability testing. The posterior probability, p, for which the process under investigation is capable, is derived. The credible interval, a Bayesian analogue of the classical lower confidence interval, is obtained. The results obtained in this paper, are generalizations of those obtained in Cheng and Spiring [IIE Trans. 21 (1), 97]. Practitioners can use the proposed procedure to Cheng and Spiring determine whether their manufacturing processes are capable of reproducing products satisfying the preset precision requirement.     

Bayesian confidence intervals for probability of default and asset correlation of portfolio credit risk

We derive Bayesian confidence intervals for the probability of default (PD), asset correlation (Rho), and serial dependence (Theta) for low default portfolios (LDPs). The goal is to reduce the probability of underestimating credit risk in LDPs. We adopt a generalized method of moments with continuous updating to estimate prior distributions for PD and Rho from historical default data. The method is based on a Bayesian approach without expert opinions. A Markov chain Monte Carlo technique, namely, the Gibbs sampler, is also applied. The performance of the estimation results for LDPs validated by Monte Carlo simulations. Empirical studies on Standard & Poor’s historical default data are also conducted.     

A graphical method for assessing multivariate normality

Summary  In this paper we suggest a simple graphical device for assessing multivariate normality. The method is based on the characteristic that linear combinations of the sample mean and sample covariance matrix are independent if and only if the random variable is normally distributed. We demonstrate the usage of the suggested method and compare it to the classical Q-Q plot by using some multivariate data sets.     

Experimental study of a Bayesian method for daily electricity load forecasting

The method of Bayesian model discrimination is investigated for the possible contributions it may provide in the area of automatically forecasting the daily electricity demand cycle. A set of differing demand models have probabilities attached to them in such a way that these would be continuously updated with the available data and the actual forecasts obtained as expectations across all the models. Simulation experiments indicate significantly improved forecasting performance over a commonly used rescaling type of approach. Some practical issues in implementation are discussed.     

Non-probabilistic Bayesian update method for model validation

Model validation is the principal strategy to evaluate the accuracy and reliability of computational simulations. A systematic model validation procedure including uncertainty quantification, model update and prediction is described based on a non-probabilistic interval model. The crucial technical challenge in model validation is limited data, thus the non-probabilistic interval model is adopted to describe uncertain parameters. To establish the model update formula, the concepts of the interval escape rate and interval coverage rate are first described. Then, not only can the possibility of failure be estimated but also the credibility of the possibility of failure based on the proposed model validation method. The data in the validation experiment are used to update the credibility of each interval model, while the data from the accreditation experiment are used to conduct a final check of the validated models. To demonstrate that the proposed method can be applied to model validation problems successfully, a validation benchmark, the static frame challenge problem, is implemented. In addition, a practical aviation structure engineering validation problem is described. The results of these two validation problems show the feasibility and effectiveness of the proposed model validation method. The theoretical framework proposed in this paper is also suitable for model validation of computational simulations in other research fields.     

The Combination of Forecasts: a Bayesian Approach

This paper develops a general Bayesian approach to the problem of combining forecasts. This approach leads to the results of Bates and Granger in certain special cases and to a geometric averaging formula in other special cases.     

Identifying future defaulters: A hierarchical Bayesian method

Traditional methods of applying classification models into the area of credit scoring may ignore the effect from censoring. Survival analysis has been introduced with its ability to deal with censored data. The mixture cure model, one important branch of survival models, is also applied in the context of credit scoring, assuming that the study population is a mixture of never-default and will-default customers.     

Bootstrap confidence interval for a correlation curve

A correlation curve measures the strength of the association between two variables locally at different values of x. The purpose of this study is to obtain point-wise confidence intervals for a correlation curve using wild bootstrap techniques. Empirical coverage probabilities are found to be close to the specified nominal level. Bootstrapping is an attractive alternative to confidence intervals based on asymptotic expressions that have slow rate of convergence.     

A multi-criteria decision model for assessing health and self-care ability

Central European Journal of Operations Research - Population ageing together with the greater prevalence of multimorbidity add to the need for and complexity of healthcare services. This makes it...     

Building confidence in models for multiple audiences: The modelling cascade

This paper reports on a model building process developed to enable multiple audiences, particularly non-experts, to appreciate the validity of the models being built and their outcomes. The process is a four stage reversible cascade. This cascade provides a structured, auditable/transparent, formalized process from “real world” interviews generating a rich qualitative model through two intermediate steps before arriving at a quantitative simulation model. There are a number of advantages of the cascade process including; achieving comprehensiveness, developing organizational learning, testing the veracity of multiple perspectives, modeling transparency, achieving common understanding across many audiences and promoting confidence building in the models. The paper, based on extensive work with organizations, discusses both the cascade process and its inherent benefits.     

An empirical method for assessing the research relevance gap

There has been much debate on the relevance to firms of the academic research produced by business schools. However, what has not received as much attention is how the relevance of the research to businesses should be measured in a systematic and empirical way. We develop a systematic method to test for the relevance of academic research to businesses. Our method models as a vector autoregressive process the interests of the academic and practitioner communities in some new topic, as expressed by the number of articles published in the academic and the practitioner literature on that topic per calendar quarter, and then studies Granger causality between the academic and practitioner interest processes. This method can be used by academics to empirically demonstrate the impact of their intellectual contributions on practitioners and thence on the business world. We employ our approach to two relatively new and important topics, Real Options and Economic Value Added.     

Chain ladder method: Bayesian bootstrap versus classical bootstrap

The intention of this paper is to estimate a Bayesian distribution-free chain ladder (DFCL) model using approximate Bayesian computation (ABC) methodology. We demonstrate how to estimate quantities of interest in claims reserving and compare the estimates to those obtained from classical and credibility approaches. In this context, a novel numerical procedure utilizing a Markov chain Monte Carlo (MCMC) technique, ABC and a Bayesian bootstrap procedure was developed in a truly distribution-free setting. The ABC methodology arises because we work in a distribution-free setting in which we make no parametric assumptions, meaning we cannot evaluate the likelihood point-wise or in this case simulate directly from the likelihood model. The use of a bootstrap procedure allows us to generate samples from the intractable likelihood without the requirement of distributional assumptions; this is crucial to the ABC framework. The developed methodology is used to obtain the empirical distribution of the DFCL model parameters and the predictive distribution of the outstanding loss liabilities conditional on the observed claims. We then estimate predictive Bayesian capital estimates, the value at risk (VaR) and the mean square error of prediction (MSEP). The latter is compared with the classical bootstrap and credibility methods.     

关于解线性系统的一个有效的贝叶斯迭代方法

This paper concerns with the statistical methods for solving general linear systems. After a brief review of Bayesian perspective for inverse problems,a new and efficient iterative method for general linear systems from a Bayesian perspective is proposed.The convergence of this iterative method is proved,and the corresponding error analysis is studied.Finally, numerical experiments are given to support the efficiency of this iterative method,and some conclusions are obtained.     

A Monte Carlo method for an objective Bayesian procedure

This paper describes a method for an objective selection of the optimal prior distribution, or for adjusting its hyper-parameter, among the competing priors for a variety of Bayesian models. In order to implement this method, the integration of very high dimensional functions is required to get the normalizing constants of the posterior and even of the prior distribution. The logarithm of the high dimensional integral is reduced to the one-dimensional integration of a cerain function with respect to the scalar parameter over the range of the unit interval. Having decided the prior, the Bayes estimate or the posterior mean is used mainly here in addition to the posterior mode. All of these are based on the simulation of Gibbs distributions such as Metropolis' Monte Carlo algorithm. The improvement of the integration's accuracy is substantial in comparison with the conventional crude Monte Carlo integration. In the present method, we have essentially no practical restrictions in modeling the prior and the likelihood. Illustrative artificial data of the lattice system are given to show the practicability of the present procedure.