- Today, 3:29 pm
- Forum: Numerical Methods Forum
- Topic: Smoothing splines (clamped spline)
- Replies:
**30** - Views:
**1601**

The data sample is not draws from a pdf, but option prices. However you approach the problem, you have to tackle the issue of extrapolation (of the risk-neutral pdf) to non-marketed strikes. Good point. Does this mean that kernel functions are not 'suitable' for data with no discernible underlying ...

- Yesterday, 6:46 pm
- Forum: Trading Forum
- Topic: Early Exercise Risk Management for American Option
- Replies:
**6** - Views:
**290**

Fair point. The OP should indicate if he/she is talking about a negative interest rate environment. (I indeed was not). If so, there is the possibility that the "non-optimal" early exercises being seen are actually optimal. Perhaps bearish will elaborate ...

- September 20th, 2019, 2:46 pm
- Forum: Numerical Methods Forum
- Topic: Smoothing splines (clamped spline)
- Replies:
**30** - Views:
**1601**

The data sample is not draws from a pdf, but option prices. However you approach the problem, you have to tackle the issue of extrapolation (of the risk-neutral pdf) to non-marketed strikes. Anyway, I am close to finishing a write-up for a larger problem (which contains this risk-neutral density est...

- September 20th, 2019, 2:32 pm
- Forum: Trading Forum
- Topic: Early Exercise Risk Management for American Option
- Replies:
**6** - Views:
**290**

Consider covered call writing: buy the stock, sell a call. There will be a maximum return if the call is optimally exercised. If the call is exercised early (non-optimally), that maximum return will be obtained *sooner*, which is a good thing. To take a trivial case, suppose there are no dividends,...

- September 16th, 2019, 9:22 pm
- Forum: Numerical Methods Forum
- Topic: Smoothing splines (clamped spline)
- Replies:
**30** - Views:
**1601**

The problem is finding a risk neutral density from option quotes. The spline approach which I abandoned was: fit a smoothing spline to option smile data, with given derivative endpoints. The purpose of the derivatives, as opposed to so-called 'natural' boundary conditions, was to be improved extrap...

- September 16th, 2019, 5:53 pm
- Forum: Numerical Methods Forum
- Topic: Smoothing splines (clamped spline)
- Replies:
**30** - Views:
**1601**

Daniel, Thanks for the book link. Sorry, don't have a good numerical example because I abandoned this approach for another. But, may get back to it.

- September 8th, 2019, 5:45 pm
- Forum: Book And Research Paper Forum
- Topic: Performance of my option pricing model
- Replies:
**5** - Views:
**1181**

Just thinking out loud, if you *don't* want the vol to explode maybe there are some modifications that don't completely wreck the nice properties. Either truncating the lattice or keeping the full lattice but somehow truncating just the volatility. Or, a rectangular grid version of your model with e...

- September 8th, 2019, 2:27 pm
- Forum: Book And Research Paper Forum
- Topic: Performance of my option pricing model
- Replies:
**5** - Views:
**1181**

I'll bite. Why is volatility exploding "a more realistic feature"?

- September 5th, 2019, 10:23 pm
- Forum: General Forum
- Topic: Negative Strike Price
- Replies:
**11** - Views:
**1110**

So they way it works is say you have 2 different oil market locations. Spot1 is $55 and Spot2 is $50. Right now we have a -$5 basis differential and I bought a put or call option with a (negative)$5 strike on the basis diff. As long as both prices stay negative the model is fine but if your bas...

- August 25th, 2019, 6:17 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

Since I give a (quasi)-analytic "absorbing" solution, all the various cases are checkable. BTW, when differentiating the latter w.r.t. r (which you asked about before), the identity [$] M'(a,b,z) = \frac{a}{b} M(a+1,b+1,z)[$] might be helpful (from Abramowitz & Stegun). Using that will keep the ...

- August 25th, 2019, 2:01 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

Alan, Regarding eq. 10.4 of your book, how can it be differentiated WRT r and the other parameters? >> Standard differentiation: "r" appears in 3 places. BTW for [$]B_r[$] the Feller condition is always satisfied at [$]r = 0[$]. Correct? >> The Feller condition, applied to the [$]B_r[$] pde, and us...

- August 25th, 2019, 1:51 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

BTW is it not so that CIR has a unique analytical solution for all values of the parameter regime? No. As explained on pg. 456 and also on pg 349 of my Vol. II, the CIR pde does *not* have a unique solution when (to use your notation) [$]0 < a < \sigma^2/2[$]. This is where the boundary is classifi...

- August 22nd, 2019, 7:29 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

The PDE for this greek does indeed contain terms in the bond price. The far-field BC sounds reasonable. At the near field the Feller condition for the 'greek PDE' is always satisfied which means we have a drift PDE at that boundary. It must be solved numerically (upwinding) in the worst case. I get...

- August 22nd, 2019, 2:17 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

The nice thing about NDSolve (in Mathematica) is that it automatically adopts a high order inward pointing derivative at boundaries, which then couples the system smoothly to the behavior of the solution at the interior points.

[Erroneous second comment deleted]

[Erroneous second comment deleted]

- August 21st, 2019, 8:11 pm
- Forum: Numerical Methods Forum
- Topic: One-liner questions of a numerical kind
- Replies:
**40** - Views:
**2860**

Consider the well-known PDE for the bind price [$]B[$] in the CIR model.Then differentiate the PDE with respect to [$]r[$] to get a PDE in [$]\frac{\partial B}{\partial r}[$] (this is a bespoke CSE equation). Some points: 1. The Fichera/Feller conditions implies that no BC is given at [$]r= 0[$] a...

GZIP: On