SVM handles classification problem only considering samples themselves and the classification effect depends on the characteristics of the training samples but not the current information of classified problem.From th...SVM handles classification problem only considering samples themselves and the classification effect depends on the characteristics of the training samples but not the current information of classified problem.From the phenomena of data crossing in systems,this paper improves the classification effect of SVM by adding the prior probability item reflecting the classified problem information into the decision function,which fuses the Bayesian criterion into SVM.The detailed deducing process and realizing steps of the algorithm are put forward.It is verified through two instances.The results showthat the new algorithm has better effect than the traditional SVM algorithm,and its robustness and sensitivity are all improved.展开更多
A Bayesian approach is proposed for the inference of the geotechnical parameters used in slope design.The methodology involves the construction of posterior probability distributions that combine prior information on ...A Bayesian approach is proposed for the inference of the geotechnical parameters used in slope design.The methodology involves the construction of posterior probability distributions that combine prior information on the parameter values with typical data from laboratory tests and site investigations used in design.The posterior distributions are often complex,multidimensional functions whose analysis requires the use of Markov chain Monte Carlo(MCMC)methods.These procedures are used to draw representative samples of the parameters investigated,providing information on their best estimate values,variability and correlations.The paper describes the methodology to define the posterior distributions of the input parameters for slope design and the use of these results for evaluation of the reliability of a slope with the first order reliability method(FORM).The reliability analysis corresponds to a forward stability analysis of the slope where the factor of safety(FS)is calculated with a surrogate model from the more likely values of the input parameters.The Bayesian model is also used to update the estimation of the input parameters based on the back analysis of slope failure.In this case,the condition FS?1 is treated as a data point that is compared with the model prediction of FS.The analysis requires a sufficient number of observations of failure to outbalance the effect of the initial input parameters.The parameters are updated according to their uncertainty,which is determined by the amount of data supporting them.The methodology is illustrated with an example of a rock slope characterised with a Hoek-Brown rock mass strength.The example is used to highlight the advantages of using Bayesian methods for the slope reliability analysis and to show the effects of data support on the results of the updating process from back analysis of failure.展开更多
One of the main difficulties in the geotechnical design process lies in dealing with uncertainty. Uncertainty is associated with natural variation of properties, and the imprecision and unpredictability caused by insu...One of the main difficulties in the geotechnical design process lies in dealing with uncertainty. Uncertainty is associated with natural variation of properties, and the imprecision and unpredictability caused by insufficient information on parameters or models. Probabilistic methods are normally used to quantify uncertainty. However, the frequentist approach commonly used for this purpose has some drawbacks.First, it lacks a formal framework for incorporating knowledge not represented by data. Second, it has limitations in providing a proper measure of the confidence of parameters inferred from data. The Bayesian approach offers a better framework for treating uncertainty in geotechnical design. The advantages of the Bayesian approach for uncertainty quantification are highlighted in this paper with the Bayesian regression analysis of laboratory test data to infer the intact rock strength parameters σand mused in the Hoek-Brown strength criterion. Two case examples are used to illustrate different aspects of the Bayesian methodology and to contrast the approach with a frequentist approach represented by the nonlinear least squares(NLLS) method. The paper discusses the use of a Student’s t-distribution versus a normal distribution to handle outliers, the consideration of absolute versus relative residuals, and the comparison of quality of fitting results based on standard errors and Bayes factors. Uncertainty quantification with confidence and prediction intervals of the frequentist approach is compared with that based on scatter plots and bands of fitted envelopes of the Bayesian approach. Finally, the Bayesian method is extended to consider two improvements of the fitting analysis. The first is the case in which the Hoek-Brown parameter, a, is treated as a variable to improve the fitting in the triaxial region. The second is the incorporation of the uncertainty in the estimation of the direct tensile strength from Brazilian test results within the overall evaluation of the intact rock strength.展开更多
Reliability-based design (RBD) is being adopted by geotechnical design codes worldwide, and it is therefore necessary that rock engineering practice evolves to embrace RBD. This paper examines the Hoek-Brown (H-B) str...Reliability-based design (RBD) is being adopted by geotechnical design codes worldwide, and it is therefore necessary that rock engineering practice evolves to embrace RBD. This paper examines the Hoek-Brown (H-B) strength criterion within the RBD framework, and presents three distinct analyses using a Bayesian approach. Firstly, a compilation of intact compressive strength test data for six rock types is used to examine uncertainty and variability in the estimated H-B parameters m and σc, and corresponding predicted axial strength. The results suggest that within- and between-rock type variabilities are so large that these parameters need to be determined from rock testing campaigns, rather than reference values being used. The second analysis uses an extensive set of compressive and tensile (both direct and indirect) strength data for a granodiorite, together with a new Bayesian regression model, to develop joint probability distributions of m and σc suitable for use in RBD. This analysis also shows how compressive and indirect tensile strength data may be robustly used to fit an H-B criterion. The third analysis uses the granodiorite data to investigate the important matter of developing characteristic strength criteria. Using definitions from Eurocode 7, a formal Bayesian interpretation of characteristic strength is proposed and used to analyse strength data to generate a characteristic criterion. These criteria are presented in terms of characteristic parameters mk and σck, the values of which are shown to depend on the testing regime used to obtain the strength data. The paper confirms that careful use of appropriate Bayesian statistical analysis allows the H-B criterion to be brought within the framework of RBD. It also reveals that testing guidelines such as the International Society for Rock Mechanics and Rock Engineering (ISRM) suggested methods will require modification in order to support RBD. Importantly, the need to fully understand the implications of uncertainty in nonlinear strength criteria is identified.展开更多
Global spread of infectious disease threatens the well-being of human, domestic, and wildlife health. A proper understanding of global distribution of these diseases is an important part of disease management and poli...Global spread of infectious disease threatens the well-being of human, domestic, and wildlife health. A proper understanding of global distribution of these diseases is an important part of disease management and policy making. However, data are subject to complexities by heterogeneity across host classes. The use of frequentist methods in biostatistics and epidemiology is common and is therefore extensively utilized in answering varied research questions. In this paper, we applied the hierarchical Bayesian approach to study the spatial distribution of tuberculosis in Kenya. The focus was to identify best fitting model for modeling TB relative risk in Kenya. The Markov Chain Monte Carlo (MCMC) method via WinBUGS and R packages was used for simulations. The Deviance Information Criterion (DIC) proposed by [1] was used for models comparison and selection. Among the models considered, unstructured heterogeneity model perfumes better in terms of modeling and mapping TB RR in Kenya. Variation in TB risk is observed among Kenya counties and clustering among counties with high TB Relative Risk (RR). HIV prevalence is identified as the dominant determinant of TB. We find clustering and heterogeneity of risk among high rate counties. Although the approaches are less than ideal, we hope that our formulations provide a useful stepping stone in the development of spatial methodology for the statistical analysis of risk from TB in Kenya.展开更多
Proper understanding of global distribution of infectious diseases is an important part of disease management and policy making. However, data are subject to complexities caused by heterogeneities across host classes ...Proper understanding of global distribution of infectious diseases is an important part of disease management and policy making. However, data are subject to complexities caused by heterogeneities across host classes and space-time epidemic processes. This paper seeks to suggest or propose Bayesian spatio-temporal model for modeling and mapping tuberculosis relative risks in space and time as well identify risks factors associated with the tuberculosis and counties in Kenya with high tuberculosis relative risks. In this paper, we used spatio-temporal Bayesian hierarchical models to study the pattern of tuberculosis relative risks in Kenya. The Markov Chain Monte Carlo method via WinBUGS and R packages were used for simulations and estimation of the parameter estimates. The best fitting model is selected using the Deviance Information Criterion proposed by Spiegelhalter and colleagues. Among the spatio-temporal models used, the Knorr-Held model with space-time interaction type III and IV fit the data well but type IV appears better than type III. Variation in tuberculosis risk is observed among Kenya counties and clustering among counties with high tuberculosis relative risks. The prevalence of HIV is identified as the determinant of TB. We found clustering and heterogeneity of TB risk among high rate counties and the overall tuberculosis risk is slightly decreasing from 2002-2009. We proposed that the Knorr-Held model with interaction type IV should be used to model and map Kenyan tuberculosis relative risks. Interaction of TB relative risk in space and time increases among rural counties that share boundaries with urban counties with high tuberculosis risk. This is due to the ability of models to borrow strength from neighboring counties, such that nearby counties have similar risk. Although the approaches are less than ideal, we hope that our study provide a useful stepping stone in the development of spatial and spatio-temporal methodology for the statistical analysis of risk from tuberculosis in Kenya.展开更多
In this work, we seek the relationship between the order of the polynomial model and the number of knots and intervals that we need to fit the splines regression model. Regression models (polynomial and spline regress...In this work, we seek the relationship between the order of the polynomial model and the number of knots and intervals that we need to fit the splines regression model. Regression models (polynomial and spline regression models) are presented and discussed in detail in order to discover the relation. Intrinsically, both models are dependent on the linear regression model. Spline is designed to draw curves to balance the goodness of fit and minimize the mean square error of the regression model. In the splines model, the curve at any point depends only on the observations at that point and some specified neighboring points. Using the boundaries of the intervals of the splines, we fit a smooth cubic interpolation function that goes through (n + 1) data points. On the other hand, polynomial regression is a useful technique when the pattern of the data indicates a nonlinear relationship between the dependent and independent variables. Moreover, higher-degree polynomials can capture more intricate patterns, but it can also lead to overfitting. A simulation study is implemented to illustrate the performance of splines and spline segments based on the degree of the polynomial model. For each model, we compute the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) to compare the optimal polynomial order for fitting the data with the number of knots and intervals for the splines model. Both AIC and BIC can help to identify the model that best balances fit and complexity, aiming to prevent overfitting by penalizing the use of excessive parameters. We compare the results that we got from applying the polynomial regression model with the splines model results in terms of point estimates, the mean sum of squared errors, and the fitted regression line. We can say that order five of the polynomial model may be used to estimate splines with five segments.展开更多
文摘SVM handles classification problem only considering samples themselves and the classification effect depends on the characteristics of the training samples but not the current information of classified problem.From the phenomena of data crossing in systems,this paper improves the classification effect of SVM by adding the prior probability item reflecting the classified problem information into the decision function,which fuses the Bayesian criterion into SVM.The detailed deducing process and realizing steps of the algorithm are put forward.It is verified through two instances.The results showthat the new algorithm has better effect than the traditional SVM algorithm,and its robustness and sensitivity are all improved.
基金supported by the Large Open Pit Ⅱ project through contract No.019799 with the Geotechnical Research Centre of The University of Queensland,Australia and by SRK Consulting South Africa
文摘A Bayesian approach is proposed for the inference of the geotechnical parameters used in slope design.The methodology involves the construction of posterior probability distributions that combine prior information on the parameter values with typical data from laboratory tests and site investigations used in design.The posterior distributions are often complex,multidimensional functions whose analysis requires the use of Markov chain Monte Carlo(MCMC)methods.These procedures are used to draw representative samples of the parameters investigated,providing information on their best estimate values,variability and correlations.The paper describes the methodology to define the posterior distributions of the input parameters for slope design and the use of these results for evaluation of the reliability of a slope with the first order reliability method(FORM).The reliability analysis corresponds to a forward stability analysis of the slope where the factor of safety(FS)is calculated with a surrogate model from the more likely values of the input parameters.The Bayesian model is also used to update the estimation of the input parameters based on the back analysis of slope failure.In this case,the condition FS?1 is treated as a data point that is compared with the model prediction of FS.The analysis requires a sufficient number of observations of failure to outbalance the effect of the initial input parameters.The parameters are updated according to their uncertainty,which is determined by the amount of data supporting them.The methodology is illustrated with an example of a rock slope characterised with a Hoek-Brown rock mass strength.The example is used to highlight the advantages of using Bayesian methods for the slope reliability analysis and to show the effects of data support on the results of the updating process from back analysis of failure.
基金supported by the Large Open PitⅡProject through contract No.019799 with the Geotechnical Research Centre of The University of Queensland and by SRK Consulting South Africa
文摘One of the main difficulties in the geotechnical design process lies in dealing with uncertainty. Uncertainty is associated with natural variation of properties, and the imprecision and unpredictability caused by insufficient information on parameters or models. Probabilistic methods are normally used to quantify uncertainty. However, the frequentist approach commonly used for this purpose has some drawbacks.First, it lacks a formal framework for incorporating knowledge not represented by data. Second, it has limitations in providing a proper measure of the confidence of parameters inferred from data. The Bayesian approach offers a better framework for treating uncertainty in geotechnical design. The advantages of the Bayesian approach for uncertainty quantification are highlighted in this paper with the Bayesian regression analysis of laboratory test data to infer the intact rock strength parameters σand mused in the Hoek-Brown strength criterion. Two case examples are used to illustrate different aspects of the Bayesian methodology and to contrast the approach with a frequentist approach represented by the nonlinear least squares(NLLS) method. The paper discusses the use of a Student’s t-distribution versus a normal distribution to handle outliers, the consideration of absolute versus relative residuals, and the comparison of quality of fitting results based on standard errors and Bayes factors. Uncertainty quantification with confidence and prediction intervals of the frequentist approach is compared with that based on scatter plots and bands of fitted envelopes of the Bayesian approach. Finally, the Bayesian method is extended to consider two improvements of the fitting analysis. The first is the case in which the Hoek-Brown parameter, a, is treated as a variable to improve the fitting in the triaxial region. The second is the incorporation of the uncertainty in the estimation of the direct tensile strength from Brazilian test results within the overall evaluation of the intact rock strength.
文摘Reliability-based design (RBD) is being adopted by geotechnical design codes worldwide, and it is therefore necessary that rock engineering practice evolves to embrace RBD. This paper examines the Hoek-Brown (H-B) strength criterion within the RBD framework, and presents three distinct analyses using a Bayesian approach. Firstly, a compilation of intact compressive strength test data for six rock types is used to examine uncertainty and variability in the estimated H-B parameters m and σc, and corresponding predicted axial strength. The results suggest that within- and between-rock type variabilities are so large that these parameters need to be determined from rock testing campaigns, rather than reference values being used. The second analysis uses an extensive set of compressive and tensile (both direct and indirect) strength data for a granodiorite, together with a new Bayesian regression model, to develop joint probability distributions of m and σc suitable for use in RBD. This analysis also shows how compressive and indirect tensile strength data may be robustly used to fit an H-B criterion. The third analysis uses the granodiorite data to investigate the important matter of developing characteristic strength criteria. Using definitions from Eurocode 7, a formal Bayesian interpretation of characteristic strength is proposed and used to analyse strength data to generate a characteristic criterion. These criteria are presented in terms of characteristic parameters mk and σck, the values of which are shown to depend on the testing regime used to obtain the strength data. The paper confirms that careful use of appropriate Bayesian statistical analysis allows the H-B criterion to be brought within the framework of RBD. It also reveals that testing guidelines such as the International Society for Rock Mechanics and Rock Engineering (ISRM) suggested methods will require modification in order to support RBD. Importantly, the need to fully understand the implications of uncertainty in nonlinear strength criteria is identified.
文摘Global spread of infectious disease threatens the well-being of human, domestic, and wildlife health. A proper understanding of global distribution of these diseases is an important part of disease management and policy making. However, data are subject to complexities by heterogeneity across host classes. The use of frequentist methods in biostatistics and epidemiology is common and is therefore extensively utilized in answering varied research questions. In this paper, we applied the hierarchical Bayesian approach to study the spatial distribution of tuberculosis in Kenya. The focus was to identify best fitting model for modeling TB relative risk in Kenya. The Markov Chain Monte Carlo (MCMC) method via WinBUGS and R packages was used for simulations. The Deviance Information Criterion (DIC) proposed by [1] was used for models comparison and selection. Among the models considered, unstructured heterogeneity model perfumes better in terms of modeling and mapping TB RR in Kenya. Variation in TB risk is observed among Kenya counties and clustering among counties with high TB Relative Risk (RR). HIV prevalence is identified as the dominant determinant of TB. We find clustering and heterogeneity of risk among high rate counties. Although the approaches are less than ideal, we hope that our formulations provide a useful stepping stone in the development of spatial methodology for the statistical analysis of risk from TB in Kenya.
文摘Proper understanding of global distribution of infectious diseases is an important part of disease management and policy making. However, data are subject to complexities caused by heterogeneities across host classes and space-time epidemic processes. This paper seeks to suggest or propose Bayesian spatio-temporal model for modeling and mapping tuberculosis relative risks in space and time as well identify risks factors associated with the tuberculosis and counties in Kenya with high tuberculosis relative risks. In this paper, we used spatio-temporal Bayesian hierarchical models to study the pattern of tuberculosis relative risks in Kenya. The Markov Chain Monte Carlo method via WinBUGS and R packages were used for simulations and estimation of the parameter estimates. The best fitting model is selected using the Deviance Information Criterion proposed by Spiegelhalter and colleagues. Among the spatio-temporal models used, the Knorr-Held model with space-time interaction type III and IV fit the data well but type IV appears better than type III. Variation in tuberculosis risk is observed among Kenya counties and clustering among counties with high tuberculosis relative risks. The prevalence of HIV is identified as the determinant of TB. We found clustering and heterogeneity of TB risk among high rate counties and the overall tuberculosis risk is slightly decreasing from 2002-2009. We proposed that the Knorr-Held model with interaction type IV should be used to model and map Kenyan tuberculosis relative risks. Interaction of TB relative risk in space and time increases among rural counties that share boundaries with urban counties with high tuberculosis risk. This is due to the ability of models to borrow strength from neighboring counties, such that nearby counties have similar risk. Although the approaches are less than ideal, we hope that our study provide a useful stepping stone in the development of spatial and spatio-temporal methodology for the statistical analysis of risk from tuberculosis in Kenya.
文摘In this work, we seek the relationship between the order of the polynomial model and the number of knots and intervals that we need to fit the splines regression model. Regression models (polynomial and spline regression models) are presented and discussed in detail in order to discover the relation. Intrinsically, both models are dependent on the linear regression model. Spline is designed to draw curves to balance the goodness of fit and minimize the mean square error of the regression model. In the splines model, the curve at any point depends only on the observations at that point and some specified neighboring points. Using the boundaries of the intervals of the splines, we fit a smooth cubic interpolation function that goes through (n + 1) data points. On the other hand, polynomial regression is a useful technique when the pattern of the data indicates a nonlinear relationship between the dependent and independent variables. Moreover, higher-degree polynomials can capture more intricate patterns, but it can also lead to overfitting. A simulation study is implemented to illustrate the performance of splines and spline segments based on the degree of the polynomial model. For each model, we compute the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) to compare the optimal polynomial order for fitting the data with the number of knots and intervals for the splines model. Both AIC and BIC can help to identify the model that best balances fit and complexity, aiming to prevent overfitting by penalizing the use of excessive parameters. We compare the results that we got from applying the polynomial regression model with the splines model results in terms of point estimates, the mean sum of squared errors, and the fitted regression line. We can say that order five of the polynomial model may be used to estimate splines with five segments.
