Questions tagged [bayesian]

Bayesian (after Thomas Bayes) refers to methods in probability and statistics that involve quantifying uncertainty about parameter or latent variable estimates by incorporating both prior and observed information. Bayesian modeling, inference, optimization, and model comparison techniques are on topic. A programming element is expected; theoretical/methodological questions should go to https://stats.stackexchange.com.

Overview

Bayesian inference is a method of statistical inference which uses Bayes' theorem - named after Thomas Bayes (1702-1761) - to quantify the uncertainty of parameters or latent variables. The statement of Bayes' theorem in Bayesian inference is

enter image description here

Here θ represents the parameters to be inferred and d the data. P(θ|d) is the posterior probability and P(d|θ) is the likelihood function. P(θ) is the prior: a function encoding previous beliefs about θ within a model appropriate for the data. P(d) is a normalization factor.

The formula is used as an updating procedure: as more data become available, the posterior can be updated successively. In the first instance, the prior must be specified by the user. In later updates, the prior is usually chosen to be the posterior from a previous updating procedure.

References

The following threads contain lists of references:

The following journals are dedicated to research in Bayesian statistics:

Tag usage

Questions on tag should be about implementation and programming problems, not about the statistical or theoretical properties of the technique. Consider whether your question might be better suited to Cross Validated, the StackExchange site for statistics, machine learning and data analysis.

1808 questions
0
votes
1 answer

pymc3: Why is everything in terms of log?

I'm trying to wrap my head around pymc3, which seems to be a beautiful package. One thing I don't understand is why all of the probabilities are in log? The description of the Beta Distribution says it is the "Beta log-likelihood". The functions for…
benjaminjsanders
  • 827
  • 8
  • 13
0
votes
0 answers

Problem adjusting a linear model with JAGS

I am trying to adjust a linear model with JAGS but I'm having trouble with the code. I'm writing: library(R2jags) library(BEST) base<-data.table::data.table(read.csv("/Users/franco/Documents/Todo/UNAM/Facultad\ de\ Ciencias/Asignaturas\…
0
votes
1 answer

Multiple questions regarding the KL term in the ELBO loss with TensorFlow Probability

I have been trying to conduct a few experiments using TensorFlow Probability (TFP), and I got a few questions. What is the proper value of the coefficient of the KL loss? In the paper by Blundell (2015), the coefficient is set to 1/M (where M is…
0
votes
0 answers

Adapting Bayesian logistic regression script to my data

I'm looking to run a hierarchical logistic regression in a Bayesian framework, but am having trouble adapting codes for my data. I have the great book "Doing Bayesian Data Analysis", but I'm not sure how to modify the script the author provided…
0
votes
0 answers

How to update different types of distribution estimating same value in R?

we are trying to create a distribution that estimates pathogens presence on vegetables. This was done using different methods, each providing a distribution: - method S (from sludge concentration) is best fitted by weibull(1.55, 8.57) - method SO…
G B
  • 11
  • 1
0
votes
0 answers

what are the parameters of bayes optimization for tuning parameter?

I am using Bayesian optimization to tune the parameters of SVM for regression problem. In the following code, what should be the value of init_grid_dt = initial_grid ? I got the upper and lower bounds of the sigma and C parameters of SVM, but dont…
Neha gupta
  • 43
  • 5
0
votes
1 answer

Missing predictor variable computing Bayesian ANOVA using stan_aov from rstanarm package

I am computing a Bayesian ANOVA to investigate how my heading variable affects FirstSteeringTime. Here is an example of my data set: x <- structure(list(FirstSteeringTime = c(0.433389999999999, 0.449999999999989, 0.383199999999988,…
codegoblin1996
  • 302
  • 2
  • 12
0
votes
1 answer

Converting likelihood equation into R code

I'm working on coding GP model in R and calculating its likelihood function. I'm having problem to convert the equation below into R code: I used the code below but I keep getting error message : Error in (t(zlt - olt)) * (zlt - olt) :…
Reta
  • 363
  • 3
  • 4
  • 15
0
votes
0 answers

Test of Significance for dataset

I'm having trouble in the data analysis for this dataset. So there are clients for pre-enrollment period (-5 to 0) and post-enrollment period (0-5) and I'd like to know the effect of cost between these periods. For example, does the cost increase or…
0
votes
1 answer

Why is my bayesian stats code not working?

This is the problem I'm writing my code off of: You are planning a picnic today, but the morning is cloudy Oh no! 50% of all rainy days start off cloudy! But cloudy mornings are common (about 40% of days start cloudy) And this is usually a dry…
0
votes
1 answer

Comparison of Different Types of Nonlinear Regression Models

Thank you for seeing this post. Various regression models are being applied to the curve estimating (actual measured ventilation rate). Comparison was made using the GLM and GAM models including polynomial regression. I use R. Are there any other…
0
votes
2 answers
0
votes
1 answer

Kalman filter with acceleration. State or Control vector?

I have a basic understanding question in Kalman filter which I haven't found an answer yet. Assume I want to implement a Kalman filter with a constant acceleration dynamic. I can either add the acceleration the state vector and F matrix - Xt =…
Roi Yozevitch
  • 197
  • 3
  • 13
0
votes
0 answers

Coding an integrated Bayesian model with a mix of stochastic and deterministic inputs

The Problem I am having trouble figuring out how to implement a Bayesian Framework for a predictive model that contains many deterministic inputs mixed with a few stochastic inputs. Conceptually the problem seems easy, but from a coding standpoint I…
0
votes
1 answer

Coupled variables in hyperparameter optimization in MATLAB

I would like to find optimal hyperparamters for a specific function, I am using bayesopt routine in MATLAB. I can set the variables to optimize like the following: a = optimizableVariable('a',[0,1],'Type','integer'); But I have coupled variables,…
Marouen
  • 907
  • 3
  • 13
  • 36