Skip to main content

All Questions

Tagged with
1 vote
0 answers
12 views

Recycling MCMC samples for another data set from the same distribution

Suppose I'm given $\theta_0$ and I want to sample data from a density $f(Y|\theta_0)$ and then sample from the posterior of $\theta|Y$ (given, obviously, some prior). I want to do this lots of times, ...
Thomas Lumley's user avatar
5 votes
1 answer
42 views

ABC (Approximate Bayesian Computation) Sampling, Simulating data from Complex models

In ABC sampling methods, Rejection, MCMC and SMC, when we sample potential parameter values from the prior/proposal, we then use those parameters on our model and simulate data values. This can be ...
AlexS123's user avatar
0 votes
0 answers
119 views

How to compute Expected Squared Jump Distance (ESJD) of a Metropolis-Hastings algorithm

The Expected Squared Jump Distance (ESJD) seems to be defined slightly differently in various papers, which makes this very confusing. For instance, Definition 2.2 of Optimal Scaling of Random-Walk ...
Euler_Salter's user avatar
  • 2,236
1 vote
0 answers
18 views

Entries of the matrix variate generalized inverse Gaussian

I'm currently dealing with a Gibbs sampler of the matrix variate generalized inverse Gaussian distribution. In order to check the correctness of the Gibbs sampler, I'd like to know what are the ...
Stéphane Laurent's user avatar
1 vote
1 answer
37 views

In what sense do Bayesian models simulate observation?

Could you verify that i understand the author correctly here? New to Bayes. [...] These assumptions [the prior and the likelihood] also allow us to simulate the observations that the model implies. ...
rtviii's user avatar
  • 111
0 votes
1 answer
89 views

Acceptance-Reject to generate a distribution proportionate to Inverse Gamma and truncate Cauchy distribution

Assume $Y_{ij} \sim N(\mu_i,\sigma^2)$, $\mu_i \sim N(\eta,\tau^2)$ for $i=1,2$ $j=1,\cdots,n_i$ and prior $\pi(\eta,\tau^2,\sigma^2) \propto Ca^+(\tau^2,0,b_{\tau}) \times Ca^+(\sigma^2,0,b_{\sigma})$...
Justin 's user avatar
5 votes
1 answer
201 views

Sampling from the posterior with a constraint on the posterior mean

Background Under certain assumptions we know that being given the posterior mean and a family of conditional distributions, we can uniquely determine the joint distribution. I quote one of the ...
treskov's user avatar
  • 540
2 votes
0 answers
73 views

Bayesian calibration of computer simulations - Likelihood function calculation

I am starting to study Bayesian calibrations of computer models. I am not a statistician and just starting to learn so bear with me if I do not use the correct terminology. The general approach is ...
Rojj's user avatar
  • 121
2 votes
0 answers
217 views

Why is it easy for the Gibbs sampler to take long time to converge to target distribution?

This is related to Gelman's Bayesian Data Analysis 3rd Edition pg 300 first paragraph of Section 12.4. The book says the following. "An inherent inefficiency in the Gibbs sampler and Metropolis ...
user45765's user avatar
  • 1,445
1 vote
0 answers
146 views

R statistics: compare bayesian bootstrap to frequentist bootstrap for statistics: univariate odds ratio for small sample [closed]

Greetings to the community, I am seeking assistance in finding a solution to the challenges I am facing. OBJECTIVES: I aim to estimate the univariate odds ratio for a binary exposure in a population. ...
FEVE's user avatar
  • 11
0 votes
0 answers
75 views

How to use Gibbs sampler to simulate normal-normal hierarchical models?

This is related to Gelman's BDA 3rd Edition Chapter 11, Sec 3. The book says the following. "The Gibbs sampler is the simplest of the Markov chain simulation algorithms, and it is our first ...
user45765's user avatar
  • 1,445
0 votes
0 answers
26 views

How to simulate multivariate posterior distribution with a flat prior in general?

If I know that the posterior $p(\theta_1,\dots,\theta_m|y)$ can be written $p(\theta_1|\theta_2,\dots,\theta_m,y)p(\theta_2|\theta_3,\dots,\theta_m|y)\dots p(\theta_m|y)$ where $p(\cdot|y)$ in each ...
user45765's user avatar
  • 1,445
1 vote
0 answers
57 views

Why is my 95% CI from a simulation being either 100% or 0%?

So what I intended doing with the code is to simulate dataset with the outcome being number of parasites that died (yi) out of every n=50 parasites on the host. There are 2 treatment groups (control ...
soothe's user avatar
  • 11
0 votes
0 answers
22 views

Identify name of sampling / simulation strategy

From what I've learned, it seems like in order to simulate draws from a distribution $X \sim p(X)$ one can take advantage of the fact that $p(X) = \int p(X, z)p(z)dz$ and use the following strategy: ...
DancingIceCream's user avatar
1 vote
0 answers
11 views

circularness of imputing draws via regression in MCAR

Suppose we have covariates $X$ and outcome $Y$. Part of $X$ is missing under MCAR. I am reading Rubin's Missing Data book. I do not see why the following statement making sense intuitively. Let $x,y$ ...
user45765's user avatar
  • 1,445

15 30 50 per page
1
2 3 4 5
8