A vast area which includes generating results from computer models.

learn more… | top users | synonyms

0
votes
0answers
9 views

A Single-Server Queueing System

Consider a single-server queueing model in which customers arrive according to a nonhomogeneous Poisson Process. Upon arriving they either enter service if the server if free or else they join the ...
3
votes
3answers
29 views

Simulation of interacting Ornstein-Uhlenbeck processes

I would like to simulate the following system of interacting OU processes on $[0,T]$: $$dX_t^1=(X_t^2-X_t^1)\,dt+\sigma_1 \,dW_t^1,\quad X_0^1=x_1$$ $$dX_t^2=(X_t^1-X_t^2)\,dt+\sigma_2 ...
0
votes
1answer
33 views

Simulating Random Vectors Based on Conditioning

I'm working on a project where I need to simulate random vectors $(Y, X_1,\dots,X_n)$ in order to understand the joint distribution $f(y,x_1,\dots,x_n)$. I wish to simulate enough random vectors so ...
0
votes
0answers
11 views

Customers and Anti-Customer Queueing Problem: What is the Customer delete probability

Hello may ask for your help? First the setting: I have got a problem with some queueing theory. The whole problem would be a grid of nodes, all nodes have an operation intensity $\mu_{i,j}$. ...
1
vote
2answers
3k views

Please help me solve this exponential distribution problem

Question 1 : The time to service a customer at a bank teller's counter is exponentially distributed with mean of 60 seconds. What is the probability that the three customers in the front of an ...
1
vote
1answer
436 views

Simulation of a Gaussian process on $R^2$ with a stationary kernel using the Karhunen-Loève expansion

Assume $X(\omega, t) \sim \mathcal{N}(0, K(\cdot, \cdot))$ is a real-valued, centered Gaussian process on $R^2$, i.e., $X: \Omega \times R^2 \to R$. Let the covariance function of the process be ...
-1
votes
0answers
18 views

Calibration geometric Brownian motion

I have a time series of monthly data that I want to calibrate GBM parameters μ and σ to. Using the distribution of the logreturns: $\ln\left(\dfrac{S(t_{i+1})}{S(t_{i})}\right) \sim N ...
2
votes
0answers
40 views

Solving Kepler's Equation

I've been working on simulating orbits. I've found that, when solving Kepler's equation, $M = E - \varepsilon\sin{E}$, I'm unsure about the solution to use. For a true anomaly $< \pi$, using the ...
-3
votes
0answers
17 views

How to calculate IMU data using position and orientation?

Can you please explain how can I simulate the IMU(Inertial Measurement Unit) data using position and orientation data. am given: p0 - starting position, q0 - starting orientation p1 - final ...
2
votes
4answers
81 views

What is the statistical equilibrium for this simulation of happy bubbles?

Happy Bubbles I hope this is not too specific or practical, but I just made a simulation of sorts and seem to have hit quite close to an equilibrium (by accident). Now I am wondering if and how you ...
0
votes
0answers
20 views

Exlain the significance of the uniform random variable for the simulation of random variables

I can think of the "Universality of the Uniform": Given an Unif(0,1) r.v., we can construct an r.v. with any cts distribution we want. Conversely, given an r.v. with an arbitrary cts ...
0
votes
0answers
115 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 ...
1
vote
0answers
17 views

How do we derive efficiency from robustness in the virtual ant solution to the traveling salesman problem?

Using virtual ants/swarm intelligence to solve the Traveling Salesman Problem is an example of using a robust system to solve an efficiency problem. We normally think of robustness and efficiency as ...
3
votes
0answers
21 views

Simulating a SDE

I have a question about simulating a SDE. I want to simulate $dS=\alpha(K-S)dt+\sigma S dZ$ with use of a Euler-marayama scheme. The numerical scheme becomes: $S_{i+1}=S_{i}+\alpha(K-S_{i})dt+\sigma ...
4
votes
2answers
121 views

Implementing Ornstein–Uhlenbeck in Matlab

I am reading this article on Wikipedia, where three sample paths of different OU-processes are plotted. I would like to do the same to learn how this works, but I face troubles implementing it in ...
0
votes
1answer
22 views

Calculation of Confidence Interval for Simulation

I wrote a simulation that does 1,000,000 trials. It returns 430,200 successes. I want to calculate a probability of success with confidence interval for that probability estimate. I use $\sigma = ...
1
vote
0answers
81 views

distribution of the length for a random walk on an infinite 2D grid

In connection with the flatland paradox, consider a 2D-random walk $(X_n)$ on $\mathbb{Z}^2$: the four moves of length one to W,E,N, and S are equaly likely at each time. For a fixed number of moves ...
0
votes
0answers
6 views

proposal distribution for metropolis algorithm

All, I'm wondering whether it is possible to use an asymmetric distribution, eg the exponential distribution as the proposal dist'n for a metropolis algorithm (wiki) (not the metropolis-hastings). ...
0
votes
1answer
12 views

Metropolis-Hastings Algorithm Clarification

All- Could you please clarify: from wikipedia, step two states at the end if the candidate is rejected, set xt+1 = xt, instead. I don't quite understand this, so you will have two of the exact same ...
1
vote
1answer
55 views

How can I solve numerically this partial differential equation?

I am reading this paper and came across a system of differential equations with 4 ODEs and 1 PDE. The PDE is given below. My question is how to solve this numerically in MATLAB , Python or ...
4
votes
0answers
127 views

How to compute this triple integral?

I am trying to do this triple integral $$\int_{0}^{\infty }\int_{0}^{\infty }\int_{0}^{\infty }(u+w)e^{-\frac{(u+w)^2}{2}}(v+w)e^{-\frac{(v+w)^2}{2}}(u+v)e^{-\frac{(u+v)^2}{2}}e^{-(\mu +\lambda ...
0
votes
0answers
41 views

Numerical Triple integral with three other parameters in R

I am trying to do this triple integral $$\int_{0}^{\infty }\int_{0}^{\infty }\int_{0}^{\infty }(u+w)e^{-\frac{(u+w)^2}{2}}(v+w)e^{-\frac{(v+w)^2}{2}}(u+v)e^{-\frac{(u+v)^2}{2}}e^{-(\mu +\lambda ...
1
vote
1answer
52 views

How to generate integer random numbers that equal to another random number?

I am running a simulation in Excel, and need to generate a group of integer random numbers summing up to another random integer, how can I possibly do it? For instance I have an integer random number ...
0
votes
1answer
38 views

Time between arrivals Distribution

I am simulating a hair parlor queue with m number of queues and 3 different types of services (queues). I was doing the time between arrival with a uniform distribution with a min value and a max ...
1
vote
0answers
29 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
votes
1answer
40 views

Finding density and distribution functions [closed]

I have been trying to understand probability by attempting past paper question and I have been stuck on this question all day and night. I am not quite sure how to go about finding the functions ...
2
votes
1answer
36 views

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

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 ...
0
votes
0answers
13 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 ...
1
vote
3answers
115 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
2answers
42 views

exponential distribution with probability about texts

It is 9:00 p.m. The time until Joe receives his next text message has an exponential distribution with mean 5 minutes. A text has not arrived for 5 minutes. Find the probability that none will arrive ...
1
vote
0answers
60 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 ...
0
votes
0answers
19 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
1answer
24 views

Difference and confidence intervals

I performed a few series of simulation to evaluate values of two parameters. Let's say, the results can be presented like this: ...
5
votes
1answer
8k views

Generate Correlated Normal Random Variables

I know that for the $2$-dimensional case: given a correlation $\rho$ you can generate the first and second values, $ X_1 $ and $X_2$, from the standard normal distribution. Then from there make $X_3$ ...
1
vote
0answers
32 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 ...
0
votes
0answers
30 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
1answer
39 views

Is it possible to estimate $e$ based on $N$?

Consider a sequence of random numbers $u_1,\dots,u_n$ obtained from a continuous distribution $F$. Let $N$ be the first one that is greater than its immediate predecessor. In othe words, ...
2
votes
1answer
67 views

Simulate correlated $\chi^2$ distribution

I understand that when one have multiple independent variable that follows $N\sim(0,1)$, denoted as $A$ if we have a correlation matrix $R$, we can generate correlated variables $B$ that are normally ...
0
votes
1answer
46 views

mean and std deviation from little's law

I have a service node with a queue The arrival rate in this service node is exponential(0.4) The service time is exponential(0.2) Running a simulation i've calculated the average population of the ...
1
vote
0answers
79 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
68 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, ...
0
votes
0answers
11 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
14 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 ...
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?
2
votes
4answers
2k views

How do I simulate a simple pendulum?

I have the equation of motion of a simple pendulum as $$\frac{d^2\theta}{dt^2} + \frac{g}{l}\sin \theta = 0$$ It's a second order equation. I am trying to simulate it using a SDL library in C++. I ...
1
vote
2answers
41 views

How to efficiently simulate successes of several trials if probabilities are inhomogeneous

If I'm doing a simulation with $n$ trials, each with probability $p$, a quick way to select the successful trials is to choose a binomially distributed random number. Then randomly choose that many ...
0
votes
1answer
27 views

What platform is best for simulating a stochastic process on a graph/network?

I'm simulating a dynamic process which was so far done only on a lattice, and Matlab was quite sufficient for that. However, I can't seem to find a convenient way to model such a process on a graph ...
0
votes
1answer
39 views

Permutation algorithm for simulating random variables $X,Y,Z \in [0,1]: X+Y+Z = 1$ and $X,Y,Z \sim U(0,1)$

Edit: Sorry, I tend to jump back and forth between math notation and computer science notation....often to the chagrin of my more rigorous colleagues (and Math.SE folks ;-) Also, I accidentally ...
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 ...
0
votes
1answer
38 views

uniqueness of svd decomposition and its role in statistical analysis

let us consider following model according to following link http://www4.ncsu.edu/~ipsen/REU09/chapter4.pdf it says that : The singular values are unique, but the singular vector matrices are ...