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

3
votes
0answers
31 views

Estimating/approximating a very high dimensional unbounded poisson's equation

Consider the poisson equation on an unbounded domain. Suppose that the solution is known to exist. $$ \Delta u=f $$ I would like to estimate the solution of the this equation at a given point $x_0$. ...
0
votes
0answers
30 views

MCMC algorithm for solving MAXSAT [closed]

Design a MCMC algorithm for solving MAXSAT Has anyone tried it if so can you elaborate on your strategy ???
0
votes
0answers
9 views

Volume of a region given by a CSP

I have a Linear Constraint Satisfaction Problem i.e. I have variables $ x_1, x_2,...,x_m$, with corresponding domains $D_1,D_2,...,D_m $ satisfying linear constraints $C_1, C_2,...,C_n$ with $n ...
2
votes
2answers
28 views

What's the average length of a random line segment in a $1 \times 1$ field?

What is the average length of a line segment in a $1 \times 1$ field? Given $$x_1, y_1, x_2, y_2 \in [0,1]$$ $$S = (x_1,y_1,x_2,y_2)$$ $$dist(x_1,y_1,x_2,y_2) = \sqrt{(x_1-x_2)^2+(y_1-y_2)^2}$$ ...
0
votes
0answers
12 views

Forecasting disputed transaction frequencies

Problem I would like to forecast credit card chargeback/dispute frequencies using historical dispute data I have recorded over time. The data I currently store includes: Disputed transaction date ...
-1
votes
0answers
18 views

two different Monte Carlo approaches

Assume that the function $f$ is integrable and maps $[0, 1]$ into $[0, 1]$. Consider estimating $\int_0^1 f(x)\,dx$ using two different Monte Carlo approaches. The standard approximation is applied in ...
0
votes
0answers
14 views

Estimating the Square of a Mean

Suppose I want to estimate $\theta = (\mathbb{E}[f(X)])^2$, where $f: \mathbb{R} \to \mathbb{R}$ and is Borel-measureable, and $X$ is a random variable. I'll use Monte Carlo, for which one ...
2
votes
1answer
30 views

Metropolis Hastings

So I have seen the Metropolis Hastings algorithm written 2 ways, and I don't quite understand how they can be equivalent: The first way is by defining the 'acceptance probability' as: ...
2
votes
0answers
72 views

Question about Random Walks and An $O^*(n^5)$ Volume Algorithm for Convex Bodies - Kannan Lovasz Simonovits 97

I've been trying to understand this paper: "Random Walks and An $O^*(n^5)$ Volume Algorithm for Convex Bodies", Ravi Kannan, Laszlo Lovasz, Miklos Simonovits. Motivation: The paper is about ...
1
vote
0answers
24 views

MCMC/E-M limitations?MCMC over E-M?

I am currently learning hierarchical bayesian models using JAGS from R, and also pymc using ...
0
votes
2answers
33 views

Approximating an integral using Monte Carlo Method

I wrote a solution for Calculate the value of the integral I = $\int_0^\pi sin^2(x)dx$ using the Monte Carlo Method (by generating $ 10^4 $ uniform random numbers within domain [0, π] × [0, ...
1
vote
1answer
37 views

Approximate an integral using Monte Carlo method

I have a question on an assignment Calculate the value of the integral I = $\int_0^\pi sin^2(x)dx$ using the Monte Carlo Method (by generating $ 10^4 $ uniform random numbers within domain [0, ...
1
vote
3answers
68 views

Evaluating Difficult Monte Carlo Integration in R

I recently posted a simple version here (Simple Monte Carlo Integration). I was able to verify that the answer was indeed close to 1/3 when I wrote the following R code, and got a mean of X of ~1/3: ...
0
votes
1answer
53 views

Monte Carlo gamma function

This question was asked before but I'd like to ask something more precise given the answer that was given. [ Estimate gamma function using monte carlo ] What is the criterion for a random point to ...
2
votes
1answer
61 views

Simple Monte Carlo Integration

I am trying to use Monte Carlo Integration, which is nicely described in the answer here (Confusion about Monte Carlo integration). I am using Monte Carlo Integration to evaluate $\int_0^1x^2\,dx$. ...
0
votes
2answers
67 views

Estimate gamma function using monte carlo

Let $\Gamma(\beta) = \int_0^\infty x^{\beta - 1} e^{-x} dx$ how to estimate the above gamma function using monte carlo? Any idea?
0
votes
0answers
19 views

Sequential Importance Sampling

What are the pros and cons of the the basic SIS-algorithm? I know there is some drawback considering the weight degeneration, but not so much about the pros. Also, is there a proof that the extra ...
0
votes
0answers
12 views

Importance Sampling Distribution

I have an infinite set of events and these event are either ture or false. I perform a monte carlo simulations to find the probability of an event being true. Now I have the knowledge that $A$ % of ...
0
votes
0answers
5 views

Why is importance sampling always framed as calculating E(h(x))?

In all of the tutorials I've seen, importance sampling is always framed as a way to calculate: $$E(h(x)) = \int h(x) f(x) dx \qquad x \sim f$$ I don't understand why it is not framed as a more ...
1
vote
0answers
46 views

Random numbers generator

If I know how to generate random numbers from Gaussian distribution (using Box-Muller method), how can I generate random numbers from distribution with pdf ...
0
votes
1answer
16 views

What is the average minimum distance between two Sobol points?

Having the first n points of a d-dimensional Sobol sequence, what is the average Euclidean distance from one arbitrarily point to its nearest neighbour?
1
vote
2answers
123 views

Random directions on hemisphere oriented by an arbitrary vector

Hy, i'm writing a raytracer, and for that I need to generate n random vectors that are inside an hemisphere oriented by the surface normal. Ideally, I would also like being able to restrict the rays ...
1
vote
2answers
59 views

Probabilities seem to be growing exponentially

We have instituted random drug testing at our company. I was charged with writing the code to generate the weekly random list of employees. We've gotten some complaints because of people getting ...
0
votes
1answer
56 views

Histogram normalization

I need to generate random numbers from Gaussian distribution and to draw an equalized histogram. I've generated them in Matlab using Box-Muller transformation. Since I wasn't sure how to equalize the ...
1
vote
0answers
28 views

Standard deviation of the mean through bootstrap resampling of dependent samples

I'm trying to do a Monte Carlo approximation of an integral where the samples are not independent (how much so can be tuned by a parameter giving how often I sample). Therefore the regular expression ...
1
vote
0answers
56 views

Monte Carlo sampling and Kolmogorov–Smirnov test in practice

I have two deterministic algorithms, Algorithm 1 and Algorithm 2. The first has $m$ inputs and one output, and the second has $n$ inputs and one output. The distributions of the inputs of the ...
0
votes
0answers
87 views

Importance Sampling of 2D constant piecewise function convertible to 1D?

So I have a constant piecewise 2D function (luminance values of pixels of an image) that I am writing an importance sampling algorithm for. I was going to write my algorithm by first sampling the 1D ...
0
votes
0answers
11 views

LDA with fixed topics?

Suppose I have a collection of "topic" probability distributions $\{\phi^{z}\}$ for LDA (Latent Dirichlet Allocation) that I have found via alternate methods; is there a closed form MLE for the ...
0
votes
1answer
55 views

Double integral estimation with Monte Carlo method by importance

I have this integral $\int_{0}^{\infty}\int_{0}^{\infty} e^{-x-y-xy}dxdy$ I have to estimate the value of it using the Monte Carlo method, using as importance function the exponential PDF. Does ...
0
votes
0answers
37 views

Monte Carlo simulations, accuracy of mean vs variance of answers

I am working on a Monte Carlo simulation where two inputs are being used, $N$ is the amount of simulations I use, and $M$ controls the detail of each simulation (specifically the amount of time step ...
0
votes
1answer
37 views

Looking for Method to evaluate the optimal node rate vs number of simulation rate in a Monte Carlo simulation

I am currently working on evaluating an American Option using a Monte Carlo simulation, and I am getting answers but they vary quite a bit. The two variables that I can alter are number of simulations ...
2
votes
0answers
44 views

Mixing time for metropolis chain on graph coloring

I'm reading the Markov Chains and Mixing Times by David Levin et al.. In section 5.4 page 71 a proof is given for a bound of mixing time for the Metropolis Chain on graph coloring. In the proof, such ...
0
votes
0answers
13 views

Additive model and simulation of the time series

It will be great if somebody could help me with following task.. Im trying to solve it about 3 days and read a lot of papers according this issue, but it doesnt work My task is to simulate a time ...
0
votes
1answer
17 views

A question about importance sample and Metropolis Algorithm

I am reading this paper by Beichl, I., & Sullivan, F. (2000) on Metropolis algorithm. I understand rejection sample. In the section "The Rejection Sample", I can understand the equation: ...
0
votes
0answers
17 views

How to choose a proper numerical optimisation method

Given a problem in numerical analysis in finance/econometrics, how to decide whther to choose Monte Carlo, Newton Raphson , Finite Difference , Gradient descent? I had this silly misconception that ...
0
votes
2answers
143 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: ...
3
votes
1answer
76 views

Applying MCMC Metropolis algorithm

I'm interested in all possible paths (on the grid $\mathbb{N}^2 $) that goes from $ (0,0) $ to $ (n, n) $. At each step there are two possibilities: go right or go up. The path is a sequence $ ...
0
votes
1answer
45 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
143 views

Are irreducible, positiv-definite Markov chains aperiodic?

If $M$ is the transition matrix of a discrete Markov chain, and $M$ is both irreducible, symmetric and positiv-definite, is the resulting Markov chain necessarily aperiodic? In my intuition, ...
0
votes
0answers
18 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
15 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
14 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 ...
2
votes
1answer
26 views

Does the perimeter of a 2-D object “count” toward its area?

I'm writing a quick Monte-Carlo simulation for a class in Matlab in order to estimate the value of pi as demonstrated in this gif: http://en.wikipedia.org/wiki/File:Pi_30K.gif However, I'm not sure ...
1
vote
0answers
26 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 ...
2
votes
0answers
11 views

Monte Carlo Markov Chain Simulation Issues

The Markov Chain is uniformly distributed across all $50$x$50$ matrices of entries $0$ and $$1 with no neighboring $1's$. I am supposed to run a MC simulation to check the probability that the ...
0
votes
0answers
33 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 ...
-1
votes
0answers
18 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. ...
1
vote
1answer
47 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
21 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
33 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 ...