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
139 views

Volume with MonteCarlo

I'm asked to find the volume of a given bounded solid in $\mathbb{R}^3$ by MonteCarlo means: since it is contained in a prism, I generate random points in the prism and see what proportion of them lie ...
1
vote
0answers
104 views

Approximation of SDE

I have been struggling with the following problem: If you want to find a numerical result by simulating the paths of a stochastic differential equation, in particular a geometric brownian motion I ...
1
vote
0answers
217 views

Computing Margin-of-Error using Monte Carlo simulations

I am interested in computing the margin-of-error for a metric computed on a random sample. The underlying distribution (finite) from which the random sampling is done is not normal (its extremely ...
1
vote
0answers
107 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
46 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
58 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
40 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
197 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
320 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
1answer
15 views

Random sampling using Metropolis vs. Accept-Reject

I'm working on a project comparing the accept-reject (von Neumann) and Metropolis methods of sampling. I'm generating a sample of size $N$ from a normal distribution $N(1,(\frac{1}{2})^2)$. What I ...
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
2answers
46 views

How to use Monte Carlo method to find volume?

So the question asks to use the Monte Carlo method to find the volume of an irregular figure defined as: ...
0
votes
1answer
36 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
441 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
97 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
426 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
115 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
366 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
80 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
89 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
62 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
143 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
46 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
145 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
120 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
36 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
432 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
16 views

How to rewrite function for squared uniform distribution

The question is as follows: I am evaluating the following integral: $$\int_o^1\frac{\exp(\sqrt{1-x^2})}{\sqrt{x}}dx$$ by assuming it equals $E[f(U)]$ for a uniform distribution. I worked it out via ...
0
votes
1answer
11 views

Logic with increasing monte carlo possible output

I am working on Monte Carlo algorithm : Given that you have an experiment MC which has a p-correct of 75%, which means it gives you the right answer 75% of the time. You run MC3 which repeats MC 3 ...
0
votes
0answers
12 views

SEnsitivity Indices are non zero

I am trying to compute the sensitivity indices (SI) of a function using Monte Carlo simulation. I had written a matlab code that perform the computation directly and just return the final answer of my ...
0
votes
0answers
18 views

Is there a such thing as a quasi-random shuffle?

I've recently experimented with Quasi-random numbers in monte-carlo applications. Is there a way to construct a quasi-random shuffle? By that I mean can I take a sequence $Q$ and shuffle it to produce ...
0
votes
0answers
28 views

Monte Carlo Markov Chain simulation

I am going to post the python code logic we used however I want someone to look at the number that are printing out. The Markov chain is uniformly distributed across all $50x50$ matrices with entries ...
0
votes
0answers
8 views

Using the rejection method to generate values

I'm trying to solve the following problem on rejection sampling: I think I have a good idea about what rectangle I should be using. In my mind, it would be a rectangle just large enough to encompass ...
0
votes
0answers
7 views

Sensitivity analysis of paramaters and input variables

I am trying to perform a sensitivity analysis of an optimization problem $f(x,\alpha)= \min_{ Q} {g(x,\alpha , Q)}$ where $x$ is an input variable for our function, and $\alpha $ is a parameter. ...
0
votes
1answer
23 views

Compute the Std. Deviation of Multiple Monte Carlo Estimation of $\pi$

For a school programming assignment, I am trying to compute the value of $\pi$ via the classic Monte Carlo estimation of $\pi$. In the experiment, we throw a variable number of darts at a circle that ...
0
votes
0answers
15 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
26 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
52 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
32 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
39 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
33 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
44 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
11 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
27 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
74 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 ...