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.     

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.     

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.

     

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.     

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.     

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.     

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.     

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.     

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.     

A measure of skewness and kurtosis and a graphical method for assessing multivariate normality

Using principal components, a measure of skewness and kurtosis is developed for multivariate populations. The sample analogues of these measures are proposed as tests of multivariate normality. Also, a graphical method is presented for assessing multivariate normality.     

Mathematical foundations for a theory of confidence structures

This paper introduces a new mathematical object: the confidence structure. A confidence structure represents inferential uncertainty in an unknown parameter by defining a belief function whose output is commensurate with Neyman–Pearson confidence. Confidence structures on a group of input variables can be propagated through a function to obtain a valid confidence structure on the output of that function. The theory of confidence structures is created by enhancing the extant theory of confidence distributions with the mathematical generality of Dempster–Shafer evidence theory. Mathematical proofs grounded in random set theory demonstrate the operative properties of confidence structures. The result is a new theory which achieves the holistic goals of Bayesian inference while maintaining the empirical rigor of frequentist inference.     

Bayesian Lasso with neighborhood regression method for Gaussian graphical model

In this paper, we consider the problem of estimating a high dimensional precision matrix of Gaussian graphical model. Taking advantage of the connection between multivariate linear regression and entries of the precision matrix, we propose Bayesian Lasso together with neighborhood regression estimate for Gaussian graphical model. This method can obtain parameter estimation and model selection simultaneously. Moreover, the proposed method can provide symmetric confidence intervals of all entries of the precision matrix.     

中等稀疏条件下基因交互作用的两步Bayes方法

植物遗传与基因组学研究表明许多重要的农艺性状有影响的基因位点不是稀疏的,受到大量微效基因的影响,并且还存在基因交互项的影响.本文基于重要油料作物油菜的花期数据,研究中等稀疏条件下的基因选择问题,提出了一种两步Bayes模型选择方法.考虑基因间的交互作用,模型的维数急剧增长,加上数据结构特别,通常的变量选择方法效果不好....