tag:blogger.com,1999:blog-17213918846986372332008-08-07T16:57:08.023+02:00the Jon Egil Strand blogJon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.comBlogger7125tag:blogger.com,1999:blog-1721391884698637233.post-43841262637243019112008-08-07T16:34:00.007+02:002008-08-07T16:57:08.035+02:002008-08-07T16:57:08.035+02:00Introducing call tilts and put tilts on FStat<a href="http://www.fstat.no/stocks/tilts"><img style="display:block; margin:0px auto 10px; text-align:center;cursor:pointer; cursor:hand;" src="http://4.bp.blogspot.com/_wwvWMP6W7uQ/SJsIFM6O8XI/AAAAAAAAAFI/AQHdNOoMBBY/s320/call_and_put_tilt.png" border="0" alt=""id="BLOGGER_PHOTO_ID_5231784277578805618" /></a><br /><br />If all call-options are expensive and all put-options are cheap, does the market expect the stock to rise?<br /><br />Prediction markets are all the buzz (<a href="http://www.econtalk.org/archives/2007/05/sunstein_on_inf.html">Cass Sunstein</a> or <a href="http://www.econtalk.org/archives/2007/05/hanson_on_healt.html">Robin Hanson</a>), ranging from sports to macro indicators. Heck, even Christopher Cox, chairman of the U.S. Securities and Exchange Commission <a href="http://www.cfawebcasts.org/cpe/what.cfm?test_id=811">is discussing wheither to police them</a>. <br /><br />Can we learn something from the derivatives market?<br /><br />On <a href="http://www.fstat.no">FStat</a>, we compare the simulated price of an option to it's market price. The ratio of the two gives a tilt. If the market price is higher than the simulated price, the tilt is above 100%. If the market price is lower, the tilt is less than 100%. <br /><br />In the above picture, we can see the calls are expensive (i.e. market price > simulated price), but the put's are fairly valued. <br /><br />Wheither or not this might be a decent predictor remains to be seen, but it'll be fun to find out.Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-62053059608313258272008-07-11T13:34:00.012+02:002008-08-07T16:32:07.415+02:002008-08-07T16:32:07.415+02:00t-copula and various measures of dependance<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://bp0.blogger.com/_wwvWMP6W7uQ/SHdIZ9Dk54I/AAAAAAAAAEM/e_EZYG1GJn0/s1600-h/t_copula_and_spearmans_rho.png"><img style="float:left; margin:0 10px 10px 0;cursor:pointer; cursor:hand;" src="http://bp0.blogger.com/_wwvWMP6W7uQ/SHdIZ9Dk54I/AAAAAAAAAEM/e_EZYG1GJn0/s320/t_copula_and_spearmans_rho.png" border="0" alt=""id="BLOGGER_PHOTO_ID_5221721903682873218" /></a><br /><br />The t-copula can use <a href="http://en.wikipedia.org/wiki/Spearman%27s_rank_correlation_coefficient">Spearman's rho</a> to describe the dependance between the variables. To build the all important intuition a graph helps clear things up.<br /><br />One really sees how the copula takes into account the relationship from a very negative dependance (when one rises the other falls, rho < 0), via a neutral relationship (one rises, the other is not incluenced, rho=0), to a strong positive relationship (one rises, the other rises as well, rho > 0).<br /><br /><br /><br /><br /><br /><br /><br /><br />Try it out in <a href="http://www.r-project.org">R</a>:<br /><pre><br /> library("copula")<br /><br /> rho_list <- c(-0.9, -0.7, -0.5, -0.3, 0, 0.3, 0.5, 0.7, 0.9)<br /><br /> par(mfrow=c(3,3))<br /> for(i in 1:length(rho_list)){<br /> rho <- rho_list[i]<br /> t.cop <- tCopula(c(rho), dim=2, dispstr="ex", df=3)<br /> persp(t.cop, dcopula, col='lightblue', main=paste("rho", rho), <br /> phi=20, xlab='U', ylab='V', zlab='tCopula', n=30)<br /> }<br /></pre>Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-53278860960060026952008-07-08T15:47:00.005+02:002008-07-09T22:13:57.466+02:002008-07-09T22:13:57.466+02:00Some copulas<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://bp0.blogger.com/_wwvWMP6W7uQ/SHN5qbSq2lI/AAAAAAAAADk/n2iNhb-j8wg/s1600-h/tCopula_pair.png"><img style="margin: 0pt 10px 10px 0pt; float: left; cursor: pointer;" src="http://bp0.blogger.com/_wwvWMP6W7uQ/SHN5qbSq2lI/AAAAAAAAADk/n2iNhb-j8wg/s320/tCopula_pair.png" alt="" id="BLOGGER_PHOTO_ID_5220650162840066642" border="0" /></a><br /><br /><br /><br />The word 'copula' originates from the Latin noun for a "link or tie" that connects two different things. Words such as 'couple' or even 'copulate' stem from the same source.<br /><br />In Mathematics, copulas are used to couple two or more univariate distributions. The univariate distributions can be specified in isolation, and then coupled afterwards. For example can a copula join together a Pareto distribution and a Gaussian distribution, or two Student-t distributions with various degrees of freedom. This kind of flexibility is not found in traditional multivariate statistics. For example, in a mulitvariate Student-t distribution, all marginal distributions must share the same tail heavyness.<br /><br />More formally, Wolfram describes copulas as follows:<br /><blockquote><br />Let H be a two-dimensional distribution function with marginal distribution functions F and G. Then there exists a copula C such that H(x,y)=C(F(x),G(y)).<br /><br />Conversely, for any univariate distribution functions F and G and any copula C, the function H is a two-dimensional distribution function with marginals F and G. Furthermore, if F and G are continuous, then C is unique.<br /><br />Weisstein, Eric W. "Sklar's Theorem." From MathWorld--A Wolfram Web Resource. <a href="http://mathworld.wolfram.com/SklarsTheorem.html">http://mathworld.wolfram.com/SklarsTheorem.html</a><br /></blockquote><br /><br />There exists a lot of Copulas, among which the Gaussian and Student-t are the most popular. These two, so called elliptical copulas, are symmetrical in that they deliver equal probabilities for both high and low joint movements.<br /><br />One significant difference exists though. A Gaussian copula has <span style="font-weight: bold;">zero probability</span> for joint extreme movements, while a Student-t distribution has <span style="font-weight: bold;">non-zero probability</span> for joint extreme movements. In for example a financial application, if the value of one stock hits rock bottom, only the Student-t copula allows the other stock to plunge as well. Therefore, in risk management, where everything that can go wrong, will go wrong at the same time, the student-t copula should be prefered.Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-70859245920590923452008-06-07T21:11:00.000+02:002008-06-07T21:38:10.825+02:002008-06-07T21:38:10.825+02:00Risk and Asset Allocation - a gem of a book<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://www.amazon.com/Risk-Asset-Allocation-Springer-Finance/dp/3540222138/ref=pd_bbs_sr_1?ie=UTF8&s=books&qid=1212865819&sr=8-1"><img style="float:left; margin:0 10px 10px 0;cursor:pointer; cursor:hand;width: 200px;" src="http://ecx.images-amazon.com/images/I/41be3NdWGEL._SL500_BO2,204,203,200_PIsitb-dp-500-arrow,TopRight,45,-64_OU01_AA240_SH20_.jpg" border="0" alt="" /></a><br /><br />Attilio Meucci's book <a href="http://www.amazon.com/Risk-Asset-Allocation-Springer-Finance/dp/3540222138/ref=pd_bbs_sr_1?ie=UTF8&s=books&qid=1212865819&sr=8-1">Risk and Asset Allocation</a> is a true gem. Clear and practical without sacrificing quality or detail. Meucci really sets the bar for scientific writing.<br /><br />Note also the <span style="font-weight:italic;">extensive</span> amount of extra material at his <a href="http://www.symmys.com/AttilioMeucci/Home/Home.html">homepage</a>. Extra appendices, presentation slides covering the <span style="font-style:italic;">whole</span> book, Matlab code. <br /><br />Last but not least, Meucci donates heavily to charities from his book project and teachings. Just go buy the book. Now.Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-50312161545522772152008-05-12T21:24:00.000+02:002008-06-07T21:37:53.211+02:002008-06-07T21:37:53.211+02:00fstat.no - even cooler<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://www.fstat.no"><img style="float:left; margin:0 10px 10px 0;cursor:pointer; cursor:hand;" src="http://bp3.blogger.com/_wwvWMP6W7uQ/SErhpnypFXI/AAAAAAAAACg/BduQNLYOT_w/s200/fstat_small.png" border="0" alt=""id="BLOGGER_PHOTO_ID_5209224024179676530" /></a><br /><br /><a href="http://www.fstat.no">fstat.no</a> - the place to be for fat-tailed finance. Complete coverage of the Norwegian derivatives market, daily updates of market prices and monte carlo simulations, Normal Inverse Gauss (NIG) distributions, volatility plots, heavyness plots in NIG-triangles. It's all there, and it's getting better by the day. <br /><br />EnjoyJon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-16008043402586072642008-04-17T20:52:00.003+02:002008-06-07T23:49:38.934+02:002008-06-07T23:49:38.934+02:0010 x speed increase in gaussian random numbers in excel<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://en.wikipedia.org/wiki/Normal_Distribution"><img style="float:left; margin:0 10px 10px 0;cursor:pointer; cursor:hand;width: 200px;" src="http://upload.wikimedia.org/wikipedia/commons/thumb/8/8c/Standard_deviation_diagram.svg/325px-Standard_deviation_diagram.svg.png" border="0" alt="" /></a><br />For generating Gaussian random numbers in Excel, one could simply do<br /><br /><blockquote>=NORMSINV(RAND())</blockquote><br /><br />but the performance is horrible. Instead, try this VBA function from <a href="http://www.vbnumericalmethods.com">www.vbnumericalmethods.com</a>. <br />It will easily give you a 10x speed increase:<br /><pre><br />'From Jack at www.vbnumericalmethods.com<br />'Simulate a Gaussian variable N(0,1)<br />Public Function Gauss()<br /> Dim fac As Double, r As Double, V1 As Double, V2 As Double<br /> 10 V1 = 2 * Rnd - 1<br /> V2 = 2 * Rnd - 1<br /> r = V1 ^ 2 + V2 ^ 2<br /> If (r >= 1) Then GoTo 10<br /> fac = Sqr(-2 * Log(r) / r)<br /> Gauss = V2 * fac<br />End Function<br /></pre>Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0tag:blogger.com,1999:blog-1721391884698637233.post-88623825845741336902008-04-03T21:47:00.004+02:002008-06-07T23:34:03.787+02:002008-06-07T23:34:03.787+02:00Young entrepreneur<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://www.universitas.no/?sak=3829"><img style="float:left; margin:0 10px 10px 0;cursor:pointer; cursor:hand;width: 200px;" src="http://static.universitas.no/bilder/elm/unilogo/logo.gif" border="0" alt="" /></a><br /><br />Do you ever google yourself? I did so today, and found <a href="http://www.universitas.no/?sak=3829">this article in Universitas</a> from 2003. Looking back five years feels good, I recognize both who I was and who I am.<br /><br />The article describes a venture we did in 2003-2005. Building your own business and making it cash-positive is a great, although hard, experience. The mindset of business rally can't be learned any other way, and the start-up experience really helped shape who I am as a professional.Jon Egilhttp://www.blogger.com/profile/17785403161734934222noreply@blogger.com0