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

0
votes
0answers
4 views

Paper about Curse of Dimensionality of Gibbs Sampling

Do you know some references from the literature providing details about the curse of dimensionality of Gibbs sampling? Thanks!
0
votes
0answers
8 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
14 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
22 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
39 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 ...
1
vote
1answer
34 views

Numerical CDF estimation for complicated random variable

Given a combination $U$ of several random variables $X,Y,Z...$ with known distributions, what is an efficient numerical algorithm to estimate PDF or CDF of $U$, if its CDF, PDF, characteristic ...
2
votes
2answers
97 views

Monte Carlo Importance Sampling

I am reading the book on Monte Carlo by Sobol (A Primer for the Monte Carlo Method). In the section on Importance Sampling, he writes: $I = \int_a^b g(x) \: dx$ "to compute this integral, we could ...
1
vote
1answer
36 views

Resolve integral with importance sample Monte Carlo

I'm trying to compute the integral $$\int_{a}^{b}(\sin( 1 + x ) + \cos( 1 + x ))e^{-x}\ dx$$ using importance sample Monte Carlo method. The exercise ask to use Cauchy Distribution to resolve the ...
1
vote
0answers
50 views

What does “sequence is equidistributed in [0, 2]” mean?

I was reading an article in which they are mentioning this sentence: "sequence is equidistributed in [0, 2]" where the sequence in question, is a sequence of real number (the article in question is ...
0
votes
1answer
18 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 ...
2
votes
1answer
37 views

Using Monte Carlo, estimate the value of $$\int_0^1\int_0^1e^{xy} \ dx \ dy.$$

Using Monte Carlo, estimate the value of $$\int_0^1\int_0^1e^{xy} \ dx \ dy.$$ I am using matlab to solve this problem. My code is the following: A = 1; N = 100000; s = 0; for i = 1:N; x = ...
0
votes
1answer
54 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
24 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; ...
1
vote
1answer
18 views

Understanding claim in Newman and Barkema's Monte Carlo book

In Newman and Barkema's Monte Carlo Methods in mathematical physics, on page 23-24, the following claim is made: "Assume we have a function f(x) and the integral $I(x)=\int_0^xf(x')dx'$. Then pick a ...
1
vote
0answers
39 views

Monte Carlo standard deviation of the mean estimate too small.

I'm doing a Monte Carlo calculation and use the standard deviation of the mean $\sigma_M$ as the error. To get an estimate of this from the regular standard deviation I use $$\sigma_M=\dfrac\sigma ...
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
0answers
21 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
2answers
27 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).
1
vote
3answers
48 views

Using loop to approximate pi (Monte Carlo, MATLAB)

I've written the following code, based on a for loop to approximate the number pi using the Monte-Carlo-method for 100, 1000, 10000 and 100000 random points. ...
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
17 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) ...
1
vote
1answer
25 views

Generating two $-1$ correlated Poisson random variables with parameter $5$

Is it possible to generate two random variables $X$ and $Y$ that are both $Poisson(5)$ with $Corr(X,Y)=-1$? Why? I was thinking about generating $3$ independent Poisson random variables $Z_1,Z_2, and ...
0
votes
1answer
64 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 ...
4
votes
1answer
79 views

Buffon's experiment with squares

Say, we'd like to make the Buffon's experiment but with squares instead of needles. Notation: $d$ is the distance between lines $b$ is the square side length $y$ is the distance from the center of ...
2
votes
2answers
87 views

variance of 26 cards chosen from a deck

Suppose I have a well shuffled deck and I am trying to find the variance of 26 cards randomly chosen without replacement from a deck, assuming the values are from 1 to 13 for the cards. Since the mean ...
2
votes
0answers
29 views

Monte Carlo with error on individual samples

I'm performing a Monte Carlo integration where the individual samples have an error, and I'm wondering how to estimate the final error. Some more detail: The integral E I'm after is estimated in the ...
0
votes
1answer
28 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 ...
1
vote
1answer
66 views

Proving that Markov Chain Monte Carlo converges

I am trying to understand how the very basic Markov Chain Monte Carlo approach works: We try to approximately calculate the expected value $E_{\pi(x)}[X]$ by drawing sequential samples from a Markov ...
0
votes
0answers
21 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 ...
3
votes
1answer
347 views

Monte carlo integration in spherical coordinates

I was playing around with writing a code for Montecarlo integration of a function defined in spherical coordinates. As a first simple rapid test I decided to write a test code to obtain the solid ...
3
votes
2answers
90 views

Algorithm to find best in class of groups with weighting?

I have widgets and a single widget will have attributes of: Name Weight (decimal from 0-1) Group (letter A-F) Price (an integer from 1 - 100) I must pick one ...
0
votes
0answers
9 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
72 views

Sample uniform direction within cone

My question is pretty much the same as this question below, however I came up with a potential solution to this problem that I didn't see an answer to in the other question and I was wondering if it ...
0
votes
1answer
31 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
1answer
32 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 ...
4
votes
2answers
204 views

Monte Carlo estimator of the number of 1's in a very long binary sequence

Preface The question below is related to a problem I am working on, which requires counting the number of times a logic-valued function evaluates "TRUE" given an input value. The size of my input set ...
0
votes
0answers
22 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
1answer
97 views

Obtaining useful information from graph obtained via Monte-Carlo Simulations

I've been running Monte Carlo Simulations on some Matlab code and then plot the graph shown below. I was just wondering what useful information I could collect from this graph? Edit: fit ...
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
21 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 ...
1
vote
1answer
31 views

How do we calculate when to shout for the optimal payoff?

For example: A non-dividend paying stock is currently priced at $20, and you hold a put that allows early exercise in 2 months and in 4 months. The option expires in 6 months. Volatility is 30%, and r ...
3
votes
3answers
188 views

Markov chain stationary probability simulation

Having a defined markov chain with a known transition matrix, rather than to calculate the steady state probabilities, I would like to simulate and estimate them. Firstly, from my understanding there ...
2
votes
2answers
52 views

How many simulations for a game?

Lately I was interested by Monte-Carlo simulations. I found many papers about this approach in the Internet but for now they are too hard for me. I just want to start understanding this method with ...
1
vote
1answer
94 views

Monte carlo estimation of maximum likelihood estimators

I'm interested in numerically finding the maximum likelihood estimator of a parameter $\theta$, as well as the confidence interval of this estimator. First I'll describe the method I've been trying, ...
0
votes
0answers
30 views

Average over all positive functions on the unit interval whose Lebesgue integral is one

I want to average over all positive functions on the unit interval whose Lebesgue integral is one. Formally, I want to compute the mean of the following probability distribution defined over function ...
1
vote
1answer
26 views

Bias Method in monte carlo integration

This is from a proof in my monte carlo course. let $h$ be a smooth function, $T_n = h(\bar{X})$ $\mu = E(X)$ then by taylor expansion $E(T_n -\tau) = E[h(\bar{X} -h(\mu)] = E[\bar{X} - \mu]h'(\mu)+ ...
7
votes
3answers
531 views

integration method

I want to calculate an integral by using the hit and miss method. I can not understand how this method works. I would be grateful if someone could explain me and help me to calculate the value, with a ...
4
votes
0answers
48 views

$\pi$ Monte-Carlo - Probability that O-Lock hit a Spoke?

(Edit: can someone please help me migrate this to physics stack? I think they would be more interested in helping me out with this problem. Thanks.) I have a bicycle with one of those O-locks on it ...
2
votes
1answer
53 views

Numerical integration of innocent-looking singular integrand

Consider the rather innocent integral: $$I=\int_{0}^{1}a x^{a-1}dx=1,\quad 0<a<1$$ Numerically, this integral converges awfully slowly, and one must use a recursive method to get anywhere near ...
1
vote
1answer
67 views

Monte Carlo estimations of e

I need to estimate $e$ with a monte carlo method. We only learned the crude monte carlo integration, so I can't use any robust monte carlo simulations. I know that $\displaystyle \int\limits_1^x ...