The production of lithium-ion battery cells is characterized by a high degree of complexit due to numerous cause-effect relationships between process characteristics. Knowledge about the multi-stage production is spread among several experts, rendering tasks such as failure analysis challenging. In this paper, a method is presented, which includes expert knowledge acquisition in production ramp-up by combining Failure Mode and Effects Analysis (FMEA) with a Bayesian Network. We show the effectiveness of this holistic method by building up a large scale, cross-process Bayesian Failure Network in lithium-ion battery production. Using this model, we are able to conduct root cause analyses as well as analyses of failure propagation. The former support operators in identifying root causes once a cell possesses a specific failure by calculating most-probable explanations matched to the individual battery cell data. The latter enable us to analyze propagation of failures and deviations in the production chain and thus provide support for placement of quality gates, leading to a significant reduction in scrap rate. Moreover, it gives an insight into which process steps are key drivers for which final product characteristics.

Using recent developments in econometrics and computational statistics we consider the estimation of the instantaneous rate of asset return process when the underlying Data Generating Mechanism (DGM) is an Ornstein-Uhlenbeck process, driven by fractional noise, and sampled at fixed intervals of length h. To address the problem we adopt throughout the paper an exact discretization approach. This enable us to exploit the fact that a flow sampling scheme arises naturally when observing the DGM. For, while the instantaneous rate of return process is unobservable at points in time, its time integral over successive observations is observable since it equals the increment of log-prices. Exact discretization delivers an ARIMA(1,1,1) model for log-prices with a fractional driving noise. Building on the resulting exact discretization formulae and covariance function, a new Markov Chain Monte Carlo (MCMC) scheme is proposed and we examine the properties of both the time and frequency domain likelihoods / posteriors through Monte Carlo. For the exact discrete model we adopt a general sampling interval of length h. This allow us to determine the optimal choice of h independent of the sample size. An empirical application using high frequency stock price data is presented showing the relevance of aggregation over time issues in modelling asset prices.

This paper presents the results of an experimental investigation about target detecting with passive sonar in Persian Gulf. Detecting propagated sounds in the water is one of the basic challenges of the researchers in sonar field. This challenge will be complex in shallow water (like Persian Gulf) and noise less vessels. Generally, in passive sonar the targets are detected by sonar equation (with constant threshold) which increase the detection error in shallow water. Purpose of this study is proposed a new method for detecting targets in passive sonars using adaptive threshold. In this method, target signal (sound) is processed in time and frequency domain. For classifying, Bayesian classification is used and prior distribution is estimated by Maximum Likelihood algorithm. Finally, target was detected by combining the detection points in both domains using LMS adaptive filter. Results of this paper has showed that proposed method has improved true detection rate about 27% compare other the best detection method.

This study explores the relationship between job satisfaction, employee loyalty and two types of flexible employment arrangements; teleworking and flexi-time. The analysis relies on data derived by the Workplace Employee Relations Survey (WERS) in 2004 and 2011. A propensity score matching and least squares regressions are applied. Furthermore, Bayesian Networks (BN) and Directed Acyclic Graphs (DAGs) are employed in order to confirm the causality between employment types explored and the outcomes of interest. Finally, an instrumental variables (IV) approach based on the BN framework is proposed and applied in this study. The results support that there is a positive causal effect from these employment arrangements on job satisfaction and employee loyalty.

One of the major concerns in structural health management (SHM) is the early detection of a growing crack. Using this, future damage due to crack propagation can be mitigated or eliminated by implementing proper repair and maintenance. Acoustic Emission (AE) is a non-destructive testing (NDT) method with potential applications for locating and monitoring fatigue cracks. The research presented in this paper focuses on SHM using AE. A novel AE signal analysis approach was proposed in order to detect crack initiation and assess small crack lengths. Experimental investigation indicated that initiation of a crack could be identified through the statistical analysis of the resulting features of the AE signals. A probabilistic AE-based model for small fatigue crack sizing was developed and the uncertainties of the model were estimated. In addition, a probabilistic model validation approach was implemented to confirm accurate estimation of the results. The outcome of this research can be used to evaluate the integrity of structures under fatigue loading. The proposed approach can also be applied as an approach to manage health and assess prognosis of structures.

Categorical variables are common in spatial data analysis. Traditional analytical methods for deriving probabilities of class occurrence, such as kriging-family algorithms, have been hindered by the discrete characteristics of categorical fields. This study introduces the theoretical backgrounds of linear Bayesian updating (LBU) approach for spatial classification through expert system. Transition probabilities are interpreted as expert opinions for updating the prior marginal probabilities of categorical response variables. The main objective of this paper is to present the solid theoretical foundations of LBU and provide a categorical random field prediction method which yields relatively higher classification accuracy compared with conventional Markov chain random field (MCRF) approach. A real-world case study has also been carried out to demonstrate the superiority of our method. Since the LBU idea is originated from aggregating expert opinions and not restricted to conditional independent assumption (CIA), it may prove to be reasonably adequate for analyzing complex geospatial data sets, like remote sensing images or area-class maps.

Plant breeding primarily focuses on improving conventional agronomic traits, e.g. yield, quality, and resistance to biotic and abiotic stress; however, genetic improvement methods are being rapidly enhanced through genomics and phenomics. In the Genomics-Phenomics-Agronomics (GPA) paradigm, diverse research approaches have been conducted to bridge any two of these elements, and recently, all of them together. This review first highlights the progress to link i) genomics to agronomics; ii) genomics to phenomics; and iii) phenomics to agronomics. Secondly, the GPA domain is dissected into different layers, each addressing the three elements simultaneously. These dissected layers include genetic dissection through gene mapping using genome-wide association studies and genomic selection using Best Linear Unbiased Prediction, Bayesian approaches, and machine learning. The objective of the review is to help readers to grasp the core developments among the exponentially growing literature in each of these fields. Through this review, the connections among the three elements of the GPA paradigm are coherently integrated toward the prospect of sustainable development of agronomic traits through both genomics and phenomics.

When extract building from high resolution remote sensing image with meter/sub-meter accuracy, the shade of trees and interference of roads are the main factors of reducing the extraction accuracy. Proposed a Bayesian Convolutional Neural Networks(BCNET) model base on standard fully convolutional networks(FCN) to solve these problems. First take building with no shade or artificial removal of shade as Sample-A, woodland as Sample-B, road as Sample-C. Set up 3 sample libraries. Learn these sample libraries respectively, get their own set of feature vector; Mixture Gauss model these feature vector set, evaluate the conditional probability density function of mixture of noise object and roofs; Improve the standard FCN from the 2 aspect:(1) Introduce atrous convolution. (2) Take conditional probability density function as the activation function of the last convolution. Carry out experiment using unmanned aerial vehicle(UVA) image, the results show that BCNET model can effectively eliminate the influence of trees and roads, the building extraction accuracy can reach 97%.

We employ the solution concept of ex post Nash equilibrium to predict the interaction of a finite number of agents competing in a finite number of basic games simultaneously. The competition is called a multi-game. For each agent, a specific weight, considered as private information, is allocated to each basic game representing its investment in that game and the utility of each agent for any strategy profile is the weighted sum, i.e., convex combination, of its utilities in the basic games. Multi-games can model decision making in multi-environments in a variety of circumstances, including decision making in multi-markets and decision making when there are both material and social utilities for agents as, we propose, in the Prisoner's Dilemma and the Trust Game. Given a set of pure Nash equilibria, one for each basic game in a multi-game, we construct a pure Bayesian Nash equilibrium for the multi-game. We then focus on the class of so-called uniform multi-games in which each agent is constrained to play in all games the same strategy from an action set consisting of a best response per game. Uniform multi-games are equivalent to multi-dimensional Bayesian games where the type of each agent is a finite dimensional vector with non-negative components. A notion of pure type-regularity for uniform multi-games is developed and it is shown that a multi-game that is pure type-regular on the boundary of its type space has a pure ex post Nash equilibrium which is computed in constant time with respect to the number of the types and is independent of prior probability distributions. We then develop an algorithm, linear in the number of types of the agents, which tests if a multi-game is pure type-regular on the boundary of its type space in which case it returns a pure ex post Nash equilibrium for the multi-game.

In his best-selling book An Introduction to Probability Theory and its Applications, W. Feller established a way of ending the St. Petersburg Paradox by the introduction of an entrance fee, and provided it for the case in which the game is played with a fair coin. A natural generalization of his method is to establish the entrance fee for the case in which the probability of head is θ (0 < θ < 1/2). The deduction of those fees is the main result of Section 2. We then propose a Bayesian approach to the problem. When the probability of head is θ (1/2 < θ < 1) the expected gain of the St. Petersburg Game is finite, therefore there is no paradox. However, if one takes θ as a random variable assuming values in (1/2,1) the paradox may hold, what is counter-intuitive. On Section 3 we determine a necessary and sufficient condition for the absence of paradox on the Bayesian approach and on Section 4 we establish the entrance fee for the case in which θ is uniformly distributed in (1/2,1), for in this case there is paradox.

Bayesian analysis for masked data under competing risk frameworks is studied for the purpose of assessing the impact of covariates on the hazard functions when the failure time is exactly observed for some subjects but only known to lie in an interval of time for the remaining subjects. Such data, known as partly interval-censored data, usually result from periodic inspection. Dirichlet and Gamma processes are assumed as priors for masking probabilities and baseline hazards. The Markov Chain Monte Carlo (MCMC) technique is employed for the implementation of the Bayesian approach. The effectiveness of the proposed model is tested through numerical studies, including simulated and real data sets.

The problem of event detection in general noisy signals arises in many applications; usually, either a functional form for the event is available, or a previous annotated sample with instances of the event that can be used to train a classification algorithm. There are situations, however, where neither functional forms nor annotated samples are available; then it is necessary to apply other strategies to separate and characterize events. In this work, we analyze 15 minute-long samples of an acoustic signal, and are interested in separating sections, or segments, of the signal which are likely to contain significative events. For that, we apply a sequential algorithm with the only assumption that an event alters the energy of the signal. The algorithm is entirely based on Bayesian methods.

The comprehensibility of good predictive models learned from high-dimensional gene expression data is attractive because it can lead to biomarker discovery. Several good classifiers provide comparable predictive performance but differ in their abilities to summarize the observed data. We extend a Bayesian Rule Learning (BRL-GSS) algorithm, previously shown to be a significantly better predictor than other classical approaches in this domain. It searches a space of Bayesian networks using a decision tree representation of its parameters with global constraints, and infers a set of IF-THEN rules. The number of parameters and therefore the number of rules are combinatorial to the number of predictor variables in the model. We relax these global constraints to a more generalizable local structure (BRL-LSS). BRL-LSS entails more parsimonious set of rules because it does not have to generate all combinatorial rules. The search space of local structures is much richer than the space of global structures. We design the BRL-LSS with the same worst-case time-complexity as BRL-GSS while exploring a richer and more complex model space. We measure predictive performance using Area Under the ROC curve (AUC) and Accuracy. We measure model parsimony performance by noting the average number of rules and variables needed to describe the observed data. We evaluate the predictive and parsimony performance of BRL-GSS, BRL-LSS and the state-of-the-art C4.5 decision tree algorithm, across 10-fold cross-validation using ten microarray gene-expression diagnostic datasets. In these experiments, we observe that BRL-LSS is similar to BRL-GSS in terms of predictive performance, while generating a much more parsimonious set of rules to explain the same observed data. BRL-LSS also needs fewer variables than C4.5 to explain the data with similar predictive performance. We also conduct a feasibility study to demonstrate the general applicability of our BRL methods on the newer RNA sequencing gene-expression data.

This paper develops Bayesian inference in reliability of a class of scale mixtures of log-normal failure time (SMLNFT) models with stochastic (or uncertain) constraint in their reliability measures. The class is comprehensive and includes existing failure time (FT) models (such as log-normal, log-Cauchy, and log-logistic FT models) as well as new models that are robust in terms of heavy-tailed FT observations. Since classical frequency approaches to reliability analysis based on the SMLNFT model with stochastic constraint are intractable, the Bayesian method is pursued utilizing a Markov chain Monte Carlo (MCMC) sampling based approach. This paper introduces a two-stage maximum entropy (MaxEnt) prior, which elicits a priori uncertain constraint and develops Bayesian hierarchical SMLNFT model by using the prior. The paper also proposes an MCMC method for Bayesian inference in the SMLNFT model reliability and calls attention to properties of the MaxEnt prior that are useful for method development. Finally, two data sets are used to illustrate how the proposed methodology works.

This article shortly introduces Gaussian processes as a new approach for modelling time series in the field of blazar physics. In the second part of the paper, recent results from an application of GP modelling to the multi-wavelength light curves of the blazar PKS 1502+106 are discussed.

An approach to the assessment of probabilistic inference is described which quantifies the performance on the probability scale. From both information and Bayesian theory, the central tendency of an inference is proven to be the geometric mean of the probabilities reported for the actual outcome and is referred to as the “Accuracy.” Upper and lower error bars on the accuracy are provided by the arithmetic mean and the -2/3 mean. The arithmetic is called the “Decisiveness” due to its similarity with the cost of a decision and the -2/3 mean is called the “Robustness”, due to its sensitivity to outlier errors. Visualization of inference performance is facilitated by plotting the reported model probabilities versus the histogram calculated source probabilities. The visualization of the calibration between model and source is summarized on both axes by the arithmetic, geometric, and -2/3 means. From information theory, the performance of the inference is related to the cross-entropy between the model and source distribution. Just as cross-entropy is the sum of the entropy and the divergence; the accuracy of a model can be decomposed into a component due to the source uncertainty and the divergence between the source and model. Translated to the probability domain these quantities are plotted as the average model probability versus the average source probability. The divergence probability is the average model probability divided by the average source probability. When an inference is over/under-confident, the arithmetic mean of the model increases/decreases, while the -2/3 mean decreases/increases, respectively.

Knowledge-based approaches are useful alternatives to predict the Failure Probability (FP) coping with the insufficient data, process integrity, and complexity issue in manufacturing systems. This study proposes a Fault Tree Analysis (FTA) approach as proactive knowledge-based technique to estimate the FP based maintenance planning with subjective information from domain experts. However, the classical-FTA still suffers from uncertainty and static structure limitations which poses a substantial dilemma in predicting FP. To deal with the uncertainty issues, a Fuzzy-FTA (FFTA) model was developed by statistical analysing the effective attributes such as experts' trait impacts, scales variation and, assorted membership and defuzzification functions. Besides, a Bayesian Network (BN) theory was conducted to overcome the static limitation of classical-FTA. The results of FFTA model revealed that the changes in decision attributes were not statistically significant on FP variation while BN model considering conditional rules to reflect the dynamic relationship between events had more impact on predicting FP. After all, the integrated FFTA-BN model was used in the optimization model to find the optimal maintenance intervals according to estimated FP and the total expected cost. As a practical example, the proposed model was implemented in a semi-automatic filling system in an automotive production line. The outcomes could be useful for upgrading the availability and safety of complex equipment in manufacturing systems.

We present a first formal analysis of specific and complete local integration. Complete local integration was previously proposed as a criterion for detecting entities or wholes in distributed dynamical systems. Such entities in turn were conceived to form the basis of a theory of emergence of agents within dynamical systems. Here, we give a more thorough account of the underlying formal measures. The main contribution is the disintegration theorem which reveals a special role of completely locally integrated patterns (what we call ι-entities) within the trajectories they occur in. Apart from proving this theorem we introduce the disintegration hierarchy and its refinement-free version as a way to structure the patterns in a trajectory. Furthermore we construct the least upper bound and provide a candidate for the greatest lower bound of specific local integration. Finally, we calculate the i-entities in small example systems as a first sanity check and find that ι-entities largely fulfil simple expectations.

Fitting parameters on spline curves produces more parsimonious models, maintaining fit quality. Smoothing the splines reduces predictive variance. Individual splines are often fit by type of variable, e.g., in age-period-cohort models. Linear and cubic splines are most common. Several smoothing criteria have been used for parameter curves, with cubic splines fit by constraining the integral of splines’ squared-second derivatives popular recently. Constraining the sum of squared second differences for linear splines is analogous. Generally the degree of smoothing is selected using cross-validation. Known spline dummy-variable matrices allow regression estimation of splines, with smoothing done via constrained regression. Smoothing criteria based on sums of squares or absolute values of parameters, as in ridge regression or LASSO, improves predictive accuracy and produces splines similar to smoothing by second-derivative constraints. Variables with very low t-statistics represent points where curve-shapes barely change. Eliminating those variables leaves knots concentrated where spline shapes change. A Bayesian version of this puts shrinkage priors on spline parameters. This yields realistic joint parameter distributions, avoids problems associated with using cross-validation for parameter estimation, and readily expands to non-linear modeling, such as interactions among variable types. Regularized regression and Bayesian spline methods are compared for two example datasets.

We proposed a so-called Beta Kumaraswamy Burr Type X distribution which gives the extension of the Kumaraswamy-G class of family distribution. Some properties of this proposed model were provided, like: the expansion of densities and quantile function. We considered the Bayes and maximum likelihood methods to estimate the parameters and also simulate the model parameters to validate the methods based on different set of true values. Some real data sets were employed to show the usefulness and flexibility of the model which serves as generalization to many sub-models in the field of engineering, medical, survival and reliability analysis.

Historical forest management practices in the southwestern US have left forests prone to high intensity, stand-replacement fires. Effective management to reduce the cost and impact of forest-fire management and allow fires to burn freely without negative impact depends on detailed knowledge of stand composition, in particular, above-ground biomass (AGB). Lidar-based modeling techniques provide opportunities to reduce costs and increase ability of managers to monitor AGB and other forest metrics. Using Bayesian Model Averaging (BMA), we develop a regionally applicable lidar-based statistical model for Ponderosa pine and mixed conifer forest systems of the southwestern USA, using previously collected field data. The selected regional model includes a mid and low canopy height metric, a canopy cover, and height distribution term. It explains 72% of the variability in field estimates of AGB, and the RMSE of the two independent validation data sets are 23.25 and 32.82 Mg/ha. The regional model developed is structured in accordance with previously described models fit to local data, and performs equivalently to models designed for smaller scale application. Developing regional models for broad scale application provides a cost-effective, robust approach for managers to monitor and plan adaptively at the landscape scale.

The performance in terms of minimal Bayes’ error probability for detection of a high-dimensional random tensor is a fundamental under-studied difficult problem. In this work, we consider two Signal to Noise Ratio (SNR)-based detection problems of interest. Under the alternative hypothesis, i.e., for a non-zero SNR, the observed signals are either a noisy rank-R tensor admitting a Q-order Canonical Polyadic Decomposition (CPD) with large factors of size Nq  R, i.e, for 1  q  Q, where R, Nq ! ¥ with R1/q/Nq converge towards a finite constant or a noisy tensor admitting TucKer Decomposition (TKD) of multilinear (M1, . . . ,MQ)-rank with large factors of size Nq  Mq, i.e, for 1  q  Q, where Nq,Mq ! ¥ with Mq/Nq converge towards a finite constant. The detection of the random entries (coefficients) of the core tensor in the CPD/TKD is hard to study since the exact derivation of the error probability is mathematically intractable. To circumvent this technical difficulty, the Chernoff Upper Bound (CUB) for larger SNR and the Fisher information at low SNR are derived and studied, based on information geometry theory. The tightest CUB is reached for the value minimizing the error exponent, denoted by s?. In general, due to the asymmetry of the s-divergence, the Bhattacharyya Upper Bound (BUB) (that is, the Chernoff Information calculated at s? = 1/2) can not solve this problem effectively. As a consequence, we rely on a costly numerical optimization strategy to find s?. However, thanks to powerful random matrix theory tools, a simple analytical expression of s? is provided with respect to the Signal to Noise Ratio (SNR) in the two schemes considered. A main conclusion of this work is that the BUB is the tightest bound at low SNRs. This property is, however, no longer true for higher SNRs.

The evolution of manual grooming and its implications have received little attention in the quest to understand the origins of simian primates and their social and technical intelligence. All simians groom manually, whereas prosimians groom orally despite comparable manual dexterity between some members of the two groups. Simians also exhibit a variable propensity for the manipulation of inanimate, non-food objects, which has culminated in tool making and tool use in some species. However, lemuriform primates also seem capable of tool use with training. Furthermore, lemuriforms appear to understand the concept of a tool and use their own body parts as “tools”, despite not using inanimate objects. This suggests that prosimian primates are pre-adapted for proprioceptive object manipulation and tool use, but do not express these cognitive abilities by default. This essay explores the paleontological, anatomical, cognitive, ethological, and neurological roots of these abilities and attempts to explain this behavioural divide between simians and prosimians. Common misconceptions about early primate evolution and captive behaviours are addressed, and chronological inconsistencies with Machiavellian Intelligence are examined. A “licking to picking” hypothesis is also proposed to explain a potential link between manual grooming and object manipulation, and to reconcile the inconsistencies of Machiavellian Intelligence. Bayesian decision theory, the evolution of the parietal cortex and enhanced proprioception, and analogies with behavioural changes resulting from artificial selection may help provide new insights into the minds of both our primate kin and ourselves.

Construction projects are usually operating in a complex and dynamic environment in which the accumulation of many interrelated factors causes high uncertainty. Construction projects are complex and frequently involve substantial uncertainties including process complicatedness, intricate organization structure, dynamic environment, and financial strain. The study aims to categorize the influencing factors into three groups, namely construction project system, economic-market climate, and external environment. It attempts to adopt a novel analysis tool to examine the relationship between the project cost and multiple influencing factors by using Bayesian SEM. While the Bayesian SEM method has been receiving increasing attention in exploring the relationship between latent variables, construction studies still heavily rely on the covariance-based SEM approach. This study introduces several advantages of Bayesian SEM that make it more flexible and powerful than covariance-based SEM and provides the foundation of Bayesian SEM estimation and inference by illustrating this method in a project cost application.

One of the most common and important destructive attacks on the victim system is Advanced Persistent Threats (APT)-attack. The APT attacker can achieve his hostile goals by obtaining information and gaining financial benefits regarding the infrastructure of a network. One of the solutions to detect a secret APT attack is using network traffic. Due to the nature of the APT attack in terms of being on the network for a long time and the fact that the network may crash because of high traffic, it is difficult to detect this type of attack. Hence, in this study, machine learning methods such as C5.0 decision tree, Bayesian network and deep neural network are used for timely detection and classification of APT-attacks on the NSL-KDD data set. Moreover, 10-fold cross validation method is used to experiment these models. As a result, the accuracy (ACC) of the C5.0 decision tree, Bayesian network and 6-layer deep learning models is obtained as 95.64%, 88.37% and 98.85%, respectively, and also, in terms of the important criterion of the false positive rate (FPR), the FPR value for the C5.0 decision tree, Bayesian network and 6-layer deep learning models is obtained as 2.56, 10.47 and 1.13, respectively. Other criterions such as sensitivity, specificity, accuracy, false negative rate and F-measure are also investigated for the models, and the experimental results show that the deep learning model with automatic multi-layered extraction of features has the best performance for timely detection of an APT-attack comparing to other classification models.

This article describes the use of the PI ProcessBook software tool for visualization and indirect monitoring of occupancy of SHC rooms from the measured operational and technical quantities for monitoring of daily living activities for support of independent life of elderly persons. The proposed method for data processing (predicting the CO₂ course using neural networks from the measured temperature indoor T_i (°C), temperature outdoor T_o (°C) and the relative humidity indoor rHi (%)) was implemented, verified and compared in MATLAB SW tool and IBM SPSS SW tool with IoT platform connectivity. Within the proposed method, the Stationary Wavelet Transform de noising algorithm was used to remove the noise of the resulting predicted course. In order to verify the method, two long-term experiments were performed, (specifically from February 8 to February 15, 2015, from June 8 to June 15, 2015) and two short-term experiments (from February 8, 2015 and from June 8, 2015). For the best results of the trained ANN BRM within the prediction of CO₂, the correlation coefficient R for the proposed method was up to 90%. The verification of the proposed method confirmed the possibility to use the presence of persons of the monitored SHC premises for rooms ADL monitoring.

The Autopoiesis and Cognition Theory (ACT), by Maturana and Varela, based on the notions of Biological Closure and Structural Coupling, is a well-known theory on how to understand biological organization [1, 2, 3]. Although, for example, the Free Energy Principle framework evokes some entailments of autopoiesis in a more formal setting [4, 5]; and ACT has been used in many fields, its impact has been restricted because it lacks quantitative analysis. Here we present a theoretical framework grounded in accepted and well-developed ideas from Mathematics and Physics which advance the understanding of the Principles of Biological Organization under the guidance of Biological Closure and Structural Coupling. The disciplines of Differential Geometry/Topology, Mechanics and Complex Dynamical Systems provide a powerful, elegant, and well-established body of knowledge to support our Biological Organization Principles (BOP) framework. In particular, Stochastic Mechanics and KAM theory (from Kolmogorov, Arnold and Moser theorem) allow us to develop, using the notions of Biological Closure and Structural Coupling, a central core of BOP termed Dynamical Closure Mechanism. Under the proposed framework, a wide variety of bio- logical phenomena can be understood, shedding new light on biological explanations. However, an understanding of biological organization may require the re-evaluation of dogmas on how we think on biology as it seems inescapable that what is needed is an integration of analysis and notions derived from mathematics, physics, and biology to generate a new landscape of ideas.

Especially for drought periods, the higher the accuracy of reservoir inflow forecasting, the more reliable the water supply from a dam. The article focuses on probabilistic forecasting of seasonal inflow to reservoirs and determines estimates from the probabilistic seasonal inflow according to drought forecast results. The probabilistic seasonal inflow was forecasted by a copula-based Bayesian network employing a Gaussian copula function. Drought forecasting was performed by calculation of the standardized streamflow index value. The calendar year is divided into four seasons; the total inflow volume of water to a reservoir for a season is referred to as the seasonal inflow. Seasonal inflow forecasting curves conforming to drought stages produce estimates of probabilistic seasonal inflow according to the drought forecast results. The forecasted estimates of seasonal inflow were calculated by using the inflow records of Soyanggang and Andong dams in the Republic of Korea. Under the threshold probability of drought occurrence ranging from 50 to 55 %, the forecasted seasonal inflows reasonably matched critical drought records. Combining the drought forecasting with the seasonal inflow forecasting may produce reasonable estimates of drought inflow from the probabilistic forecasting of seasonal inflow to a reservoir.

Advances in remote sensing have led to use of satellite-derived rainfall products to complement the sparse rain gauge data. Although globally derived and some regional bias corrected, these products often show large discrepancies with ground measurements attributed to local and external factors that require systematic consideration. Decreasing rain gauge network however inhibits continuous validation of these products. We propose to deal with this problem by the use of Bayesian approach to merge the existing historical rain gauge information to create a consistent satellite rainfall data that can be used for climate studies. Monthly Bayesian bias correction is applied to the Climate Hazards Group Infrared Precipitation with Stations (CHIRPS v2) data to reduce systematic errors using a corresponding gridded (0.05°) rain gauge data over East Africa for a period of 33 (1981–2013) years of which 22 years are utilized to derive error fields which are then applied to an independent CHIRPS data for 11 years for validation. The bias correction is spatially and temporally assessed during the rainfall wet months of March-May (MAM), June-August (JJA) and October–December (OND) in East Africa. Results show significant reduction of systematic errors at both monthly and yearly scales and harmonization of their cumulative distributions. Monthly statistics showed a reduction of RMSD (29–56)% and MAE (28–60)% and an increase of correlations (2–32) %, while yearly ones showed reductions of RMSD (9-23)%, and MAE (7–27)% and increase of correlations (4–77)% for MAM months, reduction of RMSD (15–35)% and MAE (16–41)% and increase in correlations (5–16)% for JJA months, and reduction of RMSD (3–35)% and MAE (9–32)% and increase of correlations (3–65)% for OND months. Systematic errors of corrected data were influenced by local processes especially over Lake Victoria and high elevated areas. Large-scale circulations induced errors were mainly during JJA and OND rainfall seasons and were reduced by the separation of anomalous years during training. The proposed approach is recommended for generating long-term data for climate studies where consistencies of errors can be assumed.

The present reality about the Nigerian economy calls for investment and development in the non-oil sector. This becomes necessary as a result of fall in the oil price in the global market. This paper examined the Bayesian Vector Autoregression (BVAR) modeling and forecasting of the dynamic interrelationship between Economic growth and revenue from the oil and non-oil sectors in Nigeria. To achieve this, annual data on Gross Domestic Product (GDP), revenue from oil and non-oil sectors were collected from Central Bank of Nigeria (CBN) bulletin, the sample from 1981 to 2008 was used for analysis, while sample from 2009 to 2014 was used for model validation. Six (6) versions of Sims-Zha BVAR models were compared for out-of-sample forecast, the result revealed the superiority of the BVAR6 model over the other BVAR models. Lastly, evidence from the decomposition forecast errors revealed that revenue of oil sector contributed 7.69% to GDP while revenue from non-oil sector contributed 0.12% to GDP in Nigeria. This paper therefore recommended that the present government should encourage investment that is geared toward development in the non-oil sector, of which it has the capacity to improve the Economic growth of the Nigerian economy.

Many scientific and technological problems are studied with the help of computer codes that simulate the phenomena of interest rather than via traditional laboratory experiments. Such models play an important role in neuroscience where they are used to mimic brain function from the sub-cellular to the macroscopic level. Exploration with computer models carries with it a number of statistical challenges: where to sample the input space for the simulator, how to make sense of the data that is generated, how to estimate unknown parameters in the model, how to validate a model. The simulator setting also has some unique problems and possibilities. This review paper describes statistical research on these issues and how that work might be applied to neural simulations.

A vast majority of the countries is under the economic and health crises due to the current epidemic of coronavirus disease 2019 (COVID-19). The present study analyzes the COVID-19 using time series, which is an essential gizmo for knowing the enlargement of infection and its changing behavior, especially the trending model. We have considered an autoregressive model with a non-linear time trend component that approximately converted into the linear trend using the spline function. The spline function split the COVID-19 series into different piecewise segments between respective knots and fitted the linear time trend. First, we obtain the number of knots with its locations in the COVID-19 series and then the estimation of the best-fitted model parameters are determined under Bayesian setup. The results advocate that the proposed model/methodology is a useful procedure to convert the non-linear time trend into a linear pattern of newly coronavirus case for various countries in the pandemic situation of COVID-19.

The goal of VAR or BVAR is the characterization of the dynamics and endogenous relationships among time series. Also the VAR models are known for their applications to forecasting and policy analysis. This paper compare the performance of VAR and Sims-Zha Bayesian VAR models when the multiple time series are jointly influenced by different levels of collinearity and autocorrelation in the short term (T=16, 32, 64 and 128). Five levels (-0.9,-0.5, 0,+0.5,+0.9) of collinearity and autocorrelation were considered and the results from the simulation study revealed that VAR(2) model dominated for no and moderate levels of autocorrelation (-0.5, 0, +0.5) irrespective of the collinearity level except in few cases when T=16. While the BVAR models dominated for high autocorrelation levels (-0.9 and +0.9) irrespective of the collinearity level except in few cases when T=128. The performance of the models varies at different levels of the collinearity and autocorrelated error, and also varies with the short term periods. Furthermore, the values of the RMSE and MAE criteria decrease as a result of increase in the time series length. In conclusion, the performance of the forecasting models depend on the time series data structure and the time series length. It is therefore recommended that the data structure and series length should be considered in using an appropriate model for forecasting.

This paper investigates an inverse problem of option pricing in the extended Black--Scholes model. We identify the model coefficients from the measured data and attempt to find arbitrage opportunities in financial markets using a Bayesian inference approach. The posterior probability density function of the parameters is computed from the measured data. The statistics of the unknown parameters are estimated by Markov Chain Monte Carlo (MCMC), which explores the posterior state space. The efficient sampling strategy of MCMC enables us to solve inverse problems by the Bayesian inference technique. Our numerical results indicate that the Bayesian inference approach can simultaneously estimate the unknown drift and volatility coefficients from the measured data.

Cartels cause tremendous damage to the market economy and disadvantage consumers by causing higher prices and lower quality; moreover, they are difficult to detect. We need to prevent them by scientific analysis, which includes the determination of an indicator to explain antitrust enforcement. Particularly, the probability of cartel penalization is a useful indicator for the evaluation of the competition enforcement. This study is to estimate the probability of cartel penalization by using a Bayesian approach. In the empirical study, the probability of cartel penalization is estimated by Bayesian approach from cartel data of Department of Justice in United States from 1970 to 2009. The probability of cartel penalization is seen to be sensitive to change of competition law and the results shows the usefulness of higher interpretation than other research. The result of the policy simulation shows how effective the leniency program is. From this estimation, antitrust enforcement is evaluated, and thereby, can be improved.

This study is framed on the Information Theory as a constructive criterion to generate probabilistic distributions –through the elaboration of Bayesian Networks- and to reduce the uncertainty in the occurrence and relationship between two key psychological variables associated with the sports’ performance: Self-Determined Motivation and Competitive Anxiety. We analyzed 674 universitary students/athletes who competed in the 2017 Universitary Games (Universiade) in México, from 44 universities, with an average age of 21 years old (SD = 2.07), and with a sportive experience of 8.61 years of average (SD = 5.15). Methods: Regarding the data analysis, first of all a CHAID algorithm was carried out for to know the independence links among variables, and then two Bayesian networks (BN) were elaborated. The validation of the BN revealed AUC values ranging from 0.5 to 0.92. Subsequently, various instantations were carried out with hypothetical values applied to the “bottom” variables. Results showed two probability trees that have Extrinisic Motivation and Amotivation at the top, while the anxiety/activation due to the worry for performance was at the bottom of probabilities. The instantiations carried out support the existence of these probabilistic relationships, demonstrating the little influence on the competition anxiety generated by the intrinsic motivation. In conclusion, the reduction of the uncertainty made up by the restricted BN may aloe to re-introduce Information Theory principles in psychosocial studies, allowing authors to obtain useful probabilities values upon target psychological variables related with sportive performance.

The aim of the paper is to identify the most likely factors that determine the demand for Renewa-ble Energy Consumption (R.E.C.) in European countries. Although in Europe a high environmen-tal awareness is omnipresent, countries differ in scope and share of R.E.C. due to historical ener-getic policies and dependencies, investments into renewable and traditional energetic sectors, R&D development, structural changes required by energetic policy change, and many other fac-tors. The study refers to a set of macroeconomic, institutional, and social factors affecting energetic renewable policy and R.E.C. in selected European countries in two points of time: i.e., before and after the Paris Agreement. The Bayesian Average Classical Estimates (BACE) is applied to indicate the most likely factors affecting R.E.C. in 2015 and 2018. The comparison of the results reveals that the G.D.P. level, nuclear and hydro energy consumption were the determinants significant in both analyzed years. Furthermore, it became clear that in 2015 the R.E.C. depended strongly on the energy consumption structure, while in 2018, the foreign direct investment and trade openness played their role in increasing renewable energy consumption. The direction of changes is positive and complies with sustainable development goals (S.D.G.s).

Determination of indoor location based on fine time measurement (FTM) of the round trip time (RTT) of a signal between an initiator (smartphone) and a responder (Wi-Fi access point) enables a number of applications. However, the accuracy currently attainable — standard deviations of 1–2 meter in distance measurement under favorable circumstances — limits the range of possible application. A first responder, for example, may not be able to unequivocally determine on which floor someone in need of help is in a multi-story building. The error in location depends on several factors, including the bandwidth of the RF signal, delay of the signal due to the high relative permittivity of construction materials, and the geometry-dependent “noise gain” of location determination. Errors in distance measurements have unusual properties that are exposed here for the first time. Improvements in accuracy depend on understanding all of these error sources. This paper introduces “frequency diversity,” a method for doubling the accuracy of indoor location determination using weighted averages of measurements with uncorrelated errors obtained in different channels. The properties of this method are verified experimentally with a range of responders. Finally, different ways of using the distance measurements to determine indoor location are discussed and the Bayesian grid update method shown to be more useful than others, given the non- Gaussian nature of the measurement errors.

We proposed a so-called Beta Kumaraswamy Burr Type X distribution which gives the extension of the Kumaraswamy-G class of family distribution. Some properties of this proposed model were provided, like: the expansion of densi- ties and quantile function. We considered the Bayes and maximum likelihood methods to estimate the parameters and also simulate the model parameters to validate the methods based on dierent set of true values. Some real data sets were employed to show the usefulness and exibility of the model which serves as generalization to many sub-models in the elds of engineering, medical, survival and reliability analysis.

Genomic selection (GS) can accelerate variety improvement when training set (TS) size, and its relationship with the breeding set (BS) are optimized for prediction accuracies (PA) of genomic prediction (GP) models. Sixteen GP algorithms were run on phenotypic best linear unbiased predictors (BLUPs) and estimators (BLUEs) of resistance to both fall armyworm (FAW) and maize weevil (MW) in a tropical maize panel. For MW resistance, 37% of the panel was the TS, and BS was the remainder whilst for FAW, random-based training sets (RBTS) and pedigree-based training sets (PBTS) were designed. PAs achieved with BLUPs varied from 0.66 to 0.82 for MW resistance traits, and, for FAW resistance, 0.694 to 0.714 for RBTS of 37%, and 0.843 to 0.844 for RBTS of 85%, and, these were at least two-fold those from BLUEs. For PBTS, FAW resistance PAs were generally higher than those for RBTS, except for one dataset. GP models generally showed similar PAs across individual traits whilst the TS designation was determinant since a positive correlation (R=0.92***) between TS size and PAs was observed for RBTS and, for the PBTS, it was negative (R=0.44**). This study pioneers the use of GS for maize resistance to insect pests in sub-Saharan Africa.

The hourly and daily rainfall data which is spatially distributed are required as an input for run-off rain model. Furthermore, the run-off rain model is used to detect early flooding. The daily and hourly rainfall data have characteristics that most of data are zero. Therefore we need a model which can capture the phenomenon. A time series model involving location, which is a model that can be developed to approach the daily and hourly rainfall data, we can call Vector Autoregressive (VAR) model. The VAR model allows us for modeling rainfall data in several areas. However, in certain conditions the VAR model often occurs over-parameterization and reduces degrees of freedom. The aim of this study is to compare the VAR model with Maximum Likelihood Estimator (MLE) and Bayesian to hourly and daily rainfall data in SampeanWatershed of Bondowoso. The results showed that the hourly and daily rainfall data are fitted to VAR process of orde 5 and 1 respectively. Based on the AIC and SBC values indicate that the Bayesian is better than the MLE method. The Bayesian is able to predict parameters by producing a smaller variance covariance matrix than the MLE.

The adoption of Machine Learning (ML) for building emulators for complex physical processes has seen an exponential rise in the recent years. While neural networks are good function approximators, optimizing the hyper-parameters of the network to reach a global minimum is not trivial, and often needs human knowl- edge and expertise. In this light, automatic ML or autoML methods have gained large interest as they automate the process of network hyper-parameter tuning. In addition, Neural Architecture Search (NAS) has shown promising outcomes for improving model performance. While autoML methods have grown in popularity for image, text and other applications, their effectiveness for high-dimensional, complex scientific datasets remains to be investigated. In this work, a data driven emulator for turbulence closure terms in the context of Large Eddy Simulation (LES) models is trained using Artificial Neural Networks and an autoML frame- work based on Bayesian Optimization, incorporating priors to jointly optimize the hyper-parameters as well as conduct a full neural network architecture search to converge to a global minima, is proposed. Additionally, we compare the effect of using different network weight initializations and optimizers such as ADAM, SGDM and RMSProp, to explore the best performing setting. Weight and function space similarities during the optimization trajectory are investigated, and critical differences in the learning process evolution are noted and compared to theory. We observe ADAM optimizer and Glorot initialization consistently performs better, while RMSProp outperforms SGDM as the latter appears to have been stuck at a local minima. Therefore, this autoML BayesOpt framework provides a means to choose the best hyper-parameter settings for a given dataset.

Bayesian nonlinear mixed effects models for data in the form of continuous, repeated measurements from a population, also known as Bayesian hierarchical nonlinear models, are a popular platform for analysis when interest focuses on individual specific characteristics and relevant uncertainty quantification. Due to the limitation of computational power, this framework was relatively dormant until the late 1980s, but in recent years, the statistical research community saw vigorous development of new methodological and computational techniques for these models, the emergence of software, and wide application of the models in numerous industrial and academic fields. This article presents an overview of the formulation, interpretation, and implementation of Bayesian nonlinear mixed effects models and surveys recent advances and applications.

Equilibrium climate sensitivity (ECS) is the change in global mean temperature expected to result from doubling atmospheric CO2 concentration from pre-industrial levels. Extensive research during the past 40 years has not reduced the uncertainty associated with ECS. Sherwood et al. [1] applied Bayesian statistics to evidence from climate-process physics, historical observations and earlier proxies to reduce the range of ECS from 1.5 – 4.5 K to 2.6 – 4.1 K. This paper examines their methods and many of the assumptions they made. It also evaluates two additional periods in the Holocene to show that factors other than CO2 drove recent climate change. It identifies potential systematic errors resulting from adding non-equilibrium short-term adjustments to the radiative forcing of greenhouse gases and from underestimating the effects of solar irradiance, ocean currents and aerosols. These factors have resulted in estimates of the forcing by CO2 that far exceed the apparent effects in paleoclimate data.

The general objective of the study is to analyze food insecurity determinants mainly related to food price at macro level. The study used a time series data from the year 1988 quarter one to 2018 quarter four. Our empirical finding on the macroeconomic model show the short term consequence of food price change is positive with food security (negative with food insecurity). The long term consequence of food price change is, however, negative with food security (positive with food insecurity) at macro level. From the variance decomposition, it can be learned that the speed by which the undernourished people decrease in number in the short term because of the inflation happened is higher than the speed by which the number of undernourished people increase in the long run. At macro level, concern should be made to the supply side for food security improvement.

The recent Google’s claim on breakthrough in quantum computing is a gong signal for further analysis of foundational roots of (possible) superiority of some quantum algorithms over the corresponding classical algorithms. This note is a step in this direction. We start with critical analysis of rather common reference to entanglement and quantum nonlocality as the basic sources of quantum superiority. We elevate the role of the Bohr’s principle of complementarity¹ (PCOM) by interpreting the Bell-experiments as statistical tests of this principle. (Our analysis also includes comparison of classical vs genuine quantum entanglements.) After a brief presentation of PCOM and endowing it with the information interpretation, we analyze its computational counterpart. The main implication of PCOM is that by using the quantum representation of probability, one need not compute the joint probability distribution (jpd) for observables involved in the process of computation. Jpd’s calculation is exponentially time consuming. Consequently, classical probabilistic algorithms involving calculation of jpd for n random variables can be over-performed by quantum algorithms (for big values of n). Quantum algorithms are based on quantum probability calculus. It is crucial that the latter modifies the classical formula of total probability (FTP). Probability inference based on the quantum version of FTP leads to constructive interference of probabilities increasing probabilities of some events. We also stress the role the basic feature of the genuine quantum superposition comparing with the classical wave superposition: generation of discrete events in measurements on superposition states. Finally, the problem of superiority of quantum computations is coupled with the quantum measurement problem and linearity of dynamics of the quantum state update.

Chikungunya virus (CHIKV) and Zika virus (ZIKV) have recently emerged as global infections with consequential disability adjusted life years (DALYs) and economic burden. This study aimed to explore the spatiotemporal heterogeneity in the occurrence of CHIKV and ZIKV outbreaks throughout Barranquilla, Colombia during 2014 and 2016 and explored the potential for clustering. Incidence data were fitted using multiple Bayesian Poisson models based on a suite of explanatory variables as potential risk factors and multiple options for random effects. A best fit model was used to analyse the case incidence risk for both epidemics to identify any risk factors during their epidemics. Neighbourhoods in the northern region of Barranquilla were hotspots for the outbreaks of both CHIKV and ZIKV. Additional hotspots occurred in the south-western and central regions of the CHIKV and ZIKV outbreaks, respectively. Multivariate conditional autoregressive models strongly identified higher socioeconomic strata (SES) and residing in detached houses as risk factors for ZIKV case incidences. These novel findings challenge the belief that these infections are driven by social vulnerability and merits further study both in Barranquilla and throughout the tropical and subtropical regions of the world.

Traditional flood hazard analyses often rely on univariate probability distributions; however, in many coastal catchments, flooding is the result of complex hydrodynamic interactions between multiple drivers. For example, synoptic meteorological conditions can produce considerable rainfall-runoff, while also generating wind-driven elevated sea levels. When these drivers interact in space and time, they can exacerbate flood impacts; this phenomenon is known as compound flooding. In this paper, we build a Bayesian Network based on Gaussian copulas to generate the equivalent of 500 years of daily stochastic boundary conditions for a coastal watershed in Southeast Texas. In doing so, we overcome many of the limitations of conventional univariate approaches and are able to probabilistically represent compound floods caused by riverine and coastal interactions. We calculate the resulting water levels using a 1D steady-state hydraulic model and find that flood stages in the catchment are strongly affected by backwater effects from tributary inflows and downstream water levels. By comparing with a bathtub modeling approach, we show that simplifying the multivariate dependence between flood drivers can lead to an underestimation of flood impacts, highlighting that accounting for multivariate dependence is critical for the accurate representation of flood risk in coastal catchments prone to compound events.

The constant news about the corona virus is scary. It is not possible to separate treatment for Cancer due to COVID-19. An effective treatment comparison strategy is needed. We need to have a handy tool to understand cancer progression in this unprecedented scenario. Linking different events of cancer progression is the need of the hour. It is a methodological challenge. We provide the solutions to overcome the issue with interval between two consecutive events in motivating head and neck cancer (HNC) data.

Along with the advantages associated with access to information and fast communication, screen time from increased digital media consumption has recently been associated with adverse effects on youth well-being. To get a clearer picture of its value for global youth based sustainability initiatives, this study investigates the effects of increased digital media consumption on youth's interests in ecosystem services, sustainability and science as a means for disease prevention. We achieve this, using data on 187821 adolescent students from 50 countries worldwide. Methodologically, we rely on a mixed bivariate ordered probit representation of youth's joint interest in the biosphere (ecosystem services and sustainability) and science as a means for disease prevention, which we then estimate using Markov Chain Monte Carlo (MCMC) Methods. We found that each level increase in adolescent students' reported frequency of news blogs visits and web browsing on broad science adversely affect their interests in ecosystem services, sustainability and science as a means for disease prevention. Although each level increase in youth's frequency of ecological website visits also reduces by 20% (with 95% CI [-0.36; -0.32]) their interests in the biosphere, it is found to increase however by 3% (with 95% CI [0.02; 0.05]) their interest science as a means for disease prevention. Overall, our results highlight heterogeneous effects of digital media consumption on adolescents' well-being in terms of their interests in ecosystem services, sustainability, and science as a means for disease prevention.

This pandemic of COVID-19 is tedious to control. The only lockdown is the way to stop the spread of this infection. Conventional health care is facing a real challenge to operate. Primarily the challenge is to provide health care support for COVID-19 patients with limited resources and continue the health care services like earlier.Perhaps, this challenge is the same but magnitude is different from different geographical locations around the globe. In this article, we presented a Bayesian algorithm with the Code to predict cancer death due to COVID-19. This code is possible to run at different time points and different geographical locations around the world. This code will help us to get the best strategy and shift the treatment option for cancer treatment. The model would provide physicians with an objective tool for counseling and decision making at different hotspots and small areas to implement.

Background: All Payer Claims Databases (APCD) are a rich source of health information, however, race and ethnicity (R&E) data are largely missing. Bayesian Improved Surname Geocoding (BISG) is a common R&E imputation method, yet, validation of BISG in APCDs is lacking. We used the BISG to impute missing R&E in the Oregon APCD. Methods: BISG imputed R&E for Asian Pacific Islanders (API), Blacks, Hispanics and Whites were contrasted to the gold standard (vital statistics) and sensitivity and specificity improvements were assessed. Logistic regression examined whether missing R&E was random across patient characteristics. Results: Among 85,857 individuals in the study, 32.1% (n=27,594) had missing R&E. Missing R&E was not randomly distributed. There were higher odds of missingness among males, Whites, those age 65 and older, and commercially insured individuals. Differences in the percent missing were also found by co-morbid conditions and mortality causes. Imputing the missing R&E with BISG method improved the sensitivity to identify White, Black, API, and Hispanics. Conclusions: APCDs can benefit from enhancing missing R&E with BISG imputation to perform more robust population-health level analyses and identify inequities according to R&E without losing power or dropping non-random records with missing R&E data.

The creation of computer models is often driven by the need to make predictions in regions where there is no data (i.e. extrapolations). This makes validation challenging as it is difficult to ensure that a model will be suitable when it is applied in a region where there are no observations of the system of interest. The current paper proposes a method that can reveal flaws in a model which may be difficult to identify using traditional approaches for model calibration and validation. The method specifically targets the situation where one is attempting to model a dynamical system that is believed to possess time-invariant calibration parameters. The proposed approach allows these parameters to vary with time, even though it is believed that they are time-invariant. The of such an analysis is to identify key discrepancies - indications that a model has inherent flaws and, as a result, should not be used to influence decisions in regions where there is no data. The proposed method isn't necessarily a predictor of extrapolation performance, rather, it is a stringent test that, the authors believe, should be applied before extrapolation is attempted. The approach could therefore form a useful part of wider validation frameworks in the future.

The enormous increase in the volume of Earth Observations (EOs) has provided the scientific community with unprecedented temporal, spatial, and spectral information. However, this increase in the volume of EOs has not yet resulted in proportional progress with our ability to forecast agricultural systems.This study examines the applicability of EOs obtained from Sentinel2 and Landsat8 for constraining the APSIM-Maize model parameters. We leveraged leaf area index (LAI) retrieved from Sentinel2 and Landsat8 NDVI to constrain a series of APSIM-Maize model parameters in three different Bayesian multi-criteria optimization frameworks across 13 different sites across the U.S Midwest. A time variant sensitivity analysis was performed to identify the most influential parameters driving the LAI estimates in APSIM-Maize model. Then surrogate models were develop using random samples taken from the parameter space using Latin hypercube sampling to emulate APSIM’s behavior in simulating NDVI and LAI at all sites. Site-level, global and hierarchical Bayesian optimization models were then developed using the site-level emulators to simultaneously constrain all parameters and estimate the site to site variability in crop parameters. For within sample predictions, site-level optimization showed the largest predictive uncertainty around LAI and crop yield, whereas the global optimization showed the most constraint predictions for these variables. Lowest RMSE for within sample yield prediction was found for hierarchical optimization scheme (1423 Kg ha−1) while the largest RMSE was found for site-level (1494 Kg ha−1). In out-of-sample predictions within the spatio-temporal extent of the training sites, global optimization showed lower RMSE (1627 Kg ha−1) compared to the hierarchical approach (1822 Kg ha−1) across 90 independent sites in the U.S Midwest. On comparison between these two optimization schemes across another 242 independent sites outside the spatio-temporal extent of the training sites, global optimization also showed substantially lower RMSE (1554 Kg ha−1) as compared to the hierarchical approach (2532 Kg ha−1). Overall, EOs demonstrated their real use case for constraining process-based crop models and showed comparable results to model calibration exercises using only field measurements.

This article describes techniques employed in the production of a synthetic dataset of driver telematics emulated from a similar real insurance dataset. The synthetic dataset generated has 100,000 policies that included observations about driver’s claims experience together with associated classical risk variables and telematics-related variables. This work is aimed to produce a resource that can be used to advance models to assess risks for usage-based insurance. It follows a three-stage process using machine learning algorithms. The first stage is simulating values for the number of claims as multiple binary classifications applying feedforward neural networks. The second stage is simulating values for aggregated amount of claims as regression using feedforward neural networks, with number of claims included in the set of feature variables. In the final stage, a synthetic portfolio of the space of feature variables is generated applying an extended SMOTE algorithm. The resulting dataset is evaluated by comparing the synthetic and real datasets when Poisson and gamma regression models are fitted to the respective data. Other visualization and data summarization produce remarkable similar statistics between the two datasets. We hope that researchers interested in obtaining telematics datasets to calibrate models or learning algorithms will find our work valuable.

This paper investigates the inverse option problems (IOP) in the extended Black--Scholes model arising in financial markets. We identify the volatility and the drift coefficient from the measured data in financial markets using a Bayesian inference approach, which is presented as an IOP solution. The posterior probability density function of the parameters is computed from the measured data. The statistics of the unknown parameters are estimated by a Markov Chain Monte Carlo (MCMC) algorithm, which exploits the posterior state space. The efficient sampling strategy of the MCMC algorithm enables us to solve inverse problems by the Bayesian inference technique. Our numerical results indicate that the Bayesian inference approach can simultaneously estimate the unknown trend and volatility coefficients from the measured data.

In the past few decades, probabilistic interpretations of brain functions have become widespread in cognitive science and neuroscience. In particular, the free energy principle and active inference are increasingly popular theories of cognitive functions that claim to offer a unified understanding of life and cognition within a general mathematical framework derived from information and control theory, and statistical mechanics. However, we argue that if the active inference proposal is to be taken as a general process theory for biological systems, it is necessary to understand how it relates to existing control theoretical approaches routinely used to study and explain biological systems. For example, recently, PID control has been shown to be implemented in simple molecular systems and is becoming a popular mechanistic explanation of behaviours such as chemotaxis in bacteria and amoebae, and robust adaptation in biochemical networks. In this work, we will show how PID controllers can fit a more general theory of life and cognition under the principle of (variational) free energy minimisation when using approximate linear generative models of the world. This more general interpretation provides also a new perspective on traditional problems of PID controllers such as parameter tuning as well as the need to balance performances and robustness conditions of a controller. Specifically, we then show how these problems can be understood in terms of the optimisation of the precisions (inverse variances) modulating different prediction errors in the free energy functional.

Establishing an accurate mathematical model is the foundation of simulating the motion of marine vehicles and structures, and it is the basis of modeling-based control design. System identification from observed input-output data is a practical and powerful method. However, for modeling objects with different characteristics and known information, a single modeling framework can hardly meet the requirements of model establishment. Moreover, there are some challenges in system identification, such as parameter drift and overfitting. In this work, three robust methods are proposed for generating ocean hydrodynamic models based on Bayesian regression. Two Bayesian techniques, semi-conjugate linear regression and noisy input Gaussian process regression, are used for parametric and nonparametric gray-box modeling and black-box modeling. The experimental free-running tests of the KVLCC2 ship model and a multi-freedom wave energy converter (WEC) are used to validate the proposed Bayesian models. The results demonstrate that the proposed schemes for system identification of the ship and WEC have good generalization ability and robustness. Finally, the developed modeling methods are evaluated considering the aspects required conditions, operating characteristics and prediction accuracy.

With data for SARS-CoV-2 and with many countries entering the second wave of contagion required the improvement of the forecasting model, structuring its model to forecast the peak of the first and second contagion wave in Brazil. The Model Estimation for SARS-CoV-2 Peak Contagion (MESPC) was structured, capable of estimating the peak of contagion for SARS-CoV-2 in the first and second waves, as the main objective of this work. Using the MESPC model, it was possible to estimate, with a certain reliability degree, the peak of contagion for the first and second waves in Brazil, with one day difference from the real to the forecast. It is possible to use MESPC to forecast the peak of contagion for several regions, provided that the necessary structure and calibration are respected.

Construction projects are usually operating in a complex and dynamic environment in which the accumulation of many interrelated factors causes high uncertainty. Construction projects are complex and frequently involve substantial uncertainties including process complicatedness, intricate organization structure, dynamic environment, and financial strain. The study aims to categorize the influencing factors into three groups, namely construction project system, economic-market climate, and external environment. It attempts to adopt a novel analysis tool to examine the relationship between the project cost and multiple influencing factors by using Bayesian SEM. While the Bayesian SEM method has been receiving increasing attention in exploring the relationship between latent variables, construction studies still heavily rely on the covariance-based SEM approach. This study introduces several advantages of Bayesian SEM that make it more flexible and powerful than covariance-based SEM and provides the foundation of Bayesian SEM estimation and inference by illustrating this method in a project cost application.

The eastern subterranean termite, Reticulitermes flavipes, currently inhabits previously glaciated regions of the northeastern U.S., as well as the unglaciated southern Appalachian Mountains and surrounding areas. We hypothesized that Pleistocene climatic fluctuations have influenced the distribution of R. flavipes, and thus the evolutionary history of the species. We estimated contemporary and historical geographic distributions of R. flavipes by constructing Species Distribution Models (SDM). We also inferred the evolutionary and demographic history of the species using mitochondrial (cytochrome oxidase I and II) and nuclear (endo-beta-1,4-glucanase) DNA sequence data. To do this, genetic populations were delineated using Bayesian spatial genetic clustering, competing hypotheses about population divergence were assessed using approximate Bayesian computation (ABC), and changes in population size were estimated using Bayesian skyline plots. SDMs identified areas in the north with suitable habitat during the transition from the Last Interglacial to the Last Glacial Maximum, as well as an expanding distribution from the mid-Holocene to the present. Genetic analyses identified three geographically cohesive populations, corresponding with northern, central, and southern portions of the study region. Based on ABC analyses, divergence between the Northern and Southern populations was the oldest, estimated to have occurred 64.80 thousand years ago (kya), which corresponds with the timing of available habitat in the north. The Central and Northern populations diverged in the mid-Holocene, 8.63 kya, after which the Central population continued to expand. Accordingly, phylogeographic patterns of R. flavipes in the southern Appalachians appear to have been strongly influenced by glacial-interglacial climate change.

Technology development changes society and society demands new and innovative technology development. We analyze technology to understand society and technology itself. Many researches have been introduced in various fields. Most of them were about patent analysis. This is because detailed and accurate results of research and development are patented. In this paper, we study on new patent analysis method based on count data model and Bayesian regression analysis. Using count data model, we analyze the technological keywords extracted from the collected patent documents. We use the posterior distribution of Bayesian statistics to reflect the experience and knowledge of the relevant technological experts in the analysis model. Moreover, we apply the proposed model to finding sustainable technologies. Finding and developing sustainable technologies is an important activity for companies and research institutes to maintain their technological competitiveness. To illustrate how our modeling could be applied to real domain, we carry out a case study using the patent documents related to artificial intelligence.

We discuss the influence played by different statistical models in the prediction of porosity and litho-fluid facies from logged and post-stack inverted acoustic impedance (Ip) values. We compare the inversion and classification results obtained under three different a-priori statistical assumptions: an analytical Gaussian distribution, an analytical Gaussian-mixture model and a non-parametric mixture distribution. The first model assumes Gaussian distributed porosity and Ip values, thus neglecting their facies-dependent behaviour caused by different lithologic and saturation conditions. Differently, the other two statistical models relate each component of the mixture to a specific litho-fluid facies, so that the facies-dependency of porosity and Ip values is taken into account. Blind well tests are used to validate the final predictions, whereas the analysis of the maximum-a-posteriori (MAP) solutions, the coverage ratio and the contingency analysis tools are used to quantitatively compare the inversion outcomes. This work points out that the correct choice of the statistical petrophysical model could be crucial in reservoir characterization studies. Indeed, for the investigated zone it turns out that the simple Gaussian model constitutes an oversimplified assumption, while the two mixture models provide more accurate results, although the non-parametric one yields slightly superior predictions with respect to the Gaussian-mixture assumption.

Machine learning forks into three main branches such as supervised learning, unsupervised learning, and reinforcement learning where reinforcement learning is much potential to artificial intelligence (AI) applications because it solves real problems by progressive process in which possible solutions are improved and finetuned continuously. The progressive approach, which reflects ability of adaptation, is appropriate to the real world where most events occur and change continuously and unexpectedly. Moreover, data is getting too huge for supervised learning and unsupervised learning to draw valuable knowledge from such huge data at one time. Bayesian optimization (BO) models an optimization problem as a probabilistic form called surrogate model and then directly maximizes an acquisition function created from such surrogate model in order to maximize implicitly and indirectly the target function for finding out solution of the optimization problem. A popular surrogate model is Gaussian process regression model. The process of maximizing acquisition function is based on updating posterior probability of surrogate model repeatedly, which is improved after every iteration. Taking advantages of acquisition function or utility function is also common in decision theory but the semantic meaning behind BO is that BO solves problems by progressive and adaptive approach via updating surrogate model from a small piece of data at each time, according to ideology of reinforcement learning. Undoubtedly, BO is a reinforcement learning algorithm with many potential applications and thus it is surveyed in this research with attention to its mathematical ideas. Moreover, the solution of optimization problem is important to not only applied mathematics but also AI.

Chikungunya virus (CHIKV) and Zika virus (ZIKV) have recently emerged as global infections with consequential disability adjusted life years (DALYs) and economic burden. This study aimed to explore the spatiotemporal heterogeneity in the occurrence of CHIKV and ZIKV outbreaks throughout Barranquilla, Colombia during 2014 and 2016 and explored the potential for clustering. Incidence data were fitted using multiple Bayesian Poisson models based on a suite of explanatory variables as potential risk factors and multiple options for random effects. A best fit model was used to analyse the case incidence risk for both epidemics to identify any risk factors during their epidemics. Neighbourhoods in the northern region of Barranquilla were hotspots for the outbreaks of both CHIKV and ZIKV. Additional hotspots occurred in the south-western and central regions of the CHIKV and ZIKV outbreaks, respectively. Multivariate conditional autoregressive models strongly identified higher socioeconomic strata (SES) and residing in detached houses as risk factors for ZIKV case incidences. These novel findings challenge the belief that these infections are driven by social vulnerability and merits further study both in Barranquilla and throughout the tropical and subtropical regions of the world.

Bayesian model calibration techniques are commonly employed in the characterization of nonlinear dynamic systems, as they provide a conceptual and effective framework to deal with model uncertainties, experimental errors and procedure assumptions. This understanding has resulted in the need to introduce a model discrepancy term to account for the differences between model-based predictions and real observations. Indeed, the goal of this work is to enhance model-driven Structural Health Monitoring procedures by incorporating the posterior uncertainty linked to updated model discrepancy, and thus make relevant considerations for its use in the Structural Health Monitoring. Specifically, the Bayesian inference has been applied to the calibration of nonlinear hysteretic systems to both provide: (i) most probable values (MPV) of the parameters following the calibration, and; (ii) estimates of the model discrepancy posterior distribution. The effect of the model discrepancy in the calibration is first illustrated recurring to a single degree of freedom Bouc-Wen type oscillator, and then applied for calibrating a reference nonlinear Bouc-Wen model, deriving from real data acquired on a monitored masonry building.

Understanding the spatially varying effects of demographic factors on the spatio-temporal variation of intestinal parasites infections is important for public health intervention and monitoring. This paper presents a hierarchical Bayesian spatially varying coefficient model to evaluate the effects demographic factors on intestinal parasites morbidities in Ghana. The modeling relied on morbidity data collected by the District Health Information Management Systems. We developed Poisson and Poisson-gamma spatially varying coefficient Models. We used the demographic factors, unsafe drinking water, unsafe toilet and unsafe liquid waste disposal as model covariates. The models were fitted using the Integrated Nested Laplace Approximations (INLA). The overall risk of intestinal parasites infection was estimated to be 10.9 per 100 people with a wide spatial variation in the district-specific posterior risk estimates. Substantial spatial variation of increasing multiplicative effects of unsafe drinking water, unsafe toilet and unsafe liquid waste disposal occurs on the variation of intestinal parasites risk. The structured residual spatial variation widely dominates the unstructured component, suggesting that the unaccounted-for risk factors are spatially continuous in nature. The study concludes that both the spatial distribution of the posterior risk and the associated exceedance probability maps are essential for monitoring and control of intestinal parasites.

This review covers some of the central relationships between artificial intelligence, psychology, and economics through the lens of information theory, specifically focusing on formal models of decision-theory. In doing so we look at a particular approach that each field has adopted and how information theory has informed the development of the ideas of each field. A key theme is expected utility theory, its connection to information theory, and the Bayesian approach to decision-making and forms of (bounded) rationality. What emerges from this review is a broadly unified formal perspective derived from three very different starting points that reflect the unique principles of each field. Each of the three approaches reviewed can, in principle at least, be implemented in a computational model in such a way that, with sufficient computational power, they could be compared with human abilities in complex tasks. However, a central critique that can be applied to all three approaches was first put forward by Savage in {\it The Foundations of Statistics} and recently brought to the fore by the economist Binmore: Bayesian approaches to decision-making work in what Savage called `small worlds' but cannot work in `large worlds'. This point, in various different guises, is central to some of the current debates about the power of artificial intelligence and its relationship to human-like learning and decision-making. Recent work on artificial intelligence has gone some way to bridging this gap but significant questions still need to be answered in all three fields in order to make progress on these problems.

In this study, based on our previous study, we examined the mathematical properties, especially the stability of the equilibrium for our proposed mathematical model. By means of the results of the stability in this study, we also used actual data representing transient booms and resurgent booms, and conducted parameter estimation for our proposed model using Bayesian inference. In addition, we conducted a model fitting to five actual data. By this study, we reconfirmed that we can express the resurgences or minute vibrations of actual data by means of our proposed model.

Classical methods for inverse problems are mainly based on regularization theory. In particular those which are based on optimization of a criterion with two parts: a data-model matching and a regularization term. Different choices for these two terms and great number of optimization algorithms have been proposed. When these two terms are distance or divergence measures, they can have a Bayesian Maximum A Posteriori (MAP) interpretation where these two terms correspond, respectively, to the likelihood and prior probability models.

Aiming at the problem that the performance of Adaptive Kalman filter estimation will be affected when the statistical characteristics of the process and measurement noise matrix are inaccurate and time-varying in the linear Gaussian state-space model, an algorithm of Multi-fading factor and update monitoring strategy adaptive Kalman filter based variational Bayesian is proposed. Inverse Wishart distribution is selected as the measurement noise model, the system state vector and measurement noise covariance matrix are estimated with the variational Bayesian method. The process noise covariance matrix is estimated by the maximum a posteriori principle, and the update monitoring strategy with adjustment factors is used to maintain the positive semi-definite of the updated matrix. The above optimal estimation results are introduced as time-varying parameters into the multiple fading factors to improve the estimation accuracy of the one-step state predicted covariance matrix. The application of the proposed algorithm in target tracking is simulated. The results show that compared with the current filters, the proposed filtering algorithm has better accuracy and convergence performance, and realizes the simultaneous estimation of inaccurate time-varying process and measurement noise covariance matrices.

Rapid post-earthquake damage diagnosis of bridges can guide decision-making for emergency response management and recovery. This can be facilitated using digital technologies to remove the barriers of manual post-event inspections. Prior mechanics-based Finite Element (FE) models can be used for post-event response simulation using the measured ground motions at nearby stations; however, the damage assessment outcomes would suffer from uncertainties in structural and soil material properties, input excitations, etc. For instrumented bridges, these uncertainties can be reduced by integrating sensory data with prior models through a model updating approach. This study presents a sequential Bayesian model updating technique, through which a linear/nonlinear FE model, including soil-structure interaction effects, and the foundation input motions are jointly identified from measured acceleration responses. The efficacy of the presented model updating technique is first examined through a numerical verification study. Then, seismic data recorded from the San Rogue Canyon Bridge in California are used for a real-world case study. Comparison between the free-field and the foundation input motions reveals valuable information regarding the soil-structure interaction effects at the bridge site. Moreover, the reasonable agreement between the recorded and estimated bridge responses shows the potentials of the presented model updating technique for real-world applications. The updated FE model is considered as the digital twin of the bridge and can be used to analyze the bridge and monitor the structural response at element, section, and fiber levels to diagnose the location and severity of any potential damage mechanism.

Sleep deficiency impacts the quality of life and may have serious health consequences in the long run. Questionnaire-based subjective assessment of sleep deficiency has many limitations. On the other hand, objective assessment of sleep deficiency is challenging. In this study, we propose a polysomnography-based mathematical model for computing baseline sleep deficiency severity score and then investigated the estimation of sleep deficiency severity using features available only from wearable sensor data including heart rate variability and single-channel electroen-cephalography for a dataset of 500 subjects. We used Monte-Carlo Feature Selection (MCFS) and inter-dependency discovery for selecting the best features and removing multi-collinearity. For developing the Regression model we investigated both the frequentist and the Bayesian ap-proaches. An Artificial Neural Network achieved the best performance of RMSE = 5.47 and an R-squared value of 0.67 for sleep deficiency severity estimation. The developed method is com-parable to conventional methods of Functional Outcome of Sleep Questionnaire and Epworth Sleepiness Scale for assessing the impact of sleep apnea on sleep deficiency. Moreover, the results pave the way for reliable and interpretable sleep deficiency severity estimation using a wearable device in Smart Health.

We investigate the power of time series analysis based on a variety of information-theoretic approaches from statistics (AIC, BIC) and machine learning (Minimum Message Length) - and we then compare their efficacy with traditional time series model and with hybrids involving deep learning. More specifically, we develop AIC, BIC and Minimum Message Length (MML) ARMA (autoregressive moving average) time series models - with this Bayesian information-theoretic MML ARMA modelling already being new work. We then study deep learning based algorithms in time series forecasting, using Long Short Term Memory (LSTM), and we then combine this with the ARMA modelling to produce a hybrid ARMA-LSTM prediction. Part of the purpose of the use of LSTM is to seek capture any hidden information in the residuals left from the traditional ARMA model. We show that MML not only outperforms earlier statistical approaches to ARMA modelling, but we further show that the hybrid MML ARMA-LSTM models outperform both ARMA models and LSTM models.

While Europe was beginning to deal with the resurgence of COVID-19 due to the Delta variant, the European football championship took place, June 11 - July 11, 2021. We studied the inversion in the decrease/increase rate of new SARS-COV-2 infections in the countries of the tournament, investigating the hypothesis of an association. Using a Bayesian piecewise regression with a Poisson Generalized Linear Model, we looked for a changepoint in the timeseries of the new SARS-COV-2 cases of each country, expecting it to appear not later than two to three weeks after the date of their first match. The two slopes, before and after the changepoint, were used to discuss the reversal from a decreasing to an increasing rate of the infections. For 17 out of 22 countries (77%) the changepoint came on average 14.97 days after their first match [95% CI 12.29 to 17.47]. For all those 17 countries, the changepoint coincides with an inversion from a decreasing to an increasing rate of the infections. Before the changepoint, the new cases were decreasing, halving on average every 18.07 days [95% CI 11.81 to 29.42]. After the changepoint, the cases begin to increase, doubling every 29.10 days [95% CI 14.12 to 49.78]. This inversion in the SARS-COV-2 case rate, happened during the tournament, provides evidence in favor of a relationship

We discuss what we believe could be an improvement in future discussions of the ever-changing brain. We do so by distinguishing different types of brain variability and outlining methods suitable to analyse them. We argue that, when studying brain and behaviour data, classical methods such as regression analysis and more advanced approaches both aim to decompose the total variance into sensible variance components. In parallel, we argue that a distinction needs to be made between innate and acquired brain variability. For varying high-dimensional brain data, we present methods useful to extract their low-dimensional representations. Finally, to trace potential causes and predict plausible consequences of brain variability, we discuss how to combine statistical principles and neurobiological insights to make associative, explanatory, predictive, and causal enquires; but cautions are needed to raise association- or prediction-based neurobiological findings to causal claims.

Industrial dynamical systems often exhibit multi-scale response due to material heterogeneities, operation conditions and complex environmental loadings. In such problems, it is the case that the smallest length-scale of the systems dynamics controls the numerical resolution required to effectively resolve the embedded physics. In practice however, high numerical resolutions is only required in a confined region of the system where fast dynamics or localized material variability are exhibited, whereas a coarser discretization can be sufficient in the rest majority of the system. To this end, a unified computational scheme with uniform spatio-temporal resolutions for uncertainty quantification can be very computationally demanding. Partitioning the complex dynamical system into smaller easier-to-solve problems based of the localized dynamics and material variability can reduce the overall computational cost. However, identifying the region of interest for high-resolution and intensive uncertainty quantification can be a problem dependent. The region of interest can be specified based on the localization features of the solution, user interest, and correlation length of the random material properties. For problems where a region of interest is not evident, Bayesian inference can provide a feasible solution. In this work, we employ a Bayesian framework to update our prior knowledge on the localized region of interest using measurements and system response. To address the computational cost of the Bayesian inference, we construct a Gaussian process surrogate for the forward model. Once, the localized region of interest is identified, we use polynomial chaos expansion to propagate the localization uncertainty. We demonstrate our framework through numerical experiments on a three-dimensional elastodynamic problem

Fish schooling implies an awareness of the swimmers for their companions. In flow mediated environments, in addition to visual cues, pressure and shear sensors on the fish body are critical for providing quantitative information that assists the quantification of proximity to other swimmers. Here we examine the distribution of sensors on the surface of an artificial swimmer so that it can optimally identify a leading group of swimmers. We employ Bayesian experimental design coupled with two-dimensional Navier Stokes equations for multiple self-propelled swimmers. The follower tracks the school using information from its own surface pressure and shear stress. We demonstrate that the optimal sensor distribution of the follower is qualitatively similar to the distribution of neuromasts on fish. Our results show that it is possible to identify accurately the center of mass and even the number of the leading swimmers using surface only information.

The average age of in-service bridges has increased in recent years in the United States. To address this issue, structural health monitoring and damage identification approaches can be employed to prioritize maintenance/replacement of aging bridges. Among the damage identification and operational health monitoring approaches, finite element (FE) model updating methods can offer a solution to evaluate the mechanics-based characteristics of bridges. However, in a real-world setting, unidentifiability and mutual dependency between model parameters, modeling errors, especially due to boundary conditions, as well as ill-conditioning of updating algorithms can pose challenges to the application of FE model updating methods. To address these challenges, this study presents a two-step FE model updating approach. In the first step, modal-based model updating is used to estimate linear model parameters mainly related to the stiffness of boundary conditions and material properties. In the second step, in order to refine parameter estimation accounting for nonlinear response behavior of the bridge, a time-domain model updating is carried out. In this step, boundary conditions are fixed at their final estimates using modal-based model updating. To prevent the convergence of updating algorithm to local solutions, the initial estimates for nonlinear material properties are selected based on their corresponding final estimates in the modal-based model updating. To validate the applicability of the two-step FE model updating approach, a series of forced-vibration experiments are designed and carried out on a pair of decommissioned and deteriorated prestressed bridge I-girders. After carrying out the two-step FE model updating, the final estimates of concrete compressive strength are shown to provide reasonable assessment of the damage extent in the girders.

There is a growing public health need for effective preventive interventions against dengue, and a safe, effective and affordable dengue vaccine against the four serotypes would be a significant achievement for disease prevention and control. Two tetravalent dengue vaccines, Dengvaxia (Sanofi Pasteur) and DENVax (Takeda Pharmaceutical Company), have now completed phase 3 clinical trials. While Dengvaxia resulted in serious adverse events and is restricted to individuals with prior dengue infections, DENVax has shown, at first glance, some encouraging results. Using the available data for the TAK 003 trial, we estimate, via the Bayesian approach, vaccine efficacy (VE) of the post-vaccination surveillance periods. Although better measurement over long time was expected for the second part of the post-vaccination surveillance, variation in serotype-specific efficacy needs careful consideration. Besides observing that individual serostatus prior to vaccination is determinant of DENVax vaccine efficacy, we also compare the VE estimations for 12 and 18 months and we observe that the efficacy is decreasing over time. The comparison of efficacies over time is informative and very important, bring up the discussion of the role of temporary cross-immunity in dengue vaccine trials and the impact of serostatus prior to vaccination in the context of dengue fever epidemiology.

Bootstrap resampling techinques, introduced by Efron and Rubin, can be presented in a general Bayesian framework, approximating the statistical distribution of a statistical functional φ(F), where F is a random distribution function. Efron’s and Rubin’s bootstrap procedures can be extended introducing an informative prior through the Proper Bayesian bootstrap. In this paper different bootstrap techniques are used and compared in predictive classification and regression models based on ensemble approaches, i.e. bagging models involving decision trees. Proper Bayesian bootstrap, proposed by Muliere and Secchi, is used to sample the posterior distribution over trees, introducing prior distributions on the covariates and the target variable. The results obtained are compared with respect to other competitive procedures employing different bootstrap techniques. The empirical analysis reports the results obtained on simulated and real data.

Lysine is the first limiting essential amino acid in rice because it is present in the lowest quantity compared to all the other amino acids. Amino acids are the building block of proteins and play an essential role in maintaining the human body’s healthy functioning. Rice is a staple food for large proportion of the global population, thus increasing the lysine content in rice will improve its nutritional value. In this paper, we studied the lysine biosynthesis pathway in rice (Oryza Sativa) to identify the regulators of the lysine reporter gene LYSA (LOC_Os02g24354). Genetically intervening at the regulators has the potential to increase the overall lysine content in rice. We modeled the lysine biosynthesis pathway in rice seedlings under normal and saline (NaCl) stress conditions using Bayesian networks. We estimated the model parameters using experimental data and identified the gene DAPF(LOC_Os12g37960) as a positive regulator of the lysine reporter gene LYSA under both normal and saline stress conditions. Based on this analysis, we conclude that the gene DAPF is a potent candidate for genetic intervention. Upregulating DAPF using methods such as CRISPR-Cas9 has the potential to upregulate the lysine reporter gene LYSA and increase the overall lysine content in rice.

Tephra plumes can cause a significant hazard for surrounding towns, infrastructure, and air traffic. The current work presents the use of a small and compact X-band Multi-Parameter (X-MP) radar for the remote tephra detection and tracking of two eruptive events at Merapi Volcano, Indonesia, in May and June 2018. Tephra detection was done by analysing the multiple parameters of radar: copolar correlation and reflectivity intensity. These parameters were used to cancel unwanted clutter and retrieve tephra properties, which are grain size and concentration. Real-time spatial and temporal forecasting of tephra dispersal was performed by applying an advection scheme (nowcasting) in the manner of Ensemble Prediction System (EPS). Cross-validation was done using field-survey data, radar observations, and Himawari-8 imagery. The nowcasting model computed both the displacement and growth and decaying rate of the plume based on the temporal changes in two-dimensional movement and tephra concentration, respectively. Our results with ground-based data, where the radar-based estimated grain size distribution fell within the range of in-situ data. The uncertainty of real-time forecasted tephra plume depends on the initial condition, which affects the growth-and decaying rate estimation. The EPS improves the predictability rate by reducing the number of missed and false forecasted events. Our findings and the method presented here are suitable for early warning of tephra fall hazard at the local scale.

The ISO 2859 and ISO 3951 series provide acceptance sampling procedures for lot inspection, allowing both sample size and acceptance rule to be determined, starting from a specific value either for the consumer or producer risk. However, insufficient resources often make it difficult to implement “ISO sampling plans.” In cases where the sample size is determined by external constraints, the focus shifts from determining sample size to determining consumer and producer risks. Moreover, if the sample size is very low (e.g. one single item), prior information should be included in the statistical analysis. For this reason, it makes sense to work within a Bayesian theoretical framework, such as that described in JCGM 106. Accordingly, the approach from JCGM 106 is adopted and broadened so as to allow application to lot inspection. The discussion is based on a “real-life” example of lot inspection on the basis of a single item. Starting from simple assumptions, expressions for both the prior and posterior distributions are worked out, and it is shown how the concepts from JCGM 106 can be reinterpreted in the context of lot inspection. Conceptual differences regarding the definition of consumer and producer risks in JCGM 106 and in the ISO acceptance sampling standards are elucidated and a numerical example is provided.

Next-generation cancer and oncology research needs to take full advantage of the multi-modal structured, or graph, information, with the graph datatypes ranging from molecular structures to spatially resolved imaging and digital pathology to biological networks to knowledge graphs. Graph Neural Networks (GNNs) efficiently combine the graph structure representations with the high predictive performance of deep learning, especially on the large multi-modal datasets. In this review article, we survey the landscape of recent (2020-present) GNN applications in the context of cancer and oncology research, and delineate six currently predominant research areas. Subsequently, we identify the most promising directions for future research. We compare GNNs with graphical models and "non-structured" deep learning, and devise the guidelines for cancer and oncology researchers or physician-scientists asking the question of whether they should adopt the GNN methodology in their research pipelines.

In this paper we perform numerous numerical studies for the problem of low-rank matrix completion. We compare the Bayesian approaches and a recently introduced de-biased estimator which provides a useful way to build confidence intervals of interest. From a theoretical viewpoint, the de-biased estimator comes with a sharp minimax-optimal rate of estimation error whereas the Bayesian approach reaches this rate with an additional logarithmic factor. Our simulation studies show originally interesting results that the de-biased estimator is just as good as the Bayesian estimators. Moreover, Bayesian approaches are much more stable and can outperform the de-biased estimator in the case of small samples. However, we also find that the length of the confidence intervals revealed by the de-biased estimator for an entry is absolutely shorter than the length of the considered credible interval. These suggest further theoretical studies on the estimation error and the concentration for Bayesian methods as they are being quite limited up to present.

Often, when solving forward, inverse or data assimilation problems, only a part of the solution is needed. As a model, we consider the stationary diffusion problem. We demonstrate an algorithm that can compute only a part or a functional of the solution, without calculating the full inversion operator and the complete solution. It is a well-known fact about partial differential equations that the solution at each discretisation point depends on the solutions at all other discretisation points. Therefore, it is impossible to compute the solution only at one point, without calculating the solution at all other points. The standard numerical methods like a conjugate gradient or Gauss elimination compute the whole solution and/or the complete inverse operator. We suggest a method which can compute the solution of the given partial differential equation 1) at a point; 2) at few points; 3) on an interface; or a functional of the solution, without computing the solution at all points. The required storage cost and computational resources will be lower as in the standard approach. With this new method, we can speed up, for instance, computation of the innovation in filtering or the likelihood distribution, which measures the data misfit (mismatch). Further, we can speed up the solution of the regression, Bayesian inversion, data assimilation, and Kalman filter update problems. Applying additionally the hierarchical matrix approximation, we reduce the cubic computational cost to almost linear $\mathcal{O}(k^2n \log^2 n)$, where $k\ll n$ and $n$ is the number of degrees of freedom. Up to the hierarchical matrix approximation error, the computed solution is exact. One of the disadvantages of this method is the need to modify the existing deterministic solver.

Limited accessibility in some parts of Indonesia causes difficulties in logistic distribution, especially for emergency and medical supplies. Emerging drone cargo technology is a potential solution to improve logistic distribution in those areas. However, implementing drone cargo technology involves unknown risks, both in technical and non-technical aspects. Since data on drone operations in Indonesia is limited, a new method is explored to build a Bayesian Network (BN) model for risk analysis of drone crashes in Indonesia’s outermost and underdeveloped areas. The method optimizes the modelling process, in which significant risk factors are selected based on three drone operator companies’ experiences, which include wind speed, rain intensity, and system component failure or malfunction. Real wind speed and rain probability data are then implemented in the model. The operator’s data shows that wind speed contributes to drone crashes, which can be appropriately modelled in the BN model. The model produced a probability of safe operation of 94.1%, comparable to the annual operator’s data. The result shows that most operations are safe, with a minimum case of crashing and no case of harming human life.

Impulsive noise is the main limiting factor for transmission over channels affected by electromagnetic interference. We study the estimation of (correlated) Gaussian signals in an impulsive noise scenarios. In this work, we analyze some of the existing as well as some novel estimation algorithms. Their performance is compared, for the first time, for different channel conditions, including the Markov-Middleton scenario, where the impulsive noise switches between different noise states. Following a modern approach in digital communications, the receiver design is based on a factor graph model and implements a message passing algorithm. The correlation among signal samples as well as among noise states brings about a loopy factor graph, where an iterative message passing scheme should be employed. As it is well known, approximate variational inference techniques are necessary in these cases. We propose and analyze different algorithms and provide a complete performance comparison among them, showing that both Expectation Propagation, Transparent Propagation, and the Parallel Iterative Schedule approaches reach a performance close to the optimal, at different channel conditions.

The fatigue life of concrete is affected by many interwoven factors whose effect is nonlinear. Be-cause of its unique self-learning ability and strong generalization capability, the Bayesian regu-larized backpropagation neural network (BR-BPNN) is proposed to predict concrete behavior in tensile fatigue. The optimal model was determined through various combinations of network parameters. The average relative impact value (ARIV) was constructed to evaluate the correla-tion between fatigue life and its influencing parameters (maximum stress level Smax, stress ratio R, static strength f, failure probability P). ARIV results were also compared with other factor as-sessment methods (weight equation and multiple linear regression analyses). Using BR-BPNN, S-N curves were then obtained for the combinations of R=0.1, 0.2, 0.5; f=5, 6, 7MPa; P=5%, 50%, 95%. The tensile fatigue results under different testing conditions were finally compared for compatibility. It was concluded that Smax has the most significant negative effect on fatigue life; the degree of influence of R, P, and f, which positively correlate with fatigue life, decreases suc-cessively. ARIV is confirmed as a feasible way to analyze the importance of parameters and could be recommended for future applications. The tensile fatigue performance of plain concrete under different stress states (flexural tension, axial tension, splitting tension) does not differ sig-nificantly. Besides utilizing the valuable fatigue test data scattered in the literature, insights gained from this work could provide a reference for subsequent fatigue test program design and fatigue evaluation.

Algorithms for machine learning have found extensive use in numerous fields and applications. One important aspect of effectively utilizing these algorithms is tuning the hyperparameters to match the specific task at hand. The selection and configuration of hyperparameters directly impact the performance of machine learning models. Achieving optimal hyperparameter settings often requires a deep understanding of the underlying models and the appropriate optimization techniques. While there are many automatic optimization techniques available, each with its own advantages and disadvantages, this article focuses on hyperparameter optimization for well-known machine learning models. It explores cutting-edge optimization methods and provides guidance on applying them to different machine learning algorithms. The article also presents real-world applications of hyperparameter optimization by conducting tests on spatial data collections for landslide susceptibility mapping. Based on the experiment's results, both Bayesian optimization and metaheuristic algorithms showed promising performance compared to baseline algorithms. For example, the metaheuristic algorithm improved the overall accuracy of the random forest model. Additionally, Bayesian algorithms, such as Gaussian processes, performed well for models like KNN and SVM. The paper thoroughly discusses the reasons behind the efficiency of these algorithms. By successfully identifying appropriate hyperparameter configurations, this research paper aims to assist researchers, spatial data analysts, and industrial users in developing machine learning models more effectively. The findings and insights provided in this paper can contribute to enhancing the performance and applicability of machine learning algorithms in various domains.

Chronic non-cancer pain is a highly debilitating condition affecting approximately 20% of the global population. Chronic pain may lead to significant physical disability, emotional distress, social isolation and financial burden. Whilst^. pharmacological therapies remain the cornerstone of pain management in non-cancerous chronic pain, factors including the current opioid epidemic have led to non-pharmacological techniques becoming a more attractive proposition. We explored the prevalence of medical device use and their treatment efficacy in non-cancer pain management. A systematic methodology was developed, peer reviewed and published in PROSPERO (CRD42021235384). Key words of medical device, pain management devices, chronic pain, lower back pain, back pain, leg pain and chronic pelvic pain using Science direct, PubMed, Web of Science, PROSPERO, MEDLINE, EMBASE, PorQuest and ClinicalTrials.gov. All clinical trials, epidemiology and mixed methods studies that reported the use of medical devices for non-cancer chronic pain management published between the 1^st of January 1990 and the 30^th of April 2022 were included. 13 studies were included in systematic review, of these 6 were used in the meta-analysis with 173 participants. Our meta-analysis for pain reduction in each study showed that transcutaneous electrical nerve stimulation combined with instrument-assisted soft tissue mobilization treatment and pulsed electromagnetic therapy produced significant treatment on chronic lower back pain patients. Pooled evidence revealed the use of medical device related interventions resulted in 0.7 degree of pain reduction under a 0-10 scale. Significant improvement in disability scores, with a 7.44 degree reduction in disability level compared to a placebo using a 50 score range was also seen. The application of medical devices in patients with chronic pain has gained popularity due to increasingly cost effective techniques, minimally invasiveness and greater awareness of risks associated with pharmacological management. Our analysis has shown that the optimal use of medical devices in a sustainable manner requires further extensive research, needing larger cohort studies with greater gender parity, in a more diverse range of geographical locations.

The tobacco epidemic is one of the leading public health threats the world has ever faced and public health policy that seeks to limit the problem may not only have to target the price of tobacco but also the initiation stage in a smoker’s life – the adolescent stage. This research contributes to the health economics literature by using a Bayesian hierarchical logistic model, estimated using Hamiltonian Monte Carlo (HMC) methods to empirically identify what drives the intentions to quit smoking among adolescent smokers in Zambia. Results suggest that among the junior secondary school-going adolescent smokers in Zambia, about 63% have plans to quit smoking. We find socio-demographic characteristics and several tobacco-smoking-related factors as the main drivers of adolescent smokers’ plans to quit smoking. Most importantly, we provide insights that could be useful to help adolescent smokers realize their quitting plans.

In many Bayesian computations, first, we obtain the expression of the joint distribution of all the unknown variables given the observed data. In general, this expression is not separable in those variables. Thus, obtaining their marginals for each variable and computing the expectations are difficult and costly. This problem becomes even more difficult in high dimensional quandaries, which is an important issue in inverse problems. We may then try to propose a surrogate expression with which we can do approximate computations. Often a separable expression approximation can be useful enough. The Variational Bayesian Approximation (VBA) is a technique that approximates the joint distribution $p$ with an easier, for example separable, one $q$ by minimizing Kullback–Leibler Divergence $KL(q|p)$. When $q$ is separable in all the variables, the approximation is also called Mean Field Approximation (MFA) and so $q$ is the product of the approximated marginals. A first standard and general algorithm is alternate optimization of $KL(q|p)$ with respect to $q_i$. A second general approach is its optimization in the Riemannian manifold. However, in this paper, for practical reasons, we consider the case where $p$ is in the exponential family and so is $q$. For this case, $KL(q|p)$ becomes a function of the parameters $\thetab$ of the exponential family. Then, we can use any other optimization algorithm to obtain those parameters. In this paper, we compare three optimization algorithms: standard alternate optimization, a gradient-based algorithm and a natural gradient algorithm and study their relative performances on three examples.

Bayesian Maximum Entropy (BME) is increasingly used in predicting and mapping spatio-temporal data. However, studies that have fully evaluated its robustness empirically are rare. Therefore, this research examined empirically the effect of skewness, sample size and spatial dependency level using simulated data. We considered symmetric data, data positively skewed by 1, 3, 6 and 9, data with weak, moderate, and strong spatial dependency levels, and sample sizes from 100 to 500 at the interval length of 50. The results showed that the variation of sample sizes and spatial dependency levels do not affect the Mean Square Error (MSE) and bias of BME prediction. However, skewness affects the MSE of prediction but does not affect the bias. This result indicates that BME is robust to sample size and is unbiased. Despite the significant difference due to skewness, a graphical plot showed values of MSE close to zero, suggesting that BME can be considered robust to skewness.

The main goal of this paper is to show that Bayesian optimization could be regarded as a general framework for the data driven modelling and solution of problems arising in water distribution systems. Hydraulic simulation, both scenario based, and Monte Carlo is a key tool in modelling in water distribution systems. The related optimization problems fall in a simulation/optimization framework in which objectives and constraints are often black-box. Bayesian Optimization (BO) is characterized by a surrogate model, usually a Gaussian process, but also a random forest and increasingly neural networks and an acquisition function which drives the search for new evaluation points. These modelling options make BO nonparametric, robust, flexible and sample efficient particularly suitable for simulation/optimization problems. A defining characteristic of BO is its versatility and flexibility, given for instance by different probabilistic models, in particular different kernels, different acquisition functions. These characteristics of the Bayesian optimization approach are exemplified by the two problems: cost/energy optimization in pump scheduling and optimal sensor placement for early detection on contaminant intrusion. Different surrogate models have been used both in explicit and implicit control schemes. Showing that BO can drive the process of learning control rules directly from operational data. BO can also be extended to multi-objective optimization. Two algorithms have been proposed for multi-objective detection problem using two different acquisition functions.

The description of stellar interiors remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter $\rho _{0} = 2.8\times 10^{14} {\rm\ g\ cm^{-3}}$, regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to five to twenty times $\rho_{0}$, little can be said about the microphysics of such obejects. Here, we employ an MCMC strategy in order to access the variability of polytropic three-pircewised models for neutron star equation of states. With a fixed description of the hadronic matter we explore a variety of models for the high density regimes leading to stellar masses up to 2.5 $M_{\odot}$. In addition, we also discuss the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from LIGO was used as inputs and treated as data set for testing case.

This study probes the probabilistic features of major fire hazards in enclosed spaces to establish their importance to the occurrence of fires onboard ro-ro passenger ships and in turn to raise effective operational countermeasures. Distinct from the previous studies, the present research employs Bayesian Network (BN) analysis to determine the probabilities of fire hazards more effectively. The findings of the research include five critical basic events (BE) identified namely, vehicle engine fire (fuel system fault), vehicle electric fire (electrical equipment defect or short circuit), used car electrical fire, reefer units electrical fire (electrical appliances defects or short circuit), and cargo spontaneous combustion. Additionally, the risk of fire for lithium - ion battery powered vehicles is also highlighted in the process of BN analysis, which prompts the authors to propose preventive measures for mitigating the possibility of fire occurrence on this type of electric vehicles. It is hoped that these measures can be essential justifications for establishing relevant rules regarding carrying LIB vehicles in enclosed spaces on international level.

Oral squamous cell carcinoma often arises from an oral potentially malignant disorder called oral leukoplakia (OL). With this work we aimed to develop a novel data-driven predictive model based on gene expression profiles to distinguish OL patients who underwent malignant transformation from those who did not. We used the Tree Augmented Naïve (TAN) Bayes classifier to predict the posterior probability of having oral cancer given the data. 86 patients were included with a median follow-up of 7.11 years. Fifty-one patients (51/86; 59%) underwent malignant transformation. We found that 16 genes were predictors of oral cancer in patients with OL and these included SLC7A11, SPINK6, SERPINA12, VIT, ATP1B3, CST6, FLRT2, ELMOD1, AZGP1, RNASE13, DIO2, ECM1, CYP4F11, SYTL4, AKR1C1, and AKR1C3. In conclusion, we showed that Bayesian gene networks are a data-driven approach which could be used also in other predictor models in oncology.

The persistence of geographic inequities in vaccination coverage often evidences the presence of zero-dose and missed communities and vulnerabilities to vaccine-preventable diseases. These inequities were exacerbated in many places during the COVID-19 pandemic, due to severe disruptions to vaccination services. Understanding changes in zero-dose prevalence and associated risk factors in the context of the COVID-19 pandemic is, therefore, critical to designing effective strategies to reach vulnerable populations. Using data from nationally representative household surveys conducted before the COVID-19 pandemic in 2018 and during the pandemic in 2021 in Nigeria, we fitted Bayesian geostatistical models to map the distribution of three vaccination coverage indicators: receipt of DTP1, MCV1 and any of the four basic vaccines (BCG, OPV0, DTP1 and MCV1), and corresponding zero-dose estimates independently at 1x1 km resolution and the district level during both time periods. We also explored changes in the factors associated with non-vaccination at the national and regional levels using multilevel logistic regression models. Our results revealed no increases in zero-dose prevalence due to the pandemic at the national level, although considerable increases were observed in a few districts. We found substantial subnational heterogeneities in vaccination coverage and zero-dose prevalence both before and during the pandemic, showing broadly similar patterns in both time periods. Areas with relatively higher zero-dose prevalence occurred mostly in the north and a few places in the south in both time periods. We also found consistent areas of low coverage and high zero-dose prevalence using all three zero-dose indicators, revealing areas in greatest need. At the national level, risk factors related to socioeconomic/ demographic status (e.g., maternal education), maternal access to and utilization of health services and remoteness were strongly associated with the odds of zero dose in both time periods, while those related to communication were mostly relevant before the pandemic. These associations were also supported at the regional level, but we additionally identified risk factors specific to zero dose in each region; for example, communication and cross-border migration in the northwest. Our findings can help guide tailored strategies to reduce zero-dose prevalence and boost coverage levels in Nigeria.