Questions on Monte Carlo methods, methods that require the repeated generation of (pseudo-, quasi-)random numbers for computing their results.

learn more… | top users | synonyms

1
vote
0answers
106 views

Three ideas of perfect sampling

From David J.C. MacKay's Information Theory, Inference, and Learning Algorithms 32.2 Exact sampling concepts Propp and Wilson's exact sampling method (also known as "perfect simulation" or ...
1
vote
0answers
45 views

How to sample the walk which visits each vertex of a graph specific number of times?

Is there any MCMC mathod that allow me to uniformly sample from all feasible walks where the following restrictions apply: ...
1
vote
0answers
57 views

Gibbs Sampling versus General Cases

Suppose we are given a prior distribution about an unknown parameter $\pi(\theta)$. Also we are given $f(x_{1}, \dots, x_n|\theta)$. We want to find $\pi(\theta|x_1, \dots, x_n)$. Now ...
1
vote
0answers
162 views

Markov Chains - Using Gibbs & Metropolis algorithm.

Suppose $f_x,_y$ is bivariate normal distribution. I was given the parameters $(μ_1, μ_2, σ_1^2, σ_2^2)$ and $ρ=0.95$ the correlation coefficient. I want to generate $(x_1,y_1), ...
1
vote
0answers
39 views

Integrating with respect to an MCMC Kernel

I'm currently trying to evaluate an integral with respect to an MCMC kernel, and I wanted someone else to read what I `proved' to see if it's really correct. It just doesn't feel right. Suppose I ...
1
vote
1answer
193 views

Will I have learned the prerequisites for self learning stochastic calculus and monte carlo method?

I'm an undergraduate econ major, and my main focus is in actuarial sciences, which as you may or may not know it's pretty mathematical. Some of the topics I will have to learn at some point on my own ...
0
votes
1answer
311 views

Why are random numbers necessary for a Monte Carlo simulation?

This may be somewhat of a question with an obvious answer, but I can not seem to understand the necessity of "truly" random numbers to make a Monte Carlo simulation a good one. I understand that not ...
0
votes
2answers
42 views

Numerical estimation of simple integral

Considering the problem of numerical evaluation of the integral of a 'good' function $f(x)$ over a unit interval $I = \int_0^1f(x)dx$ Why can we say $I = E[f(U)]$, where $U\sim Uniformly[0, 1]$?
0
votes
1answer
34 views

Determining Errors in Monte Carlo Simulation

I was wondering if anyone could throw light on possible errors associated with Monte Carlo sampling. I seem to be getting values that are slightly different each time despite running my model for ...
0
votes
2answers
411 views

How to recognize ellipse/ellipsoid from random points? UN-weighted average?

Suppose we are getting random points in 2D (or 3D) which tend to be on ellipse (or ellipsoid). We can't guarantee points are uniformly distributed over ellipse (ellipsoid surface). The task is to ...
0
votes
1answer
96 views

Logic Question Regarding Sample Number and Time Increase

So this question is one of those "I'd rather ask and look stupid now than never know" types of questions. It goes as follows: The error in a Monte Carlo estimate is dominated by ...
0
votes
1answer
403 views

Combining two triangular distributions to yield one distribution

I am interested in using some Monte Carlo methods to help with an estimation problem I have. I need to allow multiple estimators to estimate line-items giving a best, average and worst case estimate ...
0
votes
1answer
89 views

Double Integrals & Expected Value Monte Carlo Method

Tell me if I'm wrong Let $\Omega = [a,b]\times[c,d]\subseteq\mathbb{R}^2$, then $$ \iint_\Omega ...
0
votes
1answer
151 views

Evaluating double integrals using monte carlo methods in matlab.

I used the monte carlo method to integrate $\int_{0}^{1}x^2dx$ in matlab. My matlab code was simply the following: A=1; N=10000; s=0; for i=1:N x=rand; y=rand; if y<= x^2; ...
0
votes
2answers
28 views

Let Y be a random variable with $0\le Y\le 1.$ [duplicate]

Let Y be a random variable with $$0\le Y\le 1.$$Show that $$var(Y)\le 1/4 $$ and that $$var(Y)= 1/4 $$ if and only if P(0)=1/2=P(1).
0
votes
1answer
69 views

How to mathematically prove that we are sampling from same distributions?

The content of this question is about rigorously proving something which is otherwise considered easily correct intuitively. Let's assume we have a multivariate distribution $g(x_1,x_2,...,x_n)$ over ...
0
votes
1answer
78 views

Monte Carlo method error in Bernoulli random variables

Assume I am flipping an unfair coin. Flipping the coin will be heads with probability $p$ and tails with $1-p$. I have no idea what $p$ is (it could even be $.5$!) Let's say I decide to use the Monte ...
0
votes
1answer
58 views

Monte Carlo estimator

I have hopefully a short/simple question regarding monte carlo estimators. The expected value of a function of a random variable can be defined as: $$E[f(x)] = \int_{-\infty}^{\infty} f(x) p(x) dx$$ ...
0
votes
1answer
139 views

Alternatives to Monte-Carlo simulation

Imagine I have a model of economy of a region, which consists of several companies, importers and population. Let's assume that all local companies in question produce food and agricultural ...
0
votes
1answer
41 views

Confusion about Monte Carlo integration

I find I can not really understand the Monte Carlo integration, even I use it for many applications, like stochastic ray tracing. Let us take circle-area-calculation for an example, First, we think ...
0
votes
1answer
128 views

Monte Carlo integration, expected value of the sample mean and expected value of f(x)

I am still progressing in my learning of probability and monte carlo method. I understand a basic MC estimator can be written as: $$\bar x = { 1 \over N } \sum_{i=1}^N f(x_i) \approx E[f(x)]$$ I ...
0
votes
1answer
114 views

Monte Carlo sampling a binomial expansion

I want to figure out the following question $$ 1 = (10 - 9)^{100} = 10^{100}-100 \times 10^{99} \ 9 + \frac{100 \times 99}{2} 10^{98} \ 9^{2} - \frac{100 \times 99 \times 98 }{3}10^{97} 9^{3} ...
0
votes
1answer
35 views

Why cannot the Markov Chains used in MCMC simulations be null recurrent?

I am aware this question borderlines retardedness, but I am seeking an accurate explanation. I understand in null-recurrent cases, the expected amount of time to explore states can be infinite. Is ...
0
votes
2answers
412 views

How to decide what is the probability distribution in a Monte-Carlo simulation?

For a Monte-Carlo integration of $$\int_\Omega P(x)f(x)\ \text d x,$$ there seems to be no apriori distinction if $f$ or $P$ is the probability function. So does it matter if I consider $$P, f, P ...
0
votes
2answers
104 views

Finding a distribution of a random variable generated using a Monte Carlo method

I would greatly appreciate if somebody could confirm or negate my result to the following problem. I am especially not sure about "putting it all together" step. Generate $U_1,U_2,U_3 \sim ...
0
votes
0answers
6 views

Processing Issues [on hold]

I have to run heaps of models using Monte Carlo simulation so in essence it takes forever. Are there any websites that I can rent processing time so that I can get it done really quickly? Cheers for ...
0
votes
0answers
14 views

How to sample points from a bounded polytope?

I have a bounded polytope $C \subset \mathbb{R}^n$ characterized by the following restraints: $$ x \in C \Leftrightarrow \sum_{i=1}^n x_i = 1 \text{ and } Ax \leq b$$ for some matrix $A \in ...
0
votes
1answer
16 views

Approximation and Monte Carlo simulation.

I am a bit up over my head here, I will present an argument and then I hope you guys will say if my reasoning is correct or what should be changed, ultimately I am hoping to say something qualified ...
0
votes
0answers
23 views

Correlated variables from Latin Hypercube

Say I have a vector $\mathbf{Y}$ of $n$ normally distributed random variables. I have its mean vector $\mu$ and covariance matrix $\Sigma$. Normally if I were to generate a sample, I would decompose ...
0
votes
0answers
6 views

cholesky decomposition and generate random vector

Hi someone please explain me this code : it is about running stress scenario simulation with two hypothetical shocks on macroecomics variables (gdpgap- exchange rate). The first hypothetical ...
0
votes
0answers
21 views

Optimizing Buffon's needle to minimize the variance in $\pi$

This is a homework problem so I don't expect a full answer, I'm looking for some pointers on where to start. Problem text: Find L (stick length) and D (separation between lines) that minimizes the ...
0
votes
0answers
26 views

Calculating error bounds for linear regression fit using Monte Carlo methods

so my question is pretty simple. I have some data that has a known error in the y coordinate and I'm fitting it to a linear model using least-squares. Now normally I know we neglect the error and ...
0
votes
0answers
28 views

Multi-armed bandit optimization

I've got a variation of multi-armed bandit problem, where I need not to minimize the regret, but to find a bandit with a maximum reward. Could you please tell how this problem is called or suggest ...
0
votes
0answers
42 views

Sampling and averaging in Monte Carlo Simulation

(First of all, I apologize for the vague title. Couldn't think of rather proper one.) Let's say that we have 10 items where each item has probability distribution of one's own, say Lognormal ...
0
votes
0answers
9 views

Metropolis Monte Carlo with modified acceptance

What happens, if I change the acceptance criterion in a Metropolis Monte Carlo algorithm? I do know the classic proof of detailed balance, which for symmetric transition matrices gives a set of ...
0
votes
0answers
19 views

Simultaneous multiple perturbations in Markov chain Monte Carlo

I'm coding a McMC algorithm for geophysical applications. Using the Metropolis-Hastings scheme to accept/reject the proposed models is smth that i thought i completely understood, but i don't. To be ...
0
votes
0answers
25 views

Correlating random numbers seems to skew the data

I am trying to generate a series of correlated random numbers that represent currency exchange rates for a Monte-Carlo simulation. I am attempting to do this via a Cholesky decomposition of the ...
0
votes
1answer
63 views

Forecasting future revenue and expensces

I am trying to forecast future revenue and expenses in a company. In the past I used moving average method but later I am more inclined to try to do that by using monte carlo simulation. I am ...
0
votes
1answer
19 views

Monte Carlo Methods for non-orthogonal functions

I'm trying to approximate a function using a set of piecewise polynomials. For example, perhaps I'd like to uniformly split the domain [-1,+1] 20 times and place Wendland RBF, Gaussian, or maybe a ...
0
votes
0answers
17 views

Showing that the variance increases with the dimension of the random vector

This is actually related to a more complex question; but I want to re-ask it by trying to simplifying it as possible: 1- We have $n$ dimensional functions of the form $f_n:\mathbb{R}^{n} \mapsto ...
0
votes
1answer
31 views

Uniform convergence of Monte Carlo approximation

Usually Monte Carlo method is used to compute integration. For example, let $g(x,\theta)$ be a continuous function about $x$ and $\theta$, $f(x \mid \theta)$ is a continuous pdf with parameter ...
0
votes
0answers
19 views

How to show that the variance increases with the dimension $n$?

This can be seen as a statistics related question, but it is actually a more general mathematics related one. I am trying to understand the Particle Filter and the motivation to use it over the ...
0
votes
0answers
20 views

Monte carlo error: Combining “experimental” and statistical errors

I'm doing a slightly involved Monte Carlo approximation of a quantity $E$ where I end up with the following formula: $E=\frac{\sum_{i=1}^np_ie_iG_i}{\sum_{i=1}^np_iG_i}\ .\ \ \ \ \ \ \ \ \ $ (1) ...
0
votes
1answer
34 views

Calculate expectation under risk neutral measure: $\mathbb{E_Q}(\max(S-1,0))$

I am busy with a numerical simulation and I want the calculate the following expectation under the risk neutral measure: $\mathbb{E_Q}(\max(S-1,0))$. $S$ is some variable that I calculated using ...
0
votes
0answers
30 views

Monte carlo formula to compute the approximation of variance of MLE

In the book of "Monte Carlo Statistical Methods", the book gives an approximation formula for the variance of MLE, Later on, the book mentions that this approximation formula can be written as ...
0
votes
0answers
10 views

control variates - estimating definite integrals

What are the good techniques to find $g(x)$ so that $f (x) - g(x)$ is minimal, in order to evaluate $\int_a^b [(f(x) - g(x)) + g(x)]dx$?
0
votes
1answer
32 views

Analytic approach to find probability and total value of a set of independent events

I have a forecasting worksheet which describes a set (worksheet) of independent events, all of which have a likelihood of happening given as a probability (e.g. 0.7). Every event also has a yield ...
0
votes
0answers
23 views

Antithetic pair of non-independent normal random variables

Suppose that I have two non-independent normal random variables, X and Y such that $(X,Y)$ has mean 0 and the following variance covariance matrix: \begin{bmatrix} 1 & \rho ...
0
votes
0answers
26 views

Understanding graph obtained from Monte-Carlo simulations

I am running a Monte Carlo Simulation where I sample from about 65 Normal Distributions. I also keep track of the probability associated with each sample by approximating a thin area in the Normal ...
0
votes
0answers
26 views

low discrepancy of halton sequences

I want to proof, that the Halton sequence is low discrepancy. I have to show that $$D_N^*(\mathcal{S})\le C\frac{\ln(N)^s}{N}$$ where $D_N^*$ ist the star discrepancy, $\mathcal{S}$ is the Halton ...