The tag has no wiki summary.

learn more… | top users | synonyms

2
votes
1answer
110 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 ...
3
votes
0answers
51 views

Maximize or find an upper bound of the function $kx^{k-1}\exp(-\mu(x^k-x))$

I was programming some random variable simulation using the acceptance-rejection method and I encounter with the Weibull$(k,\lambda)$ distribution. This random variable is posible to simulate with ...
3
votes
0answers
69 views

Classify knots in a closed bead-spring like polymer simulation

my problem is to detect the crossing number (or another knot invariant) of a simulated polymer. A polymer is a closed bead spring, which mean that it is represented by a set of points connected by ...
3
votes
0answers
182 views

Simulation and Stochastic Processes

Supposing we want to take a sample from the distribution $p(x)=cp^*(x)$ where $c$ is the normalization constant and $p^*(x)$ is given by $$p^*(x)=0.5\exp(-(x-\mu_1)^2)+0.5\exp(-(x-\mu_2)^2).$$ ...
2
votes
0answers
25 views

Are these two approaches to calculating return rate mathematically consistent?

I have coded two C# programs, which use two different approaches to evaluate the outcome of a certain casino-style game (casino-style in the sense that the user pays points to take a turn, and ...
2
votes
0answers
21 views

orderstatistics of uniform distributions on different ranges

During a simulation I discovered an interesting phenomenon: Given you have 3 agents. 2 are uniformly distributed between [0,1] and one between [0,2]. The question is how often do the smaller agents ...
2
votes
0answers
113 views

Positive eigenvalues in differential-algebraic equations not appearing in time-domain simulation

I am solving a system of equations derived from power system applications. It consists of index-1 differential and algebraic equations in the form: $$\dot{x}=f(x,y) \\ 0=g(x,y)$$ To get the ...
2
votes
0answers
133 views

Simulating first passage times

I have a Brownian motion $X_t = \mu t+\sigma W_t$, where $W_t$ standard Brownian motion. I know that the first passage time $\tau = \min\{t|X_t\leq\alpha\}$, is Inverse Gaussian distributed i.e., ...
1
vote
0answers
25 views

Simulating r.v.'s from a joint density by rejection sampling in R. Continued

I wish to sample variables $v$ and $w$ from the joint density $$(v+w)e^{-\frac{(v+w)^{2}}{2x_{0}}-2\mu v-(\mu -\lambda )w},$$ where $x_0$, $\mu$ and $\lambda$ can be seen as positive constant. Since ...
1
vote
0answers
38 views

Personal Experiences with Probability Simulation

Simulations methods are increasingly used in theoretical and (especially) applied probability. Personally, I have used simulation for purposes that range from recreational Q&A to applications of ...
1
vote
0answers
29 views

Is Principal Component Analysis applicable to this type of situation?

I'm trying to model the response of ant populations to pheromones in this way: The ants are simulated as Self Propelled Particles with internal energy. They undergo acceleration due to their internal ...
1
vote
0answers
54 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 ...
1
vote
0answers
56 views

Statistics Marble Simulation

Given the jar containing 20 red, 5 white, and 10 blue marbles, Joey thinks a more interesting problem would be to find the number of marbles you would have to draw, without replacing them in the jar, ...
1
vote
0answers
26 views

Can you simulate from a cantor distribution?

Has someone run across a method for generating random variates from a Cantor Distribution? It seems like its abstract definition prevents this. In essence, can one "invert" the Cantor Function?
1
vote
0answers
79 views

Backward Euler method with a cross-product.

I want to solve the following differential equation with the backward Euler method ...
1
vote
0answers
20 views

Stability of simulation of brownian noise

As I understand, Brownian noise can be simulated by the process $$x_{n+1}=x_n+R_n$$ where $R\sim U[-a,a]$. The expected value for $x_n$ is then $x_0$. But $\text{Var} x_n\to\infty$ as $n\to\infty$ ...
1
vote
0answers
35 views

Distribution of phone calls during 24h

I would like to model the amount of phone calls at each time of the day. The phone calls should follow a poisson distribution and at 12:00 there should be the peak. So, semantically what I would like ...
1
vote
0answers
45 views

Forest fire simulation; analytically constructing a function for tree residual after fire

Consider a Cellular Automaton with an $n \times n$ grid, where each cell corresponds either to a tree or dirt. We assign a tree to cell $(i,j)$ by probability $p$. Next, we initiate a fire in some ...
1
vote
0answers
85 views

Kalman Filter application to non-linear system.

I want to use the Kalman filter to have a better estimate of the state of a system which I know its equations of motion: ...
1
vote
0answers
44 views

simulate hitting time of Brownian motion

let's say I have a brownian motion $W_t$, and I know the value of $W_1$. Is there a way to simulate the hitting time of $W_t$ and a given function $f(t)$ ? For instance I know that if $f(t)$ is a ...
1
vote
0answers
43 views

Random graphs simulation

Reading the article "Emergence of scaling in random network, by Barabasi and Albert" I faced a lot of results obtained by simulations of the A-B random graph model. I always wanted to do such ...
1
vote
0answers
57 views

Samples from the Dirichlet measure

In Ferguson, 1973, Definition 2, he defines a sample of size $n$ from a random probability measure $G$ on $(\mathcal{X}, \mathcal{B})$ as: $$ P(X_1 \in C_1, \cdots, X_n \in C_n | G(A_1), \cdots, ...
1
vote
0answers
31 views

Single evaluation for using exponential sampling until past a point

I am trying to improve an algorithm that looks like the following (and am getting stumped): I am provided with a starting time, rate, and a target time. I then use an exponential distribution to ...
1
vote
0answers
46 views

Simulating of GBM

I have a question regarding the simulation of a GBM. I have found similar questions here but nothing which takes reference to my specific problem: Given a GBM of the form $dS(t) = \mu S(t) dt + ...
1
vote
0answers
112 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
107 views

Maximizing noisy unknown function

I'm interested in maximizing a function $f(\mathbf \theta)$, where $\theta \in \mathbb R^p$. The problem is that I don't know the analytic form of the function, or of its derivatives. The only thing ...
1
vote
0answers
182 views

Expected value involving a conditional multinomial distribution

$X|z$ has a multinomial distribution $MN(m, \mathbf{q}(z))$ where $z$ is either 0 or 1 with probability $1/2$. I need to find: $E_X[\max\{\Pr(z=1|X), \Pr(z=0|X\}]$. Is there an analytical form to ...
1
vote
0answers
263 views

Semi implicit integration - stability issues

I am trying to decide whether to use semi-implicit integration vs. explicit integration (particularly Position Verlet over Semi implicit Euler). Although the Verlet approach is widely used and is ...
1
vote
0answers
181 views

Optimization via Simulation

I want to minimize and objective function $\hat{B_i}$ $i\in l$, which can be computed by a matlab code (assume $\operatorname{findB}(a, b, c)$ returns $B$. I have the following optimization problem: ...
1
vote
0answers
163 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), ...
0
votes
0answers
8 views

Stability of boostrap confidence intervals

As a word of background, I want to show that certain result is stable when averaging over a large number of simulations, but could be just a lucky draw with a small number of simulations. I have a ...
0
votes
0answers
14 views

Numerical Triple integral with three other parameters in R

I am trying to integrate this function $f(u,v,w; t,x_{0},z)$ with respect to three variables, $u$, $v$, $w$, although the function also have other three parameters $t$, $x_0$, and $z$. Question: How ...
0
votes
0answers
18 views

Sinusoid with zero boundary conditions on $[0,1] \times [0,1]$ grid

I want to make a sinusoidal plot (any shape is welcome) on a $[0,1] \times [0,1]$ grid, with boundary conditions equal to zero. It should resemble a membrane fixed along its edge. I tried out some ...
0
votes
0answers
26 views

Simulation of an ODE model with non-constant parameter

I have a model, I can formulate the model using ordinary differential equation with parameter $P$. I want to simulate the model, but instead of using a fixed constant $P$ for the parameter, I want to ...
0
votes
0answers
8 views

Simulation model, reality check on sample data

I wrote a simple simulation in C (a next-event simulation) and I'd like an opinion on the consistence of the sample data. The model is that of a two nodes service, the user enter on the first one and ...
0
votes
0answers
12 views

How can I efficiently optimize stochastic multi-modal functions?

I'm looking for methods for optimizing stochastic functions. I'm probably abusing the notation here, since this is a new field for me. By stochastic functions I mean functions whose output is a ...
0
votes
0answers
24 views

Strictly local martingales: what is the intuition behind them?

I did post this on the Quant Finance exchange a while back, but without any luck A process $X_t$ is a local martingale if for each increasing sequence of stopping times $\{τ_k,\ k=1,2,...\}$ the ...
0
votes
0answers
34 views

Understanding simulation of Brownian Motion

I am trying to understand the simulation of Brownian Motion given at http://www.math.uah.edu/stat/applets/BrownianMotion.html. There are four boxes in this simulation. For the purpose of this question ...
0
votes
0answers
19 views

Skew in black scholes model

We are modeling Foreign exchange rates using Black Scholes model given below: $F_t = F_{t-1} + (r_d - r_f)F_{t-1} dt + \sigma\cdot F_{t-1}\cdot dW_t$ Where $F_t$ and $F_{t-1}$ are FX rates at time ...
0
votes
0answers
45 views

Uniform Sampling over Convex Polytope (not full-dimensional)

I want to simulate a uniform distribution on a convex polytope that is not full-dimensional for optimization purposes (to generate random points on the set I want to minimize over). The polytope is ...
0
votes
0answers
22 views

Sampling from a random distribution with a margin of error

A survey of a population is taken using sampling. It is determined that 70% prefer option A and 30% prefer option B with a margin of error being 5%. Typically when simulating the process with the ...
0
votes
0answers
54 views

Transform this IVP into a first-order autonomous system with $\dot{x} = Ax + b$

I'm fairly new to differential equations and need some guidance with this following problem: Transform the following IVP into a first-order autonomous system in the form $\dot{x} = Ax + b$: $$ ...
0
votes
0answers
41 views

Expected Probability of a Random Agent and a Probabilistic Agent

I'm running simulations on two agents: random agent and probabilistic agent. The world they are running in is the Wumpus World where the agent is dropped in a 4x4 grid where each cell has a 20% chance ...
0
votes
0answers
16 views

Individual particle tracking simulation

I want to do a simulation of a stochastic system. I have 4 types of cell, each will divide or die with a certain probability. Let's say : A-> A+A A -> A+B A -> A+C B-> B+B B-> die and so on... ...
0
votes
0answers
74 views

Water swallowing in Matlab

I want to simulate some water passing through a vertical cylinder in Matlab, and I would like to implement a 3d animation of it. I built the cylinder using the patch function, but I do not know how to ...
0
votes
0answers
23 views

Simulation Lévy process

I need to simulate a Lévy process from its characteristic triple $(\gamma,\Sigma,\nu)$ where $\nu$ is the Lévy measure. I know that I can simulate it by summing a brownian motion and a compound ...
0
votes
0answers
99 views

Estimation covariance of the Kalman filter state

I implemented Kalman filtering for a simplest 1D coordinate+velocity model. The prediction worked, but I wanted to estimate the prediction probability distribution. I.e. how likely it is that the ...
0
votes
0answers
12 views

Generating a given length sample and skewness whose normality is verified by one normality test but not by an other

I just want to generate 1 sample of length$=n>30$, |skewness$=S|<0.3$ and for which normality is not rejected by Shapiro wilk test of normality but rejected by Anderson darling test of ...
0
votes
0answers
15 views

spatial-partitioning based physical simulation

I've learnt that spatial-partitioning based physical simulation is kind of "approximate" computation. Is it because: since the whole space is partitioned into cells, and only the interactions of ...
0
votes
0answers
42 views

STINT Approximate stochastic integrals

This is a matlab code to simulate stochastic integrals: ...