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.

Background: Aedes albopictus, the Asian Tiger Mosquito, has become a prevalent pest in Italy, causing severe nuisance and posing a threat of transmission of arboviruses introduced by infected travellers. In this study, we investigated the influence of weather parameters on the seasonal population density of Aedes albopictus. Methods: A Bayesian approach was employed to identify the best meteorological predictors of species trend, using the eggs collected monthly from 2010 to 2022 by the Emilia-Romagna regional monitoring network. Results: The findings show that the winter-spring period (January to May) plays a crucial role in the size of the first generation and seasonal development of the species. Conclusions: A temperate winter and a dry and cold March, followed by a rainy and hot spring and a rainy July, seem to favour the seasonal development of Ae. albopictus.

This paper investigates the role of English language proficiency (ELP) as a measure of cultural integration on U.S. immigrant families’ economic integration, as measured by annual earnings, and health care market integration, as captured by annual medical care spending. Bayesian Markov Chain Monte Carlo (MCMC) estimation of a mixed bivariate ordered probit (MBOP) model of annual earnings and medical care spending is conducted using a nationally representative pooled cross-sectional sample data from the U.S. National Health Interview Surveys (NHIS). The results reveal that limited English proficiency is associated with lower marginal earnings propensities (-$44.83, -$50.66, and -$56.97), and higher marginal medical care spending propensities ($8.09, $6.09, and $4.32) for each of the three decreasing levels of ELP. Furthermore, we find a small yet statistically significant positive 7.2% correlation between immigrant families’ annual earnings and medical care spending propensities. Therefore investing to raise immigrants’ families ELP in the U.S. would not only contribute to more effective socio-economic integration, but also strengthen the U.S. workforce and economy for decades to come.

This paper proposes the use of Bayesian inference techniques to search for and obtain valid instruments in dynamic panel data models where endogenous variables may exist. The use of Principal Component Analysis (PCA) allows for obtaining a reduced number of instruments in comparison to the high number of instruments commonly used in the literature, and Monte Carlo Markov Chain (MCMC) methods enable efficient exploration of the instrument space, deriving accurate point estimates of the elements of interest. The proposed methodology is illustrated in a simulated case and in an empirical application, where the partial effect of a series of determinants on the attraction of international bank flows is quantified. The results highlight the importance of promoting and developing the private sector in these economies, as well as the importance of maintaining good levels of credit worthiness.

This research paper examines a cosmological model in flat space-time via $f(R,G)$ gravity where $R$ and $G$ are respectively the Ricci scalar and Gauss-Bonnet invariant. Our model assumes that $f(R,G)$ is an exponential function of $G$ combined with a linear combination of $R$. We scrutinize the observational limitations under a power law cosmology that relies on two parameters - $H_0$, the Hubble constant, and $q$, the deceleration parameter, utilizing the 57-point $H(z)$ data, 8-point BAO data, 1048-point Pantheon data, joint data of $H(z)$ + Pantheon, and joint data of $H(z)$ + BAO + Pantheon. The outcomes for $H_0$ and $q$ are realistic within observational ranges. We have also addressed Energy Conditions, $Om(z)$ analysis and cosmographical parameters like Jerk, Lerk and Snap. Our estimate of $H_0$ is remarkably consistent with various recent Planck Collaboration studies that utilize the $\Lambda$CDM model. According to our study, the power law cosmology within the context of $f(R,G)$ gravity provides the most comprehensive explanation for the important aspects of cosmic evolution.

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.

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.

Semi-continuous data are very common in the social science and economics. In this paper, a Bayesian variable selection procedure is developed to assess the influence of exogenous factors including observed and unobserved on the semi-continuous data. Our formulation is based on the two-part latent variable model with polytomous response. We consider two schemes for the penalties of regression coefficients and factor loadings: the Bayesian spike and slab bimodal prior and the Bayesian lasso prior. Within the Bayesian framework, we implement Markov Chains Monte Carlo sampling method to conduct posterior inference. To facilitate posterior sampling, we recast the logistic model in part one as the norm-like mixture model. Gibbs sampler is designed to draw observations from the posterior. Our empirical results show that with suitable hyper-parameters, the spike and slab bimodal method slightly outperforms the Bayesian lasso in the current analysis. Finally, a real example related to the China household finance survey is analyzed to illustrate the application of the methodology.. .

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.