When did I ever say it was "dirty"? Anyways, it's not exactly rocket science right here: pp = spline(strikes, prices); ppd = diffPP(pp); ppdd = diffPP(ppd); splitIncr = strikeStep / 2; splitStrikes = [strikes(1)-splitIncr:splitIncr:strikes(end)+splitIncr]'; splitPrices = ppval(ppdd, splitStrikes); for j = 1:length(strikes)-1 RNProbs(j) = trapz(splitStrikes(j*2-1:j*2+1),splitPrices(j*2-1:j*2+1)); end This should work, if I haven't screwed anything up. Now I hope I don't go to jail for violating that crazy patent.
It's useless to say that it needs to be rescaled to get the sum of probabilities equal to unity. Masteratwork
Certainly, "needless to say", but definitely not "useless to say", MAW. You're absolutely correct, of course. I left out a few steps required to massage the results (such as discarding the insignificant points, normalizing, etc), as I thought they should be relatively self-explanatory.
Sorry mixed up the wording; I wanted to make sure I was on the same page....I can usually see what someone is talking about much better by just seeing what they are trying to do in code!
I still dont understand the difference guys, I understood that the distribution curve is always the same but according to the suggestions I guess it is not. I understand the implied deviation but not the implied distribution, can somebody please explain, thx.
Implied risk neutral distribution is the risk neutral distribution extracted from market option prices. This risk neutral probability distribution, that is not the real one which can't be known, is just a tool to simplify pricing stuff for exotic options for example. If there are different implied volties, that means that log returns are not normally distributed, and risk neutral distribution would exhibit skew and kurtosis. Because options are traded on a market, implied volties move hence there're implied expected deviation and implied risk neutral distribution variations. Masteratwork
I managed to draw the annualized lognormal distribution curve of one strike with underlying at 351 with 0.20 IV, see attachment but I am still not shure how to get something like this http://www.elitetrader.com/vb/attachment.php?s=&postid=2668308 wich I assume represents the full chain and is what we are talking about, how would exp(-rT)(C1-2C2+C3)/d² lead to such? Could you give an example maybe to clear it up a bit.