Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [simulations]

The tag has no usage guidance.

1
vote
1answer
51 views

Monte Carlo simulations in Python using quasi random standard normal numbers using sobol sequences gives erroneous values

I am trying to perform Monte Carlo Simulations using quasi random standard normal numbers. I understand that we can use sobol sequences to generate uniform numbers, and then use probability integral ...
0
votes
1answer
42 views

Simulating stock prices with and without intermediate paths

So I am simulating stock prices with what I believe to be geometric Brownian motion using parameters from the usual Black-Scholes framework (Please correct me if I am wrong) with the following formula:...
1
vote
0answers
31 views

Negative theta in Log-linear stochastic volatility model

I was asked to simulate the following geometric Brownian motion to get paths for the SPX stock price. the process follows a Log-Linear stochastic volatility. $dS_t = \mu S_tdt+e^VS_tdW_1 $ where ...
3
votes
1answer
209 views

Why is it more accurate to simulate ln(S) rather than S?

Let's take a process $S$ that satisfies: \begin{equation} dS = \mu S dt + \sigma S dz \end{equation} with $dz$ a Wiener process, $\sigma$ the volatility of $S$, $\mu$ the expected return of $S$. From ...
2
votes
0answers
78 views

Longstaff-Schwartz, special american option simulation using Python (numpy package)

I got a put option, which can be exercised 3 times, all at different times, which are each month of a year $$t_1 = \frac{1}{12}, t_2 = \frac{2}{12} ... t_{12} = 1$$. Respectively, if exercised at $$...
0
votes
1answer
71 views

Simulated Sharpe Ratio Calculation for Leveraged Portfolio

I've written some VBA code to simulate the effect of borrowing money, investing it, and repaying the loan daily. PseduoCode: Start with a portfolio value of P = 1 Each day borrow P, invest 2*P, ...
1
vote
1answer
126 views

Shifted Log-Normal model

I am trying to understand how the shifted log-normal model works, in which we shift a log-normal model by a factor before the simulation so that interest rates don't turn negative during the ...
1
vote
0answers
26 views

Why no prepayment fee for the reverse mortgage?

I am currently studying the costs (to lender) of adding certain additional options to the reverse mortgage, including the option of prepayment. Would there be any scenarios of housing price/mortgage ...
0
votes
1answer
117 views

negative values in geometric brownian motion

A GBM $ \frac{dx}{x} = \mu dx + \sigma dW $ solves to $x_t = x_o e^{(\mu - \sigma^2)t + \sigma W_t}$ From the solution, it is clear that $x_t$ cannot become negative. However, it is not so clear ...
0
votes
1answer
57 views

Single vs Multi factor interest rate model

How do we explain the difference beween a single and multi factor interest rate model. Short term interest rate is one of the factor which is used in drift and vol calculation but what are other ...
2
votes
1answer
101 views

Merton's Jump diffusion model: Specify poisson rate

Currently applying the Merton's jump diffusion to test how Option price change as parameters change. However, I am struggling to specify the poisson rate $\lambda$. We know that: $P(\text{There is a ...
1
vote
0answers
76 views

Multivariate Hawkes Process Simulation

I am trying to implement Ogata's thinning algorithm to simulate multivariate Hawkes Processes in Python (the algorithm can be found here: https://www.math.fsu.edu/~ychen/research/Thinning%20algorithm....
3
votes
0answers
59 views

Solving BSDE in R

I was wondering how to implement a BSDE approximation in R. For example, if I have the toy BSDE $$ dX_t = \mu dt + \sigma dW_t ; X_T\sim N(\mu_1,\sigma_1), $$ for fixed real numbers $\mu,\mu_1,\sigma,...
1
vote
1answer
97 views

Simulating assets of different currencies

I have a situation as follows: One year call option on a Euro stock with a Euro denominated strike. Knock in feature as follows - The option can only pay out if the growth in the Euro stock over ...
0
votes
1answer
38 views

In CVA simulation, timesteps vs number of simulations?

On a CVA system with limited computational power. For pricing, What is best, More timesteps and less number of simulations or less timesteps and more number of simulations? for example with a whole ...
1
vote
0answers
22 views

Simulating Taxed Equity Return Series (U.S.)

I'm looking to learn how to correctly simulate taxes on dividends and capital gains on simulated return series for U.S. Equities with dividend reinvestment. I understand I will have to keep track of ...
1
vote
2answers
94 views

Distribution of data for GBM

I am running some Monte Carlo simulations with GBM on time series of commodity prices. First of all, the price data is annual between 1900-1950. I would firstly like to know if it is bad practice to ...
0
votes
1answer
71 views

Simulating trading strategies

I want to simulate real time trading strategies. For simplicity, let's say I only want to simulate a long-only portfolio on S&P500. I have a couple of questions: Is there a place online where ...
0
votes
1answer
157 views

How to simulate this Gamma expansion in a Python way

Here is the simulation that I want to do: For each of the 10 million simulation paths, I have n = 100 lambda values in sequence (the lambda vector is the same for all paths), Using each of the lambda ...
2
votes
0answers
62 views

Generating Correlated Quasi Random Numbers

Hi I am trying to generate correlated quasi random numbers using a sobol sequence in matlab. My Problem is the Following: Using "standard" random numbers it is easy to generate the 6 correlated random ...
2
votes
1answer
147 views

How to model High/Low prices for Stocks with Monte Carlo

I'm using monte carlo simulation to model stock paths and measure risk, but I was wondering if there is a way to simulate the full bar/candle chart with open, high, low and close prices , as I'm only ...
2
votes
1answer
139 views

How to interpret and define statistics of GBM output

I am trying to model the future prices of a number of commodities. For this, I am applying geometric Brownian motion, writing a Monte Carlo code in Python. Given that I want to estimate tommorows ...
1
vote
0answers
82 views

Monte Carlo Simulation of correlated returns based on different frequencies

I am simulating through Monte Carlo, multivariate correlated returns of different products composing an Oil&Gas portfolio. The historical prices (from which I computed the log-returns) of the ...
4
votes
1answer
591 views

Monte Carlo - Multivariate Simulation of Returns

I am implementing a Monte Carlo simulation in R to generate multivariate correlated returns. In doing this I have used the Cholesky decomposition, applied to the covariance matrix. However, I saw that ...
1
vote
0answers
35 views

Bootstrapping to Judge the Fit of a Sampled Return Distribution

Consider the following: I have sampled yearly stock returns from a specified distribution. What I want to do is compare how well my sampled distribution fits the empirical distribution of yearly ...
1
vote
0answers
77 views

Implementing Pykthin Multi-factor adjustment

I've made a mistake in the implementation of Pykthin Multi-factor adjustment which I'm fairly certain comes from me not understanding the model completely. The model was developed to drastically ...
4
votes
1answer
144 views

Simulate double exponential process with correlated jumps?

So, I'm trying to simulate a correlated double exponential jump process for two assets, and I understand the pure exponential jump process ($\eta_1$ and $\eta_2$, the probability of an upward jump ...
1
vote
0answers
162 views

Open Source library for calculating exposures?

I would like to know an open source quantitative library/ies that can calculate exposures out of the box (I have investigated a bit on OpenGamma/Strata libraries with no luck and the website of ...
1
vote
2answers
866 views

Simulating a path of bond yields by Monte Carlo (Python)

I have a number of given time series for bond yields (given in a dataframe in pandas package in Python). I need to do the following task in Python: "1. Simulate 1000 path 30 steps ahead for any yield ...
1
vote
1answer
235 views

How to sample from a copula in matlab

I have two random variables (say, X and Y). Each of these rv's are defined by their CDFs (CDF_X and CDF_Y). These CDFs were obtained empirically, so they are a "stair" graph. I also have a copula C ...
1
vote
1answer
256 views

Terminal Variance in the Heston Model

I am trying to understand the basics of financial models. Random Walk as a model for asset prices. We use gaussian random numbers to generate a Gaussian Random walk. The variance of the terminal ...
0
votes
1answer
96 views

time step choice impact in Vasicek model simulations

I am trying to make some computations using Vasicek short rate model. Especially I a trying to compare exact expectation(obtained with the formula) and the expectation from Monte Carlo simulation. ...
4
votes
0answers
60 views

Equity Options - “How do I build a forward simulation model with regards to shocks in spot pricing and IV?”

I am trying to build a "What-If" Portfolio, consisting of a total of 20 options, across different tenors, strikes (delta), but on the same security. Simply put, the objective is for me to test the ...
1
vote
1answer
336 views

Books about Monte Carlo Simulation on derivatives with Python

I am looking for a good reference for Monte Carlo simulation applied to derivatives with Python. Most books I found until now deal with C++... I have found "Derivatives Analytics with Python" by Yves ...
3
votes
1answer
602 views

How are Brownian Bridges used in derivatives pricing in practice?

A similar question has already been asked in the past, unfortunately the 2nd question of the OP was never really addressed. Most references found on internet on Brownian Bridge and Monte-Carlo ...
1
vote
2answers
278 views

Simulation curves; PRIIPS category 3

Once the yield matrix has been computed, the eigenvectors must be calculated to project the yield matrix on the 3 main dimensions. Tehen is wasted to calculate the yield matrix to be used for the ...
1
vote
2answers
115 views

Choosing a proxy for asset credit event correlations

I'm interested in modeling the joint likelihood for rating changes and default events across a portfolio of bonds. To estimate the correlation between these assets, I can use a third-party factor ...
3
votes
1answer
165 views

Projecting a Thiele differential equation with Black Scholes returns

I am trying to solve the equation $\frac{d}{dt}V(t)=r(t)V(t)+\pi-\mu(x+t)(b_d-V(t))$ numerically using the R function 'ode'. This is a Thiele differential equation for a life insurance reserve with ...
1
vote
0answers
97 views

Simulating asset returns: (Academia) state of the art

I want to run some simulation studies of (linear) factor models and for that reasons I am wondering about the features such a simulation should contain - every suggestion is welcome, I'll do my best ...
0
votes
2answers
224 views

Log normal price simulation

I'm trying to figure out a spreadsheet I have which simulates 50000 returns in excel using the following function: LOGNORM.INV(RAND(),0,0.35)-1 Question: How ...
0
votes
1answer
344 views

Potential Future Exposure (PFE): Is there any Rigorous Walk Through with Data?

I have searched on the Internet and in several books (including John C. Hull and Jon Gregory) for concrete examples of Potential Future Exposure (PFE), but haven't had any success so far. I would ...
1
vote
0answers
88 views

Simulation of Traders [closed]

I am writing a simulation of how traders will behave in an emergent market. The idea is to see how traders can use information from other traders to make a decisions as to whether they will buy, sell, ...
3
votes
0answers
87 views

Electricity Prices: Change of measure in practice

I'm working on a model of electricity prices. I have empirical data $X(t)$ and managed to find a reasonable fit given by a Levy process $\hat{X}(t)$. I understand in theory what a risk-neutral ...
3
votes
1answer
182 views

How to price the American style Asian option with recent N day average

How to price the American style Asian option with recent N day average, for example, we exercise at t day, then the payment is $$...
2
votes
1answer
78 views

A bug in delta hedging, when for a certain step dS=0

Suppose we are doing a delta hedging simulation according to Black Scholes, where the initial condition are [stockPrice, strike, timeToExpire ,riskFreeRate, dividend, sigma, isCall] = [100, 100, 1, 0, ...
1
vote
1answer
52 views

Extreme cases of normal random numbers and NaN

While trying to implement my version of Euler's method for simulating a SDE in C++, I came up with a problem. It occurs in some cases that the path generated by my method ends up giving values which ...
-1
votes
0answers
58 views

Methods for modelling price shocks

I am doing stress-testing of central counterparties, how a price shock affects them. The central counterparties calculate the required collateral based on a "normal" market, the distribution of the ...
1
vote
0answers
34 views

Simulating Asset Prices by Independently Simulating Supply and Demand

If I have an asset, whose supply is generally mean-reverting and whose demand is generally cyclical, could I somehow simulate / project the supply and demand levels across multiple discrete time ...
1
vote
1answer
237 views

In a Monte Carlo simulation, will a delta hedge control variate necessarily reduce the standard error more than an antithetic variate?

I have four Monte Carlo simulations and will list them in order of highest standard error to lowest. Plain MC MC with delta hedge control variate MC with antithetic variate MC with antithetic and ...
2
votes
0answers
98 views

Smooth ornstein uhlenbeck process

I want to simulate paths for a commodity price. I use the historic data in the following way: $X_t$ is the price. $\ln\left(\frac{X_t}{X_{t-1}}\right)$ is the daily return. I calculate the slope of ...