Sparse direction-of-arrival (DOA) estimation of wideband signal has attracted widespread researches for its unique high-resolution performance. Numerous existing methods based on sparse Bayesian learning (SBL) don’t possess the ability to enhance sparsity even if they have already enjoyed high favor among sparse recovery approaches. In view of this, we propose a novel hierarchical Bayesian prior framework that enhance sparsity evidently and derive its corresponding iterative algorithm. On analysis, the computational complexity of the iterative are less than most existing other state-of-the-art algorithms. Not only that, proposed method possesses high angular estimation precision and sparsity performance by utilizing joint sparsity of multiple measure vector (MMV) models. Last but not least, the method obtains the ability to stabilize the estimated values between different frequencies or snapshots such that flat spatial spectrum. Extensive simulation results are presented to prove the superior performance of our methods

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.

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.

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.

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.

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 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.

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.

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.

Most studies modelling population mobility and the spread of infectious diseases, particularly using meta-population-multi-patched models, tend to focus on theoretical properties and numerical simulations of such models. There is relatively scanty literature published on fit, inference and uncertainty quantification on epidemic models with population mobility. In this research, we have used three estimation techniques to solve an inverse problem and quantify its uncertainty on a human mobility-based multi-patched epidemic model, using mobile phone sensing data and COVID-19 confirmed positive cases in Hermosillo, Mexico. First, we have utilized a Brownian bridge model using mobile phone GPS data to estimate residence and mobility parameters of the epidemic model. In the second step, we have estimated the optimal model epidemiological parameters by deterministically inverting the model using a Darwinian inspired evolutionary algorithm (EA) known as the genetic algorithm (GA). The third part of the analysis involves performing inference and uncertainty quantification on the epidemic model using two Bayesian Monte Carlo sampling methods: t-walk and Hamiltonian Monte Carlo (HMC). The results show that the estimated model parameters and incidence adequately fit the observed daily COVID-19 incidence in Hermosillo. Moreover, the estimated parameters from HMC result into large credible intervals, improving their coverage for the observed and predicted daily incidences. We also observe improved predictions when using multi-patch model with mobility against the single-patch model.

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 present work presents and discusses experimental performance and gaseous emission results of the use of ternary mixtures of diesel oil, biodiesel and ethanol in an ASTM-CFR engine. Eight different mixtures with biodiesel content equal to 8%, 20% and 60%, and ethanol content equal to 0%, 5% and 8% were analyzed. After the selection and preparation of the mixtures, the determination of the cetane number for each one was performed following the ASTM D613 standard. After cetane number evaluation, two off-standard operating conditions were selected, and performance and emissions were obtained. A thermodynamic analysis and a Bayesian inverse approach were applied and physical quantities of interest, such as indicated work, maximum pressure, rate of heat released by the fuel, ignition delay, and combustion duration, were obtained. In-cylinder pressure data were obtained through an AVL Indicating system installed in the ASTM-CFR engine. The heating values of each mixture were determined using a bomb calorimeter and thermal and indicated efficiencies were evaluated based on the total heat released by the fuel and the indicated work. Finally, a gaseous emissions analysis was conducted to evaluate the influence of the blend compositions, as well as of their cetane numbers, on the formation of NOx, CO and CO₂.

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.

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.

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.

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.

The Ocean Color - Simultaneous Marine and Aerosol Retrieval Tool (OC-SMART) is a robust data processing platform that supports a large array of multi-spectral and hyper-spectral sensors. It provides accurate aerosol optical depths and remote sensing reflectances (Rrs estimates) that can be used to generate products such as absorption coefficients due to phytoplankton and detritus/Gelbstoff as well as backscattering coefficients due to particulate matter. The OC-SMART platform yields improved performance in complex environments by utilizing scientific machine learning (SciML) in conjunction with comprehensive radiative transfer computations. This paper expands the capability of OC-SMART by quantifying uncertainties in ocean color retrievals. Bayesian inversion is used to relate measured top of atmosphere radiances and a priori data to estimate posterior probability density functions and associated uncertainties. A framework of the methodology and implementation strategy is presented and uncertainty estimates for Rrs retrievals are provided to demonstrate the approach by applying it to MODIS, OLCI Sentinel-3, and VIIRS sensor data.

We explore, for the first time in the literature, how the revenues of the ship-management companies respond to macroeconomic exogenous shocks. Using data for ship-management companies in Cyprus, we find evidence that a demand shock has the largest impact on revenues, exhibiting an almost one-for-one relationship. If the demand shock is permanent, we observe a ceteris paribus permanent effect on revenues. Similarly, this occurs irrespective of the final effect that demand has on the relevant freight rate, proxied via the Baltic dry and tanker (dirty and clean) indices. The BDI and the BDTI indices have smaller effect on revenues, standing at approximately 0.05% for every 1% shock, while the clean tanker index does not have an effect, most likely due to their fleet composition. In accordance with the literature, we find that a shock in the price of Brent oil increases revenues. Our results bear importance not only for ship-management companies per se, but also for countries that are ship-management hubs.

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.

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 propose a novel adaptive importance sampling scheme for Bayesian inversion problems where the inference of the variables of interest and the power of the data noise is split. More specifically, we consider a Bayesian analysis for the variables of interest (i.e., the parameters of the model to invert), whereas we employ a maximum likelihood approach for the estimation of the noise power. The whole technique is implemented by means of an iterative procedure, alternating sampling and optimization steps. Moreover, the noise power is also used as a tempered parameter for the posterior distribution of the the variables of interest. Therefore, a sequence of tempered posterior densities is generated, where the tempered parameter is automatically selected according to the actual estimation of the noise power. A complete Bayesian study over the model parameters and the scale parameter can be also performed. Numerical experiments show the benefits of the proposed approach.

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.

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.

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 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.

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.

Background: Cardiovascular disease (CVD) has been one of the leading causes of death and disability-adjusted life years lost worldwide. Blood pressure, lipid, and cholesterol are good predictors of CVD risk and correspond upon age and physical fitness. However, few studies have explored the variation trend of CVD risk factors across different populations upon age and their muscle strength. Objective: to analysis the variation tendency of CVD risk factors in blood according to age and relative grip strength among different populations. Method: 25363 participants were recruited in this cross-sectional study and 24709 were included in the analysis. A logistic regression and a Bayesian probabilistic analysis based on Markov Chain Monte Carlo (MCMC) Modeling is conducted to build probability prediction models of hypertension, hyperlipidemia, and hypercholesterolemia according to age, relative grip strength, body weight conditions, and physical activity levels. Results: 1) age might be the main influence factor of hypertension, which is regarded as one of the primary CVD risk factors. However, although keeping a high level of physical activity might have positive effect on preventing hypertension because that individuals with normal body weight and higher physical activity shows a lower probability of being diagnosed with hypertension, it might could not prevent individuals from getting hypertension with age. 2) After 60, individuals of normal body weight seem more likely to have hyperlipidemia than those are overweight or obese. 3) Larger relative grip strength might not be able to offset the negative effects of obesity, overweight and physical inactivity on hyperlipidemia. 4) The probability of getting hypercholesterolemia varies less with age and relative grip strength. Conclusion: Body weight management and keeping high levels of physical activity are recommended at any age. It might benefit to increase some bodyweight after 60 years old.

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.

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 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.

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.

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.

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.

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.

When there is uncertainty in the value of parameters of the input random components of a stochastic simulation model, two-level nested simulation algorithms are used to estimate the expectation of performance variables of interest. In the outer level of the algorithm (n) observations are generated for the parameters, and in the inner level (m) observations of the simulation model are generated with the value of parameters fixed at the value generated in the outer level. In this article, we consider the case in which the observations at both levels of the algorithm are independent, showing how the variance of the observations can be decomposed into the sum of a parametric variance and a stochastic variance. Next, we derive central limit theorems that allow us to compute asymptotic confidence intervals to assess the accuracy of the simulation-based estimators for the point forecast and the variance components. Under this framework, we derive analytical expressions for the point forecast and the variance components of a Bayesian model to forecast sporadic demand; and we use these expressions to illustrate the validity of our theoretical results by performing simulation experiments using this forecast model.

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.

Vegetation status assessment is crucial for agricultural monitoring and management. Vegetation indices derived from high resolution image time series can be used to derive key phenological parameters for annual crops. In this work, we propose a procedure for the estimation of these parameters and their associated uncertainties. The approach uses Bayesian inference through Markov Chain Monte Carlo in order to obtain the full joint posterior distribution of the phenological parameters given the satellite observations. The proposed algorithm is quantitatively validated on synthetic data. Its use on real data is presented together with an application to real-time within season estimation allowing for phenology forecasting.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

Quantile regression model is widely used in variable relationship research of general size data, due to strong robustness and more comprehensive description of the response variables' characteristics. With the increase of data size and data dimension, there have been some studies on high-dimensional quantile regression under the classical statistical framework, including higher-efficient frequency perspective, which is however at cost of randomness quantification, or lower-efficient Bayesian method based on MCMC sampling. To overcome these problems, we propose the high-dimensional quantile regression with Spike-and-Slab Lasso penalty based on variational Bayesian (VBSSLQR), which can not only improve the computational efficiency but also measure the randomness via variational distributions. The simulation studies and real data analysis illustrate that the proposed VBSSLQR method is superior to or equivalent to other quantile and non-quantile regression methods (including Bayesian and non-Bayesian methods), and its efficiency is higher than any other method.

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.

FY-3E WindRAD (Fengyun-3E Wind Radar) is a dual-frequency rotating fan-beam scatterometer. Its data characteristics, NOC (NWP Ocean Calibration), and wind retrieval performance are investigated in this paper. The radar view geometry diversity varies across the swath, with maximum diversity in the sweet swaths and limited diversity in the outer and nadir swaths. When NOC backscatter calibration coefficients are computed as a function of incidence angle only (NOCint), a smooth correction is found. However, when relative antenna azimuth angle is included (NOCant), it appears that the corrections as a function of relative azimuth angle vary harmonically and substantially for a specific incidence angle. NOCant corrections yield a better fit of the measurements to the GMF (Geophysical Model Function). Hence NOCant is applied for the wind retrieval analysis from Ku-band, C-band. An extra engineering correction of 0.15dB and 0.2dB is applied on Ku-band and C-band backscatter values, respectively, to reduce the wind speed bias without increasing the standard deviation. Overall, NOCant is the best option for both channels. Also, the instrument backscatter data stability over time is good and the retrieved winds can fulfill operational requirements.

Background. COVID-19 efforts were often ineffective in controlling the spread of the pandemic. Identifying ineffective controls during a pandemic is thus vital. Method. Utilizing publicly available data on COVID deaths in the counties of US states, we create an index to capture and interpret ineffectiveness in the efforts to reduce the spread of the pandemic in US counties. This index is based on the Intervened Poisson Distribution (IPD) introduced originally by Shanmugam. Motivation for the research idea occurred while we noticed the data dispersion of the COVID deaths is smaller than the average only in some counties. Under-dispersed data is common in statistical modeling. A novel approach we adapted in this article includes the estimation of an intervention parameter estimated through iterative non-linear optimization. Results. Twenty-five counties in California, Idaho, Minnesota, Mississippi, Montana, Nebraska, North Carolina, North Dakota, Texas, and Utah were found to be ineffective in controlling for fatalities based on the expected probability distribution. A review of the policies enacted in these areas would provide insight into ineffective prevention efforts, and some of these issues are documented in current literature. Conclusion. The IPD index an innovate way to document efficacy of interventions during pandemics. The IPD may identify ineffective efforts prior to statistical models intended to evaluate efficacy of efforts.

Bayesian models allow us to investigate children’s belief revision alongside physiological states like “surprise”. Recent work finds that pupil dilation (or the “pupillary surprise response”) following expectancy-violations may be predictive of belief revision. How can probabilistic models inform interpretations of “surprise”? Shannon Information considers the likelihood of an observed event, given prior beliefs – suggesting stronger surprise occurs following unlikely events. In contrast, Kullback-Leibler divergence considers the “dissimilarity” between prior beliefs and updated beliefs following observations – with greater surprise indicating more change between belief states to accommodate information. To assess these accounts under different learning contexts, we use Bayesian models that compare these computational measures of “surprise” to contexts where children are asked to either predict or to evaluate the same evidence during a water displacement task. We find correlations between the computed Kullback-Leibler divergence and children’s pupillometry responses only when children actively make predictions, and no correlation between Shannon Information and pupillometry. This suggests that when children attend to their beliefs and make predictions, pupillary responses may signal the degree of divergence between a child’s current beliefs and updated, more accommodating beliefs.

In this paper study, we develop a Bayesian Network model selection principle that address addresses the incommensurability of network features obtained from incongruous datasets and overcomes performance irregularities of the Minimum Description Length model selection principle. This is achieved (i) by approaching model evaluation as a classification problem, (ii) by estimating the effect that sampling error has on the satisfiability of conditional independence criterion, as reflected by Mutual Information, and (iii) by utilizing this error estimate to penalize uncertainty in the Minimum Uncertainty (MU) model selection principle. We validate our findings numerically and demonstrate the performance advantages of the MU criterion. Finally, we illustrate the advantages of the new model evaluation framework on a tRNA structural biology example.

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.

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.

Mistrust between scientists and non-scientist stakeholders is a key challenge in fishery management. This problem is exacerbated with the use of complex models to support management: these models suffer from difficulties in communicating their results and a lack of confidence from end users. The European eel is an illustrative example; its complex life cycle raises problems of coordination and discussion among the multiple actors involved in the management of the species. The GEREM model has been proposed as a tool for estimating recruitment, but its complexity, which is essential for addressing the characteristics of the species, makes it difficult to understand and accept by all stakeholders. In the context of the SUDOANG project, we proposed a co-parametrization of this assessment model to tackle this mistrust. Through the use of various questionnaires and appropriate statistical analyses, stakeholders were involved in two important choices for the model (zone definition and prior construction). Regular workshops and presentations were organised to explain the model rationales and to gather feedback and expectations. The results show that stakeholders have very similar perceptions of the potential definitions of sub-areas of recruitment in south-western Europe, and these perceptions are consistent with the underlying environmental conditions. On the other hand, the stakeholders have contrasting opinions about the exploitation rates of fisheries in different river basins, and the use of their knowledge currently has little effect on GEREM estimates. More importantly, the overall approach of this study is thought to have reconstructed the trust and confidence among participants.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We performed Bayesian network meta-analysis (NMA) to suggest frontline treatments for patients with high PD-L1 expression (at least 50%). A total of 5,237 patients from 22 studies were included in this NMA. In terms of progression-free survival, immune checkpoint inhibitors (ICIs) plus bevacizumab plus chemotherapy had the highest surface under the cumulative ranking curve (SUCRA) value (98.1%), followed by ICI plus chemotherapy (82.9%). In terms of overall survival (OS), dual immunotherapy plus chemotherapy had the highest SUCRA value (79.1%), followed by ICI plus bevacizumab plus chemotherapy (73.4%). However, there was no significant difference of survival outcomes among treatment regimens combined with immunotherapy. Moreover, ICI plus chemotherapy failed to reveal a significant OS superiority to ICI monotherapy (hazard ratio = 0.978, 95% credible internal: 0.771-1.259). In conclusion, this NMA indicates that ICI plus chemotherapy with/without bevacizumab might to be the best options in terms of OS for NSCLC with high PD-L1 expression. Considering there was no significant difference of survival outcomes among treatment regimens incorporating immunotherapy and ICI plus chemotherapy failed to show significant survival benefits over ICI monotherapy, however, ICI monotherapy may be reasonable as first-line treatment for advanced NSCLC with high PD-L1 expression and no targetable aberrations.

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.

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.

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.

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.

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.

In this paper, under the symmetric entropy and the scale squared error loss functions, we consider the maximum likelihood (ML) estimation and Bayesian estimation of the Shannon entropy and Rényi entropy of the two-parameter inverse Weibull distribution. In the ML estimation, the dichotomy is used to solve likelihood equation. In addition, the approximation confidence interval is given by the Delta method. Because the form of estimation results is more complex in the Bayesian estimation, the Lindley approximation method is used to achieve the numerical calculation. Finally, Monte Carlo simulations and a real data set are used to illustrate the results derived. By comparing the mean square error between the estimated value and the real value, it can be found that the performance of ML estimation of Shannon entropy is better than that of Bayesian estimation, and there is no significant difference between the performance of ML estimation of Rényi entropy and that of Bayesian estimation.

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.

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.

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.

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.

Background The purpose of this study was to clarify the relationship structure of quality of life (QOL) in survivors of breast cancer, including difficulty in daily life and negative experience as health-related indicators. Methods Participants were survivors of breast cancer for more than 2 years after primary breast cancer surgery and belonged to self-help groups. The assessment used FACT-B (QOL), HADS (anxiety and depression), SOC (sense of coherence), WHODAS 2.0 (difficulties in daily life), and CAOD (negative experiences). Bayesian structural equation modeling (BSEM) was performed to analyze the hypothesized model. If the causal model was significant, multiplication of the path coefficient from anxiety and depression to QOL, and from SOC to anxiety and depression, was considered a direct effect on QOL, and from SOC to difficulty in daily life, from difficulty in daily life to negative experiences, and from negative experiences to anxiety and depression were considered indirect effects on QOL. Results Goodness of fit of the model in the BSEM was satisfactory. The direct effect was 0.274, and indirect effect was 0.163. Conclusions These direct and indirect approaches to difficulties in daily life are very significant because they can provide a basis for development of support strategies to increase QOL.

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).

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.

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.

Exact sequential state estimation of orbiting objects in space can be found by predicting the state using Fokker Planck Kolmogorov Equation (FPKE) and measurement correction using Bayes’ conditional/posterior Probability Density Function (PDF). Posterior PDF can be expressed as Gaussian density expanded in terms of Hermite polynomials named as Gram Charlier Series (GCS). This research is an extension of an earlier work on second order linearized solution to nonlinear Bayesian filtering using Taylor series and third order GCS expansion of posterior PDF. In this new extension, Bayes’ posterior PDF is approximated by a mixture of GCS functions for which the parameters are propagated using linear propagation theory. The update of weights of different components of GCS mixture model uses the FPKE error as feedback to adapt for the amplitude of different GCS components while solving a quadratic programming problem earlier used for Gaussian Mixture Model (GMM) PDF. Proposed filtering method is applied on tracking of space debris. The simulation results for the filter shows performance is moderately better than single GCS filter, Extended Kalman Filter (EKF) and Gaussian Sum Filter (GSF) under space debris’ highly uncertain initial conditions and sparse measurement availability.

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.

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.

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

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.

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.

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.

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.

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.

In actuarial practice, the modeling of total losses tied to a certain policy is a non-trivial task. Traditional parametric models to predict total losses have limitations due to complex distributional features such as extreme skewness, zero inflation, multi-modality, etc., and the lack of explicit solutions for log-normal convolution. In the recent literature, the application of the Dirichlet process mixture for insurance loss has been proposed to eliminate the risk of model misspecification biases; however, the effect of covariates as well as missing covariates in the modeling framework is rarely studied. In this article, we propose novel connections among covariate-dependent Dirichlet process mixture, log-normal convolution, and missing covariate imputation. Assuming an individual loss is log-normally distributed, we develop a log skew-normal Dirichlet process to approximate the log-normal sum. As a generative approach, our framework models the joint of outcome and covariates, which allows to impute missing covariates under the assumption of missingness at random. The performance is assessed by applying our model to several insurance datasets, and the empirical results demonstrate the benefit of our model compared to the existing actuarial models such as the Tweedie-based generalized linear model, generalized additive model, or multivariate adaptive regression spline.

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.

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.

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.

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.

We present a finite element model updating technique for soil-structure system identification of the Millikan Library building using the seismic data recorded during the 2002 Yorba Linda earthquake. A detailed finite element (FE) model of the Millikan Library building is developed in OpenSees and updated using a sequential Bayesian estimation approach for joint parameter and input identification. A two-step system identification approach is devised. First, the fixed-base structural model is updated to estimate the structural model parameters (including effective elastic modulus of structural components, distributed floor mass, and Rayleigh damping parameters) and some uncertain components of the foundation-level motion. Then, the identified structural model is used for soil-structure model updating wherein the Rayleigh damping parameters, the stiffness and viscosity of the soil subsystem (modeled using a substructure approach), and the foundation input motions (FIMs) are estimated. The identified model parameters are compared with state-of-practice recommendations. While a specific application is made for the Millikan Library, the present work offers a framework for integrating large-scale FE models with measurement data for model inversion. By utilizing this framework for different civil structures and earthquake records, key structural model parameters can be estimated from the real-world recorded data, which can subsequently be used for assessing and improving, as necessary, state-of-the-art seismic analysis and structural modeling techniques. This paper presents an effort towards using real-world measurements for large-scale FE model updating in the soil and structure, uniform soil time domain for joint parameter and input estimation, and thus paves the way for future applications in system identification, health monitoring and diagnosis of civil structures.

Previous studies have shown, in the United States (U.S.), that non-White communities are exposed to significantly higher temperatures in urban environments than complementary White populations. Studies highlighting this disparity have usually been cross-sectional and are therefore “snapshots” in time. Using surface urban heat island (SUHI) intensity data, U.S. Census 2020 population counts, and a measure of residential segregation, this study performs a comparative analysis between census tracts identified as prevalent for White, Black, Hispanic and Asian populations and their thermal exposure from 2003 – 2018. The analysis concentrates on the top 200 most populous U.S. cities. SUHI intensity is shown to be increasing through time for all examined tracts. However, the increase is only statistically significant for White and Black prevalent zones. There is a 1.25K to ~2.00K higher degree of thermal exposure on average for non-White relative to White prevalent areas. When examined on an inter-city basis, White and Black prevalent tracts had the largest disparity, as measured by SUHI intensity, in New Orleans, LA, USA by &lt; 6.00K. Hispanic (&gt;7.00K) and Asian (&lt;6.75K) prevalent tracts were greatest in intensity in San Jose, CA, USA. To further explore temporal patterns, two models were developed using a Bayesian hierarchical spatial temporal framework. One models the effect of varying the percentages of each population group relative to SUHI intensity within all examined tracts. Increases in percentages of Black, Hispanic, and Asian populations contributed to statistically significant increases in SUHI intensity. White increases in population percentage lowered SUHI temperature intensity. Throughout all modeled tracts, there is a statistically significant 0.01K per year average increase in SUHI intensity. A second model tests the effect of residential segregation on thermal inequity across all examined cities. Residential segregation, indeed, has a statistically significant positive association with SUHI intensity from this portion of the analysis. Similarly, there is a statistically significant 0.01K increase in average SUHI intensity per year for all cities. Results from this study can be used to guide and prioritize intervention strategies and furthers urgency related to social, climatic, and environmental justice concerns.

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.

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.