Colin Gan
2013-Dec-19 19:04 UTC
[R] ATTN: Problem with Portfolio Optimization converging with Expected Shortfall as risk measure on Asymmetric Generalized Hyperbolic Distribution
Hi there, attached are the necessary code and data. For me, the portfolio optimization with expected shortfall as risk measure only converges on gaussian not asymmetric generalized hyperbolic, normal inverse gaussian, variance gamma and student t distributions. I like to know why. Thank you. Sincere Blessings, -- Colin Gan Cheong Kiat -------------- next part -------------- BACreturn,BRKBreturn,GEreturn,KOreturn,LOWreturn,MCDreturn,MMMreturn,RFreturn,TGTreturn,WMTreturn,TBillYield,EURUSDreturn,SPYReturn 0.0615128,0.050409396,0.02373628,0.007538838,0,0.007262196,0.006451073,0.008832865,0.028791648,0.018180683,0.00E+00,-0.012921829,0.014173553 0.022919261,-0.010416761,0.010165272,0.011641054,0.01262584,0.022274095,0.031991354,-0.00631714,0.03518713,-0.001089325,2.30E+00,-0.012126488,0.023426134 -0.031941879,-0.008063259,-0.019602272,-0.003820659,0.005663064,-0.011199063,-0.006660091,-0.037087069,0.004641816,-0.004345471,3.07E+00,0.018745023,-0.00287241 0.059378171,-0.008948606,-0.019668735,0.002952199,0.006139335,0.01033389,0.001081666,0.015883046,0.01244184,0.002709295,3.44E+00,0.001929744,0.005786383 -0.059378171,-0.017660503,0.022154752,-0.015570663,-0.011304878,-0.0122902,-0.016635747,-0.009804,-0.025646392,-0.012581043,3.66E+00,0.000572492,-0.010391503 -0.01975915,-0.004055535,-0.026406201,0.005163329,-0.054239769,0.017746111,-0.004425178,0.015980672,-0.048819217,0.004475081,3.81E+00,0.002149768,-0.012896042 -0.011527505,-0.020949626,0.033274789,-0.016739594,-0.007073416,-0.015952974,-0.001588282,-0.035304945,-0.028638842,0.005937939,3.91E+00,-0.047354414,-0.012395922 -0.046194357,-0.006079046,-0.085809736,-0.006618158,-0.012259348,-0.021627699,-0.032442779,-0.047878379,-0.00871404,0.003978305,3.98E+00,0.038854059,-0.018860111 -0.032304289,-0.008449054,-0.030022079,-0.021319058,0.019332764,-0.004619682,0.004448254,0.008013781,-0.011409613,-0.000905551,4.04E+00,0.026597121,-0.009733767 0.067810978,0.049902562,0.055649574,0.032596127,0.026980054,0.036085138,0.047222464,0.051898438,0.040948824,0.00964082,4.09E+00,0.013235487,0.045991515 -0.056488846,-0.023721462,-0.009395254,0.008934614,-0.00227635,0.001649077,-0.007343097,-0.070122177,-0.032613127,0.001463326,4.12E+00,0.012663111,-0.014032142 0.001340483,-0.029780384,0.003512884,0.007883812,0.021605069,-0.00461058,0.016190316,0.06529125,-0.04660879,-0.002924513,4.14E+00,0.009414492,0.01186924 -0.112777609,-0.030953097,-0.054221554,-0.024350583,-0.057339108,-0.01289498,-0.042257308,-0.077630161,-0.034034699,-0.008360645,4.16E+00,0.022189128,-0.024376382 -0.015458096,-0.028049921,0.012297528,0.010650201,0.031196826,0.032638467,0.024106032,-0.015487035,0.022648371,-0.010084724,4.19E+00,0.006753898,0.006792521 -0.051175293,0.000988352,-0.060028931,-0.025328961,-0.034263641,-0.020564572,-0.021669129,-0.091210379,-0.008912248,-0.030874817,4.20E+00,-0.001245524,-0.016620881 0.157519577,0.0520445,0.056099652,0,0.023009865,-0.029118848,-0.042145621,0.042995363,0.091767402,-0.011229279,4.21E+00,0.014812763,0.034318463 0.060870715,0.024097552,0.016949558,0.027115075,0.039256524,0.030432909,0.011933918,-0.015568555,0.026668247,0.054725973,4.21E+00,0.003559426,0.030366589 -0.048325156,-0.028700039,-0.032514075,-0.030588808,0.010294897,0.021431172,-0.026200915,-0.023410738,-0.012843142,0.013334749,4.20E+00,-0.00023769,-0.023405077 0.046235291,0.011770977,0.029101102,0.049082795,0.073673882,0.042888219,-0.002467715,0.071943469,0.089439972,0.017249802,4.18E+00,-0.005845198,0.023756839 0.111798889,0.016929062,0.127626899,-0.009064191,0.038714512,0.01554527,-0.024187528,0.088093043,0.065194196,0.008266062,4.17E+00,0.007034212,0.037762881 -0.234749097,-0.064633791,-0.102323651,-0.055716163,-0.083760485,-0.044908346,-0.070883625,-0.184876339,-0.133531393,-0.052725959,4.16E+00,-0.009394188,-0.092749359 0.051779242,0.040536239,0.058769907,0.032306333,-0.028627909,0.014573768,0.015357122,0.018821754,-0.039778797,-0.014391294,4.12E+00,-0.01350159,0.012509936 0.041686486,0.044424707,0.026284367,0.006854643,0.035914801,0.019005018,0.035047972,0.060812139,0.093007112,0.036099816,4.10E+00,0.009632639,0.037913262 0.014290818,0.085711097,0.032874375,-0.003543747,0.052551495,0.023149182,0.005988986,0.003193191,0.050290821,0.035555164,4.09E+00,0.002970609,0.007381837 0.240601431,0.012278463,0.084037132,0.02757095,0.126362903,0.006506439,0.031473375,0.165407834,0.080699954,-0.002838492,4.09E+00,0.013478037,0.066995573 0.019366802,0.10635239,0.088632596,0.068911706,0.019429519,0.041464389,0.076419589,0.037159945,0.002853069,0.043644695,4.12E+00,0.004772519,0.052538377 -0.149185995,-0.060355165,-0.118129116,-0.028799687,-0.07002286,-0.046708998,-0.034271267,-0.129607254,0.037850224,-0.006688988,4.10E+00,-0.00730858,-0.077132835 -0.151083193,-0.115904801,-0.105637294,-0.028453412,-0.05709682,-0.019116259,-0.056241048,-0.095101651,-0.108841098,-0.033169257,4.06E+00,-0.001502749,-0.066224216 0.010589113,-0.004788517,-0.003108489,-0.010976557,0.00577581,0.016054572,0.022282779,0.056851769,-0.052490183,0.017411709,4.02E+00,-0.000553076,0.018661799 -0.0884519,-0.076312093,0.005602256,-0.023599252,0.040843936,-0.009307388,-0.036496865,-0.065159189,-0.041730672,-0.01722201,4.00E+00,-0.001184133,-0.013365299 -0.040578359,-0.010566579,-0.051106011,-0.027603067,-0.079071731,-0.002135612,-0.021492709,-0.080440369,-0.067316222,-0.041254456,3.95E+00,0.011904903,-0.05119516 0.005865119,0.019160498,0.03439253,0.054391388,0.072510746,0.040451693,0.080611963,0.08354757,0.059775845,0.042963372,3.93E+00,-0.000399122,0.060473816 -0.094775277,-0.029327615,-0.089208675,-0.01382145,-0.037437527,-0.041162551,-0.048277311,-0.068162651,-0.074386951,-0.039681093,3.99E+00,-0.017170181,-0.044998764 -0.050548664,-0.040697658,-0.021978907,-0.007596347,-0.012282652,0.018044259,0.001849284,-0.025879448,-0.006317835,0.014420255,3.91E+00,0.010527644,-0.013191247 0.018222617,-0.00379096,0.02795881,0.038797019,0.02889778,0.023715084,0.038844415,0.030572089,0.029172787,0.016685593,3.86E+00,-0.001096921,0.03248445 -0.077899412,-0.050076338,-0.083141667,-0.005379973,-0.071224237,-0.033224441,-0.034051209,-0.0777969,-0.062298207,-0.01189384,3.86E+00,-0.007800352,-0.057005542 -0.120283101,-0.015306421,-0.041283504,-0.02864894,-0.05270624,-0.031291004,-0.01675287,-0.085051233,-0.091086946,-0.036281868,3.83E+00,0.00389257,-0.042936402 0.038226507,0.023048395,0.073404542,0.012680537,0.04396832,0.01306525,0.023127765,0.060974403,0.050339126,0.002854596,3.79E+00,-0.000155994,0.033417447 -0.023441844,0.005194817,-0.010822059,0.031060479,-0.047661769,-0.01565794,-0.002959273,0.05951012,-0.020141668,0.002862768,3.76E+00,0.005082305,0.002887493 0.059229237,0.029064133,0.008234735,-0.013525905,0.077587966,0.00103627,0.027593394,0.021878721,0.044077954,0.019175696,3.71E+00,-0.021557907,0.005488545 0.020399821,0.009427679,0.00778214,0.018303528,0.036515213,0.011994929,0.00079968,-0.045933333,-0.058922363,-0.004919387,3.68E+00,0.020539377,0.03400898 -0.030880266,0.051364517,-0.014991234,-0.031427951,-0.013849926,0.009841907,-0.028548583,-0.020906685,0.055801424,-0.002722572,3.64E+00,0.019292186,-0.007278987 0.114476172,-0.002844952,0.094643393,0.098386731,0.152708351,0.089744629,0.063222847,0.153432547,0.163848727,0.105018198,3.56E+00,0.005282971,0.110516942 -0.025962997,-0.053379186,-0.005624312,-0.014524584,-0.020558727,-0.024805699,0.012668948,0.107036962,-0.007011155,-0.034237043,3.52E+00,-0.01085579,-0.0361464 -0.087643988,-0.036232973,-0.052974438,-0.034253969,-0.051554888,-0.046400076,-0.031864025,-0.090825867,-0.030219203,-0.026115155,3.45E+00,-0.016846813,-0.052045539 0.01489296,-0.025262612,-0.008474627,-0.052476987,-0.013710113,0.025511588,0.055973851,-0.0916821,0.015444322,0.009330735,3.43E+00,-0.0025728,0.011517705 -0.056201793,-0.027847214,-0.070749415,-0.014221878,-0.045777018,-0.017382196,-0.031297488,-0.050873822,-0.120804838,-0.026431593,3.40E+00,-0.026205057,-0.055993177 -0.017780084,-0.016852273,0.010373025,-0.00929837,-0.027202088,-0.030546271,0.043052161,0.059296894,-0.009766476,-0.014061286,3.32E+00,-0.018040177,-0.03031002 0.0487082,0.037758958,0.025648879,0.049866938,0.022530395,0.054966836,0.017895214,-0.012144007,-0.033577535,0.012199782,3.34E+00,0.00149098,0.058343802 -0.042541685,0.054528428,-0.013158085,-0.035122525,0.050993522,-0.01219302,-0.013713295,-0.06209822,0.037789175,-0.015684425,3.54E+00,-0.0076549,-0.005989963 0.017891053,-0.008200019,0.032706032,0.034896307,0.01151013,0.054914277,0.046278915,-0.018158735,0.059504174,0.08737761,3.57E+00,-0.008698271,0.040813108 -0.107750604,-0.041790966,-0.079842888,0.010916642,-0.055227112,-0.083156559,-0.082291036,0.026043139,-0.107653517,-0.084076672,3.67E+00,-0.009277928,-0.103637192 0.151954325,0.003289893,-0.007168489,-0.077629906,-0.042385716,-0.021718873,-0.00420204,0.251691431,-0.025562876,-0.001101524,3.60E+00,0.008250932,-0.014910813 0.061253713,-0.019646997,0.123087779,-0.042459208,-0.019364367,0.02409298,0.014105653,0.184512138,-0.014222701,-0.008598841,3.66E+00,-0.007556619,-0.024554805 -0.1185186,-0.019268419,-0.082750163,-0.090680542,-0.046569828,-0.025216996,-0.069839865,-0.249119038,-0.095721887,-0.059674111,3.45E+00,-0.00357495,-0.07239758 -0.072846673,-0.048426594,0.017094433,-0.039090904,-0.021216837,-0.028789588,-0.044322919,-0.087353199,0.035489185,-0.005302142,3.48E+00,0.00723608,-0.025515204 -0.304163096,-0.062263207,-0.05183502,-0.032331693,-0.058181955,-0.038891808,-0.049290289,-0.131061239,-0.060595106,-0.05429753,3.48E+00,-0.000146703,-0.045819981 -0.067792068,-0.057537993,-0.008847555,-0.031889841,-0.038849035,-0.052491177,-0.025079684,-0.108488348,-0.03046393,-0.031117022,3.50E+00,-0.014562657,-0.052276386 -0.053364813,0.009073296,-0.026534036,-0.019403453,-0.041416613,-0.027026459,-0.015361285,0.131842203,-0.02354115,0.014842573,3.31E+00,-0.004975314,-0.01359218 -0.047257115,0.009375409,-0.100835621,-0.010393561,-0.035045198,-0.028359396,-0.026920943,0.125725889,-0.05526749,-0.013669945,3.26E+00,-0.012721727,-0.036948611 0.085653312,0.037950983,-0.040005335,0.023917433,-0.019608471,0.031272491,-0.01355953,0.08483888,-0.037599195,-0.003847767,3.01E+00,-0.015642935,0.000603318 0.145851877,-0.014232704,0.098845835,0.0356115,0.031299447,0.025276251,0.029564658,0.151549898,0.035273366,0.024337858,2.94E+00,-0.003210948,0.040556202 -0.193280571,0.007429957,-0.088993538,-0.02801845,-0.063285559,-0.049296622,-0.046115558,-0.528844129,-0.083432294,-0.037936741,3.08E+00,-0.020077255,-0.081602329 0.065592664,-0.003384098,-0.016886323,0.019627338,0.023997842,0.014824532,-0.004023572,0.010772097,0.02458576,0.009765864,2.84E+00,-0.004089985,0.000496607 0.038557564,-0.007851975,0.043372623,0.024989003,0.008832865,0.013581454,0.01677552,0.061798596,0.019300961,0.020161973,2.67E+00,0.000680504,0.015518416 -0.00693087,0.039891329,-0.014533966,-0.009715552,0.002114613,-0.005681372,-0.006251383,-0.052367986,-0.005062277,0.008864702,2.72E+00,-0.002787127,0.003200272 -0.025205189,-0.012709588,-0.046975368,-0.013717636,-0.007171514,-0.012706243,-0.024339068,-0.128416292,-0.005840315,-0.008355407,2.60E+00,0.010920865,-0.02301444 -0.093044871,-0.053644146,-0.017813622,-0.025745899,-0.070588907,-0.022284605,-0.027620671,-0.238411023,-0.068296511,-0.013660723,2.53E+00,0.023259284,-0.022899587 0.203461257,0.078049773,0.071222467,-0.012628572,0.058887979,0.000312647,0.030453534,0.304660409,-0.04061192,-0.029379898,2.48E+00,-0.018512603,0.033173934 0.117129227,0.020202707,0.058131765,0.013767308,0.059041953,0.022771382,0.023682485,0.295014828,0.035560932,0.030385429,2.99E+00,0.019355892,0.029239873 -0.082866771,-0.002398083,-0.068964369,-0.040746033,-0.072658834,-0.027917593,-0.012700416,-0.076995113,-0.076836632,-0.041063414,3.00E+00,-0.003018499,-0.046005348 0.107053996,0.039077011,0.018526507,0.001824818,-0.004905979,0.008905613,0.012264781,0.054256526,0.043619522,0.008241139,2.76E+00,0.008163885,0.016598919 -0.239654329,0.008755528,-0.08378803,0.004576667,-0.036047155,-0.005321659,-0.018269946,-0.049988487,-0.033022706,-0.012576754,2.81E+00,0.006024755,-0.048754476 0.020359985,0.017488724,-0.051368058,0.006442728,0.005522696,-0.00777488,0.002998075,0.049988487,-0.006440966,-0.012103986,2.37E+00,0.009428641,0.004610501 0.02016567,0.004099416,0.00248889,0.00629865,0.001979807,0.021290723,0.018179696,-0.010733556,0.015563836,0.018372591,2.36E+00,-0.011417278,0.014364489 -0.003696862,-0.014907583,-0.001067426,-0.001300028,-0.016509809,0.00015824,-0.011149033,-0.027206593,0.009206864,0.01445417,2.29E+00,-0.003541331,0.004049572 -0.065748582,-0.010692598,-0.033913406,0.009323208,-0.064538521,0.012260324,-0.024604811,-0.05802384,-0.029272093,-0.014131641,2.27E+00,-0.004936889,-0.030136968 0.074706171,0.021245109,0.042484887,0.027539486,0.060647466,0.034553633,0.028355537,0.098665476,0.050666167,0.020531925,2.18E+00,-0.002318964,0.020445404 0.051897687,-0.008653655,0.006477173,0.004245471,-0.000776398,0.004654264,-0.007917698,0.054668544,0.022047558,0.015931303,2.18E+00,-0.016774362,0.003139467 -0.074294696,-0.017210399,-0.030924803,0.000967399,-0.035827002,-0.027930187,-0.015508599,-0.046520016,-0.021684054,-0.000167266,2.19E+00,0.003249337,-0.030568852 0.010062594,0.003743453,0.001401051,-0.005790404,0.014329054,-0.019097106,-0.004085951,-0.004535155,-0.000545207,0.002344274,2.12E+00,-0.005180106,-0.000859812 0.046738944,0.024680111,0.015186589,-0.002114775,0.066337835,0.014570749,-0.006726483,0.175647048,0.037390993,0.009770967,2.06E+00,-0.007207362,-0.006231035 -0.009269684,-0.007150184,-0.025646935,-0.019964538,-0.010496665,-0.014252784,-0.012491488,0.038487433,-0.016090151,-0.013619378,1.91E+00,-0.002440183,-0.010811751 0.058300725,0.019268419,0.021385565,-0.012534071,0.012121361,0.012963293,0.024858144,0.051765438,0.019301639,0.009901902,1.87E+00,0.000270838,0.012054856 0.021476878,0.000778513,-0.001770225,0.00503216,0.00162734,-0.000805088,0,0.013072082,-0.00189036,0.004903214,1.88E+00,0.011510428,0.009686819 0.00206968,-0.003757211,-0.001767097,-0.000933794,0.017248888,0.000805088,-0.000141393,0.021766276,0.020607484,0.007656354,1.82E+00,-0.011510428,0.002908693 -0.042257428,-0.000904802,-0.027856955,-0.014823321,-0.02938987,-0.021355393,-0.021677655,-0.036018299,-0.015115566,-0.01508631,1.89E+00,-0.002772614,-0.020493958 0.039498805,0.009085072,0.012787456,0.015943979,0.0117339,0.017657343,0.019135998,0.039692523,0.030854742,0.015940646,1.76E+00,-0.000405104,0.01437197 -0.008571971,-0.005201572,0.003484324,-0.015576061,0.011049836,-0.001924003,0.003107347,-0.047913356,0.019180024,0.002396031,1.83E+00,0.005414931,0.001722923 0.042188582,-0.011605546,-0.002092051,0.000368053,0.018692133,-0.000480423,-0.00873245,0.030881314,0.007414121,0.002745369,1.87E+00,0.003876106,0.004635276 -0.042529937,-0.006134989,-0.022047052,0,-0.028520222,-0.013200987,-0.015860099,-0.019116472,-0.006615239,-0.010766574,1.78E+00,-0.001837923,-0.010964163 -0.046675139,-0.019182819,-0.015215847,-0.013346939,0.001631322,-0.005357718,-0.014252674,-0.074107972,-0.026615667,-0.009137119,1.71E+00,-0.001698774,-0.013768782 0.017083201,0.027114292,0.009440394,0.000726744,0.028567457,0.000471587,0.006142106,0.024472857,0.034639782,0.021623384,1.68E+00,-0.012012562,0.004851575 0.0447229,0.022864077,0.007139245,-0.003447341,0.031139219,0.012976919,0.001095891,0.06034313,0.032340109,0.003793762,1.63E+00,0.000268366,0.007516211 -0.075715431,-0.008635406,-0.014564186,0.007636401,-0.041584269,-0.018308592,-0.00873368,-0.080410693,-0.026482098,-0.023393878,1.54E+00,-0.000268366,-0.006048406 -0.069784934,-0.018585976,-0.007036384,-0.010168964,0.01044505,-0.030951736,-0.010542063,-0.09671397,-0.037184918,0.01171391,1.54E+00,-0.00701898,-0.01045922 0.034438875,0.015983841,0.010404524,-0.001083424,0.044222702,0.004254679,0.012173828,0.034661982,0.04982525,0.01202557,1.39E+00,-0.004121524,0.010304613 0.022895778,-0.001558037,0.036767584,0.025590868,0.084221753,0.059768433,0.03180138,0.017792201,0.063089701,0.015677032,1.44E+00,-0.025995898,0.018410692 -0.059429522,-0.006207987,-0.014938614,-0.024688097,-0.014224991,-0.004516137,-0.008811862,-0.068359091,-0.047998975,-0.064582411,1.44E+00,-0.000452357,-0.015003798 -0.003878865,0.007506183,-0.008583744,0.00816183,-0.009838449,-0.003052947,-0.00097432,0.02092752,0.000416667,0.006936444,1.29E+00,-0.000581301,0.004430805 0.029005054,0.005732168,0.037621992,0.018565003,0.04823088,0.028973381,0.030080641,0.055915185,0.064107969,0.032165777,1.27E+00,-0.005088086,0.026605082 -0.021532295,-0.015556445,-0.001418943,0.014200538,0.006872879,0.012965607,-0.005718386,-0.014799424,0.007358713,0.011706107,1.23E+00,-0.000513809,-0.009317209 0.012985235,0.015033979,-0.00283186,0.031348133,-0.000984737,-0.00033456,-0.003415401,0.005260401,-0.012234609,-0.014952586,1.20E+00,-0.002372937,-0.005296664 -0.021350984,0.006287681,-0.023752426,0.003696143,-0.009306949,-0.010978154,-0.01284512,-0.027056801,-0.014486645,0.001024066,1.14E+00,0.001410257,-0.013314733 0.042236264,0.023937313,0.019871667,-0.005054443,-0.010186845,0.012484557,0.011425472,0.022846269,0.004586663,0.019136865,1.03E+00,-0.007413596,0.017660676 0.138242303,-0.004029555,0.025317808,0.008177616,0.053531114,0.03251627,0.012275347,0.09237332,0.038380436,0.025206234,9.62E-01,-0.002607232,0.021127546 -0.052753383,-0.010666768,-0.036174157,-0.017634407,-0.01716347,-0.014771825,-0.018779895,-0.045052664,-0.013780857,-0.014355585,1.06E+00,0.000762389,-0.014772266 -0.035207888,-0.019958646,0,0.003849119,-0.020312811,0.004785516,-0.001408451,-0.11238425,-0.006931274,-0.002460458,8.84E-01,0.003629306,-0.000239053 -0.087446443,-0.020332704,-0.021365389,-0.012265389,-0.038956597,-0.02185972,0.003524854,-0.149641089,-0.029420344,-0.019468719,1.00E+00,-0.004074103,-0.020972071 0.033138767,-0.002798628,0.028706796,0.022148233,0.026770241,-0.007680786,0.007087202,0.034486176,0.005641152,-0.016560344,7.58E-01,-0.011244591,0.005398018 0.124606883,-0.000254033,0.02883275,0.034674103,0.044077449,0.006173374,0.022586014,0.091807549,0.006548812,0.030078748,8.46E-01,0.00383926,0.011280836 0.038184835,0.005859143,-0.011133173,-0.008631992,-0.009072643,-0.011481947,-0.003050776,0.028280773,-0.026155171,-0.01058762,8.71E-01,0.002652019,0.000555489 0.036677569,0.012598187,0,-0.006177164,-0.008493681,0.001158845,-0.007370508,0.032163576,-0.014824499,0.004152255,8.68E-01,-0.002084321,0.006210706 0.156102268,0.019595663,0.011494379,-0.038960907,0.039061522,0.003484613,0.011731614,0.200670695,0.029655507,0.012561226,8.46E-01,-0.0024577,0.009953526 0.202191235,-0.011020841,0.03792098,0.010563817,0.055850727,0.02695581,0.013936997,0.152072109,0.055656472,0.01272102,9.89E-01,-0.00639942,0.024251035 -0.084353059,-0.006502819,-0.019692258,0.016156069,-0.038631413,0.008058338,-0.019457859,-0.035692582,-0.005479466,-0.001243892,9.85E-01,0.008983319,-0.01419744 -0.072724527,-0.013391908,-0.017505918,0.013632534,0.005277057,0.013343929,0.00435604,-0.180449773,-0.018498159,0.00035524,8.29E-01,0.000757528,-0.009085072 -0.031344833,-0.016491557,0.000723327,0.004585792,-0.027141541,-0.019350984,-0.002906134,-0.104423308,-0.003347097,-0.016211809,6.44E-01,0.008052543,-0.011720452 0.013507634,-0.002764168,0.01641474,-0.016843755,-0.040368751,-0.013765193,0.004070956,-0.012578782,-0.052717986,-0.008008399,7.92E-01,-0.000445534,0.004078537 -0.064935088,0.004527171,-0.031495816,-0.011332673,-0.02924185,0.012397204,-0.034226868,-0.131576358,-0.028129785,-0.024663007,7.53E-01,0.001783326,-0.019442742 0.08925379,0.006575643,0.034811347,0.003503312,0.035691884,0.02002138,0.016752237,0.127418348,0.031304727,0.037929042,6.45E-01,0.002297531,0.017601342 -0.039613647,0.011740819,0.007035763,-0.003697581,-0.006450035,0.002968138,0.005599038,-0.067193189,0.005740422,0.005462087,6.22E-01,-0.001277058,-0.010265478 -0.00623055,-0.013515442,0.014975944,0.002139038,-0.000494438,-0.004187757,0.014209315,-0.026795862,0.004487666,0.001768347,4.61E-01,0.004985653,0.001029744 -0.054814152,-0.012835208,-0.022749443,0.007228712,-0.018609744,0.002440726,-0.018805704,-0.090232498,-0.017831082,-0.009336807,3.49E-01,0.001963206,-0.017285157 -0.00251678,-0.00324554,0.0159825,-0.019033396,-0.006769852,0.018853523,0.002726557,0.060462858,0.022123855,0.014660688,2.71E-01,0.000697549,0.003120614 -0.029717416,-0.001121007,0.016242095,0.002696977,-0.005766474,-0.004968073,0.001150252,0.011987234,-0.029251768,-0.001777778,2.84E-01,0.001015486,0.003522371 -0.008906941,-0.005214812,-0.010229299,-0.026272214,-0.015688461,0.000885348,-0.01073202,-0.032849669,0.004816259,-0.00936982,2.32E-01,0.001080085,-0.005473895 -0.070040287,-0.007279938,-0.053571135,-0.01639992,-0.042932972,-0.027949418,-0.032348966,-0.038737238,-0.045947376,-0.022445483,1.35E-01,0.008426478,-0.027535965 -0.000375728,-0.006128221,0.014393915,0.013585396,0.009991,0.017902625,0.001931035,0.019181059,0.014336439,0.013860236,-5.68E-02,0.001989284,0.004714838 0.028192343,0.006620105,0.006910375,0.00206321,0.005951036,0.000175392,-0.007290768,0.044091712,0.008783633,0.011934154,-1.75E-01,0.00302358,-0.001979897 -0.046063258,-0.021657982,0.000730194,-0.007482265,-0.047083343,-0.006817612,-0.00082203,-0.031691199,-0.025122877,0.006731646,-2.76E-01,-0.005717789,-0.00098848 -0.037658324,-0.002164243,-0.019172228,0.005606443,-0.031899256,-0.02051973,-0.014952025,0.003571432,-0.035767224,-0.025100045,-1.99E-01,0.008298269,-0.021354195 -0.008140197,-0.008610433,-0.010691477,0.003755168,0.017979936,0.006506872,-0.00242555,-0.019486888,0.018510027,0.000346741,-3.97E-01,-0.000774843,0.001265493 -0.030205386,-0.001903403,-0.022780028,-0.014936797,-0.012449184,-0.017033285,-0.008977082,-0.110348057,-0.030980212,-0.017532247,-3.26E-01,0.001033258,-0.009784379 -0.036269926,0.002141583,-0.003804258,-0.004068805,-0.025691032,-0.012253615,-0.003063197,-0.077077412,-0.012128252,-0.010508571,-4.48E-01,0.001163693,-0.004856522 0.017970534,-0.01160463,-0.006194101,-0.022628703,-0.003735219,-0.00016682,-0.012159811,0.03477968,-0.008440451,0.002194279,-6.69E-01,0.007988352,0.000587415 0.011482733,0.005193591,0.003436429,-0.030563857,0.036275593,0.010227267,0.007516354,-0.015689526,0.025154991,0.001183532,-6.70E-01,-0.005267791,0.012564836 0.020244301,0.001184133,-0.026496115,0,0.018643506,0.009992462,0.00371304,-0.014716969,0.011558595,0.010031539,-6.56E-01,-0.006336501,0.003800444 -0.026339816,-0.001184133,-0.016622723,-0.015111013,-0.019931336,-0.01721271,-0.009651691,-0.050562567,-0.030154957,-0.021302581,-3.81E-01,-0.00070874,-0.014821673 0.000337667,-0.004015595,0.008941937,0.037766849,0.003437906,0.007725937,0.00237123,0.004175371,0.014010074,0.003519065,-5.27E-01,-0.016480731,-0.004989737 -0.029614541,-0.010319041,0.001331558,0.001074691,-0.012406576,0.040604236,0.012742408,-0.074556168,-0.001333715,0.020350032,-6.60E-02,0.011918181,0.002418381 -0.04769667,-0.015970713,-0.034056992,-0.024605189,-0.045707998,-0.019131018,-0.034658234,-0.068172686,-0.039389133,-0.024203602,2.15E-01,0.012581217,-0.03241342 0,0.001148765,0.019834824,0.008604847,0.023840011,0.001206585,0.014818938,0.006672758,0.007902274,0.036095168,1.32E-01,-0.004146695,0.019799786 -0.040434298,-0.027212564,-0.000328353,0.005305052,0.002916061,0.007270241,0.001703913,-0.015101465,-0.014652277,-0.001559117,1.32E-01,-0.008178029,-0.000507044 -0.008072999,0.018059182,0.001642846,-0.010582109,-0.003747662,-0.017053208,-0.000262329,-0.035917903,0.040443278,0.009915711,1.10E-01,0.00462815,-0.005848603 -0.012723947,-0.024303626,-0.010142407,-0.004551032,0.00249688,-0.013065427,-0.017034415,-0.030193122,-0.010171495,-0.009396276,2.63E-02,0.000837926,-0.010385053 -0.017199083,0.021801066,0.002607563,-0.010596816,-0.016119382,-0.002693604,-0.00257533,-0.041225953,-0.007282265,-0.003457818,-1.45E-01,-0.002768925,0.00249688 0.021324015,0.014418375,0.002941659,0.007458191,0.009474839,0.013370767,0.011121284,0.02674102,0.02696496,0.014779022,-3.86E-02,-0.006729285,0.005012542 -0.008818399,0.009961857,0.004922077,-0.019823109,0.027272511,0.009931588,0.009538192,-0.051735674,0.003254526,0.012159811,-1.67E-01,-0.006620854,0.004604987 0.007048487,0.039174148,-0.000986355,-0.000682477,0.008972502,0.006560797,0.004737474,0.020159834,0.002111528,0.011591748,-5.99E-02,-6.34E-05,0.007383313 -0.023304288,0.005949869,-0.018880769,0.00615913,-0.006843483,-0.013762474,-0.010888924,-0.025398191,-0.018467921,-0.005366739,5.75E-02,-0.000825475,-0.013494702 0.002883508,0.011144565,0.000645161,0.028195394,0.003415887,-0.004262218,0.006808091,0.012618464,0.002644505,0.014737866,-1.67E-01,0.000126952,0.000143369 -0.021708936,-0.010047873,-0.023282825,-0.007736983,-0.016122538,-0.015195295,-0.018353916,-0.042978398,-0.026503978,-0.012952151,-1.09E-02,0.007263029,-0.017059195 -0.019863571,-0.00450095,-0.021211026,-0.007503743,-0.020413123,-0.012819621,-0.019415605,0.001014199,-0.011537533,-0.008010724,1.31E-02,0.003972327,-0.00814211 -0.001937181,0.010370373,0.00836825,0.007854129,-0.026049499,-0.001322533,0.004945797,-0.023070231,0.000728597,-0.011283618,-5.20E-02,0.005084484,0.002730043 -0.014819152,0.002947681,-0.007441895,0.001051894,-0.009994087,-0.005437035,0.010222425,-0.055449247,-0.004000733,-0.001401542,-1.59E-01,0.001549787,0.000911673 -0.002448648,-0.013318054,-0.004315666,0.008278339,0.0095924,0.006925007,0.001671059,0.007060514,0.022947313,-0.005760686,-1.57E-01,0.00226442,0.012425151 0.005176418,-0.006894494,0.005552143,0.000707714,0.020694591,-0.012005736,0.007100927,-0.012206724,0.011204599,0.014023017,-2.46E-01,-0.002199794,0.002062223 -0.022418224,0.012859573,-0.002162831,0.00177148,-0.012225091,0.004916431,-0.006328726,-0.014362072,-0.004309948,-0.023895791,-1.98E-01,0.002782903,0.000142379 0.021326223,-0.020305266,0.004020417,0.004442477,0.02004567,0.027143453,0.017534055,0.026568797,0.023459267,0.014583592,-5.77E-02,-0.001813002,0.011168501 -0.018383865,-0.023881026,-0.009867487,-0.002490218,0.002896754,-0.008906873,-0.008220833,-0.013605652,-0.001912412,0.000349834,1.16E-01,0.007206653,-0.001870101 -0.017788867,0.004097167,0.000613874,-0.013060546,-0.032616018,0.006546394,0.003514942,-0.020752416,-0.020798632,0.005789998,9.70E-02,-0.005393651,-0.00258361 -0.032110475,-0.019169916,-0.013115942,-0.021509936,-0.0190709,-0.018520593,-0.013986242,-0.042445505,0,0.008482115,7.32E-03,-0.006331179,-0.017971148 0.006904515,0.002650229,-0.005439723,-0.001885984,-0.001572946,-0.00626238,0.002704271,-0.001311476,0.00468912,-0.010942573,-5.93E-02,0.004388803,0.008625619 -0.020823505,-0.026079027,-0.004810593,-0.006658155,-0.026368337,-0.001313413,-0.012812475,-0.025022876,-0.013445581,-0.009260134,-1.26E-01,0.000129374,-0.004816897 0.010103647,-0.002694209,0.006620548,-0.005430185,-0.005343524,0.000492328,0.009079927,0.018491174,-0.006102654,-0.009864232,-1.61E-01,-0.005290676,0.002759794 -0.008488115,0.016742473,-0.003663008,0.002551239,-0.013015368,-0.010352405,0.000520291,-0.019426859,-0.018092445,-0.010807207,9.06E-02,-0.001992865,-0.005913337 -0.008416673,0.02588543,-0.010607767,-0.004418019,-0.004691174,0.010184578,-0.003506725,0.003136906,0.002961868,0.021732487,-1.47E-01,-0.003653967,-0.00394676 -0.003138078,0.009645999,-0.004812039,-0.005916677,-0.014711841,-0.001341607,-0.008906158,0.002696631,-0.00074129,-0.005217403,-1.91E-01,0.002049706,0.000214877 0.011290654,0.001182732,0.015724545,-0.015387492,0.039189968,0.001509434,-0.016693602,0.016334301,0.020210269,0.003475242,-2.98E-01,-0.011031516,0.009211349 0.02703229,-0.002600167,0.013810255,-0.001161151,0.02386362,0.015050593,-0.018778719,0.067654074,0.0098785,0.009268231,-1.69E-01,-0.010785686,0.004347203 -0.00379096,-0.001179663,0.0009275,0.007154177,0.004512828,0.005809996,-0.006798123,-0.018909654,-0.00476146,0.006521571,1.18E-02,0.000564777,-0.00159617 -0.01635206,-0.007985012,-0.00401297,-0.007651401,-0.020753542,-0.005469164,-0.014188105,-0.008132071,-0.035831413,0.003010717,-3.14E-02,0.002073581,-0.004412463 -0.024945497,-0.003036321,-0.007060659,0.00398473,-0.000805153,0.00632643,-0.006778046,-0.028647181,-0.00091617,0.001242346,-1.00E-01,0.00745,0.00050536 0.028674857,0.001400234,0.020708513,-0.009437937,0.012960894,0.016778046,0.027516258,-0.00415801,0.022035888,-0.006372831,-1.30E-01,-0.005813235,0.010380088 0.013163387,-0.008835212,-0.006536988,-0.007715708,0.002040401,0.00261986,-0.000495847,0.066230439,-0.004110619,-0.007033614,5.60E-02,-0.003522016,0.001460387 0.038593544,0.006968669,0.007786988,0.003440089,0.011503824,0.021567037,0.02319417,0.008906541,0.020533829,0.014117051,3.46E-01,0.006298043,0.026733422 0,0.000932836,0.007217982,-0.000984091,0.007049583,0.007533668,0.007383866,0.080705618,0.011100592,0.020104788,5.62E-01,-0.002586997,0.002329339 -0.037511293,-0.015741066,-0.009404458,0.011707612,0.001665973,0.002523888,-0.002679769,-0.06315691,-0.000192437,0.006366552,6.24E-01,0.002271151,-0.003379523 -0.020099833,-0.005725425,-0.136841042,-0.017102867,-0.022259781,-0.008089932,-0.023675706,-0.009564634,-0.022639605,0.002558013,6.04E-01,-0.002649175,-0.019599737 -0.006347548,-0.002508839,0.008471153,0.009008332,0.050156772,0.001971503,0.009503635,-0.010961744,0.015545337,0.009558896,3.53E-01,0.009430151,0.00139783 -0.011794135,-0.005678607,-0.013898575,0.004452147,-0.029146731,-0.002329541,-0.003637052,-0.016708005,-0.017424683,-0.005709565,3.63E-01,0.00203692,-0.007262089 -0.028764181,0.002267575,-0.008354718,0.002316348,-0.007878954,0.003945486,-0.00387331,-0.030239885,0.009242733,0.004048589,1.36E-01,6.37E-05,-0.001022719 0.002281081,0.008206109,-0.008019289,0.00198972,-0.002063133,-0.001257071,-0.004106782,0.011412392,0.001896814,-0.003129893,1.27E-02,-0.001846252,0.000511229 -0.024067344,0.006199079,-0.00742709,-0.005132037,-0.002470154,-0.009289099,0.004356221,-0.042141067,-0.010951771,-0.009695493,2.03E-01,0.012544967,-0.00109517 0.001735466,-0.005168564,-0.004745593,0.00165262,-0.010634047,-0.007440247,-0.009188046,0.004595596,-0.003748833,0.00859316,5.21E-02,-0.006804057,0.00248411 -0.01380012,-0.017150741,-0.010726068,-0.016079099,0.003668233,-0.004402578,-0.005054563,0.012048339,0.003373317,0.007002056,6.24E-02,-0.001789595,0.000658593 0.07493666,-0.007404948,0.037650258,0.009320646,0.065382759,0.020234982,0.027047202,0.082571154,0.049841965,0.026228555,1.82E-01,-0.009659455,0.034555527 -0.004211641,0.004705891,0.010866719,-0.001149331,0.011839647,0.005213495,0.015277191,0.005076153,0.01972765,0.010687125,3.56E-01,0.001012402,0.00349173 -0.014861439,0.018590424,-0.005991303,-0.006053188,-0.015320939,0.000540882,0.001155253,-0.045756365,-0.025629693,-0.004785156,2.82E-01,0.000633272,-0.009610729 -0.030583423,0.020342819,-0.008112538,0.0026131,-0.008219817,-0.002701487,-0.010222425,-0.029223637,-0.028045548,-0.010069431,2.39E-01,0.004826019,-0.003158135 -0.027968659,-0.004194833,-0.003763445,-0.004242826,-0.02004759,-0.010733556,-0.004059889,-0.042520268,-0.019080583,-0.002831526,2.51E-01,0.009015695,-0.012311291 -0.035486821,-0.009489714,-0.003481991,0.00375235,-0.025433933,0.007680666,0.00431418,-0.001808319,-0.003735529,-0.010873748,1.31E-02,0.009421856,0.000964499 0.082115247,0.019156562,0.052009531,0.017851714,0.054856742,0.014063888,-0.018927208,0.062698594,0.043132845,0.047119571,8.10E-01,-0.017291127,0.013414835 -0.009549692,-0.009333868,-0.015335586,0.009215111,-0.017132971,-0.021934342,-0.018452088,0.019302752,-0.016810442,-0.003930823,7.97E-01,-0.005021782,-0.025081822 0.079357952,0.014151455,0.051380678,0.029378714,0.085833544,0.010446352,0.030605675,0.025817683,0.050065737,0.020410894,5.84E-01,6.34E-05,0.040703515 0.007536672,-0.017050512,0.014967264,0.002777297,-0.028333228,-0.009169329,0.012561047,-0.037740328,-0.027053414,0.002605995,7.99E-01,0.013341617,-0.01015867 -0.039824018,-0.003698571,-0.012634221,-0.024043748,-0.025246783,-0.002006019,-0.01979504,-0.06036315,-0.014543538,-0.015535066,6.15E-01,-0.001027683,-0.015616945 0.002966161,-0.005292844,0.008503203,-0.00440231,-0.006901339,0.021918181,0.008381001,0.006529874,-0.010428034,0.006145326,1.08E-01,0.006440416,0.00220524 -0.018462063,-0.009820794,0.016627461,-0.004719375,-0.018174058,-0.010189991,0.006011401,-0.025411393,-0.009545217,0.000795703,-2.41E-02,0.006352108,-0.009395432 0.066024248,0.008444647,0.052239219,0.01371392,0.027386467,0.008329525,0.032199875,0.093669191,0.006808703,0.028256168,-3.03E-01,0.002539147,0.035306814 -0.03969987,-0.022437202,-0.016581015,-0.00323377,-0.035904788,0.028850874,-0.013553731,-0.019841921,-0.014918421,-0.021266624,2.97E-01,-0.005007001,-0.013270925 0.006006024,-0.013357278,-0.019358417,-0.004239091,-0.010215501,-0.014623753,-0.018646144,0.019341044,-0.005178874,-0.001601923,-1.38E-02,0.006376909,-0.010354044 -0.027813322,-0.01905654,-0.024351121,-0.009766031,-0.048306071,-0.016825972,-0.012493788,-0.047430058,-0.010088598,0.008640665,-1.48E-01,0.008061653,-0.020915104 -0.032231272,0.007404215,0.004464958,0.005207908,-0.013383721,0.005763705,0.006609913,0.006227565,-0.004912158,-0.006437359,-5.52E-01,-0.000657851,0.0062964 -0.010261763,0.014531299,0.003586376,0.005912679,-0.010743905,0.008802379,-0.001019758,0.010628119,0.002830992,-0.000400962,-5.69E-01,0.00019731,-0.00382754 -0.014191825,-0.036046806,0.007814888,0.00953361,0.014907108,-0.017712778,0.001147301,-0.029195661,0.005686141,0.006031381,-6.46E-01,0.002370762,-0.002394136 -0.041405619,0.003106757,-0.021197987,-0.016961058,-0.029190151,-0.009381104,-0.017699577,-0.034310273,-0.02273893,-0.02213671,-7.90E-01,0.003037509,-0.022535957 -0.034641661,-0.005136998,-0.005009588,-0.009873731,-0.032293019,-0.012915141,-0.008485202,-0.056702521,-0.018592344,-0.014295749,-1.26E+00,0.005105266,-0.009815048 0.003738322,0.012459882,0.002354327,0.008865158,0.019810474,-0.020040034,-0.010014301,0.001292825,-0.001093693,0.000583487,-1.59E+00,0.011364518,-0.001012365 -0.005136598,-0.006464147,-0.007923739,0.01268087,0.019802627,0.007111141,0.003573853,0.003455727,0.030333378,0.021036852,-1.76E+00,0.003839552,0.007472196 0.007949539,0.001719691,0.019481136,0.008715769,0.037850224,0.006443552,0.016556051,0,0.022408837,0.007178495,-1.97E+00,-0.002090006,0.012529931 0.009197093,-0.011338232,-0.004164194,0.00464517,0.015378347,0.008294314,0.008951708,0.027639399,0.013729308,0.003608665,-1.88E+00,0.007571697,0.006138846 -0.017845038,-0.015618722,-0.01969204,-0.004130104,-0.003437906,-0.017410473,-0.012334958,-0.032818522,-0.038951087,0.001809228,-1.80E+00,0.005443672,-0.008348466 0.007006101,-0.003136437,0.002331003,0.002235022,0.007319731,0.009475355,0.001126479,0.031484892,0.022345183,0.000805153,-2.17E+00,-0.005850755,0.002947246 -0.000702823,-0.001043297,-0.002621997,-0.011295687,-0.008177364,0.006668494,-0.001251565,-0.018486442,-0.016146234,0.004439967,-2.03E+00,0.007216316,0.002807953 0.010831281,0.005437069,-0.000581734,0.003580261,-0.009385734,-0.005050516,-0.004243641,0.004808752,-0.001318143,-0.010663012,-1.83E+00,0.003276453,-0.000221967 -0.025477596,0.00020971,-0.017010645,-0.014244772,-0.025986374,-0.018892025,-0.013484462,-0.033185549,-0.025454231,-0.013713819,-1.95E+00,0.002738602,-0.008838534 0.011839953,0.023126097,0.017592379,-0.008884477,-0.00082713,0.004778347,0.028416668,-0.005074007,0.011068182,0.008126098,-1.90E+00,0.003296253,0.010171078 0.016007875,0.003009461,0.010529492,0.007034027,-0.003301695,0.006942618,-0.002273015,0.029967996,-0.0051804,0.013019713,-2.02E+00,-0.000275103,0.009228316 -0.000474496,0.011475713,0.005011064,0.004210533,0.027570669,0.006092117,0.009123226,-0.008654316,0.03168322,0.017080529,-2.04E+00,0.001996215,0.00509708 -0.028296054,0.011829273,-0.011458918,0.013594092,-0.015549799,0.02143588,-0.007103023,-0.041769413,-0.030021016,-0.02190757,-2.07E+00,-0.003851184,-0.00644197 0.024271903,-0.000440626,0.000876808,0.012741219,0.013011726,0.013494983,-0.000758054,0.021299615,0.059008553,0.020472996,-2.19E+00,-0.003562868,0.006592278 -0.00094451,-0.001100716,-0.000292355,0.005386165,-0.017588393,-0.002045181,0.001137297,-0.022951143,-0.03373814,-0.014839179,-2.11E+00,-0.004571982,-0.008084481 -0.038430655,0.009728133,-0.033346004,-0.021202208,-0.02824006,-0.000742666,-0.018539927,-0.034863685,-0.024345772,-0.01023797,-2.16E+00,-0.009553889,-0.027139007 -0.022457715,-0.004875896,-0.022089532,-0.010688031,-0.03060044,-0.006290493,-0.007912021,-0.02478978,-0.054018532,-0.021926804,-2.54E+00,-0.004037962,-0.012689428 0.019736007,-0.005951742,0.022372298,0.004397098,-0.033862412,0.01187398,0.019396356,0.022004817,0.028987537,0.008634277,-2.48E+00,0.001276925,0.015959898 0.044935767,-0.014401296,0.011948933,0.021067834,0.042115554,0.049154172,0.018116667,0.033531722,0.0564157,0.031634269,-2.43E+00,0.004043132,0.018070126 0.006417134,-0.006262843,0.005165004,-0.01341148,-0.017602645,0.005110074,0.002560165,-0.02395081,-0.007986352,0.003055926,-3.09E+00,0.002501438,-0.007385013 0.017801769,0.027279378,0.001151411,-0.014581472,0.015237646,-0.006285623,0.007461761,0.046806573,-0.012609558,0.006140011,-3.21E+00,0.001219182,0.004941924 0.042644042,-0.024044874,0.020955365,0.016975454,0.027541154,-0.057636945,0.025238345,0.042943481,0.035599076,0.012810093,-3.13E+00,0.003394436,0.016401143 -0.010582109,0.002378637,-0.018937201,-0.024855451,-0.027934777,0.001850139,-0.002909285,0.014140757,-0.026958903,-0.016907619,-3.34E+00,0.002860258,-0.014550867 -0.016652553,0.004121928,0.001733103,0.017010937,-0.02449102,0.018503472,-0.004873882,0.028437935,-0.025699738,-0.025236022,-3.56E+00,0.006088195,0.008406224 0.081625588,0.006761941,0.01573459,0.001359851,0.059329378,0.01062024,0.01242909,0.081531911,0.010349381,0.019324273,-2.87E+00,0.005504349,0.023736847 0.038718031,0.036554103,-0.007606824,-0.032630687,0.100662702,0.000953744,0.003331782,0.04782193,0.071809924,0.033481119,-3.24E+00,0.000828272,-0.010198734 -0.025797259,0.01302429,0.032585821,-0.021176879,0.018190588,0.008047562,-0.000667245,-0.019608471,0.003198723,0.002103935,-4.71E+00,-0.001157841,-0.010320631 -0.047098652,0.028696708,-0.039845909,-0.009784343,0.007832337,-0.008238384,-0.030739318,-0.078114826,-0.004794256,-0.000210593,-5.15E+00,-0.006851433,-0.026257934 0.0211579,-0.009657359,0.000868433,-0.015207009,0.040109174,-0.025432283,0.001553599,0.039523216,0.006196921,0.010584349,-5.35E+00,-0.00633468,-0.008649881 -0.034763552,-0.014313388,-0.026576789,-0.010788941,-0.003364579,-0.022075952,-0.016958269,-0.003435586,-0.012355677,-0.014367456,-5.42E+00,-0.000604412,-0.022258926 0.018528579,-0.01137355,0.008211854,0.008276771,0.025760841,0.011713165,0.011788963,-0.002935423,0.010752792,-0.001048328,-5.55E+00,0.006939092,0.008030459 -0.020566278,-0.019009398,-0.021100797,-0.027682894,-0.02287761,-0.068477275,-0.033338533,0.020232617,-0.031529705,-0.014149216,-5.67E+00,0.008827381,-0.008101238 0.014351861,0.033275632,0.00334635,0.004433241,0.005791522,0.012106685,0.000249377,0.024219145,0.032130486,0.031482138,-6.16E+00,-0.00122691,0.006532722 0.008554816,-0.02291051,0.011236073,0.026390852,0.006799443,0.00680928,-0.000249377,0.004094172,0.019622355,0.020028666,-6.19E+00,-0.001701549,0.010455552 -0.038058482,-0.00519599,-0.02179466,0.004098366,-0.028346842,-0.016506326,-0.011158063,-0.066939483,-0.000816827,-0.012752794,-6.40E+00,-0.001223325,-0.01628026 0.001253919,-0.024923408,0.003877047,0.023331196,0.001422138,0.017032042,-0.007981868,0.002400962,0.018749905,0.018205964,-6.50E+00,-0.000271647,-0.000849558 -0.011229064,-0.021739987,-0.020868412,0.001942062,-0.031748698,-0.015307309,-0.011431476,-0.112641636,-0.030519522,-0.014332493,-6.42E+00,-0.001763908,-0.0248117 -0.00643089,0.00972455,0.001087548,0.01042185,-0.01911756,-0.002930279,-0.000120912,0,0.001211143,-0.011149344,-6.70E+00,0.004415618,-0.000483108 -------------- next part -------------- # check code again to ensure it is right # bond return calculation is wrong except t-bill # security market line should only have stocks # Loading generalized hyperbolic distribution facilities library(ghyp) library(fPortfolio) library(stats) library(timeSeries) library(fBasics) library(fAssets) library(fCopulae) library(timeDate) library(robustbase) library(stabledist) library(mnormt) library(sn) library(quadprog) library(Rglpk) library(slam) library(zoo) library(PerformanceAnalytics) library(utils) library(graphics) library(ggplot2) library(plyr) mydata <- read.csv(file.choose(), header = TRUE, sep = ",") print(mydata, zero.print = ".")# nicer to read View(mydata) ghypmv.fit <- fit.ghypmv(data = mydata, opt.pars = c(lambda = FALSE), lambda = 2, control = list(rel.tol = 1e-5, abstol = 1e-5), reltol = 0.01) summary(ghypmv.fit) nigmv.fit <- fit.NIGmv(data = mydata, opt.pars = c(alpha.bar = T, mu = T, sigma = T, gamma = !0), symmetric = F) summary(nigmv.fit) vgmv.fit <- fit.VGmv(data = mydata, lambda = 1, opt.pars = c(lambda = T, mu = T, sigma = T, gamma = !0), symmetric = F) summary(vgmv.fit) tmv.fit <- fit.tmv(data = mydata, nu = 3.5, opt.pars = c(lambda = T, mu = T, sigma = T, gamma = !0), symmetric = F) summary(tmv.fit) # only this converges during portfolio optimization gaussmv.fit <- fit.gaussmv(data = mydata, na.rm = T, save.data = T) summary(gaussmv.fit) # Using expected shortfall as risk measure with target return from 0.0001 to 0.0020 g1_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0001, distr = "return") print(g1_es) g2_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0002, distr = "return") print(g2_es) g3_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0003, distr = "return") print(g3_es) g4_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0004, distr = "return") print(g4_es) g5_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0005, distr = "return") print(g5_es) g6_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0006, distr = "return") print(g6_es) g7_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0007, distr = "return") print(g7_es) g8_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0008, distr = "return") print(g8_es) g9_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0009, distr = "return") print(g9_es) g10_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0010, distr = "return") print(g10_es) g11_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0011, distr = "return") print(g11_es) g12_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0012, distr = "return") print(g12_es) g13_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0013, distr = "return") print(g13_es) g14_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0014, distr = "return") print(g14_es) g15_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0015, distr = "return") print(g15_es) g16_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0016, distr = "return") print(g16_es) g17_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0017, distr = "return") print(g17_es) g18_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0018, distr = "return") print(g18_es) g19_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0019, distr = "return") print(g19_es) g20_es <- portfolio.optimize(ghypmv.fit, risk.measure = "expected.shortfall", type = "target.return", level = 0.95, target.return = 0.0020, distr = "return") print(g20_es) plot(c(g1_es$risk, g2_es$risk, g3_es$risk, g4_es$risk, g5_es$risk, g6_es$risk, g7_es$risk, g8_es$risk, g9_es$risk, g10_es$risk, g11_es$risk, g12_es$risk, g13_es$risk, g14_es$risk, g15_es$risk, g16_es$risk, g17_es$risk, g18_es$risk, g19_es$risk, g20_es$risk), c(mean(g1_es$portfolio.dist), mean(g2_es$portfolio.dist), mean(g3_es$portfolio.dist), mean(g4_es$portfolio.dist), mean(g5_es$portfolio.dist), mean(g6_es$portfolio.dist), mean(g7_es$portfolio.dist), mean(g8_es$portfolio.dist), mean(g9_es$portfolio.dist), mean(g10_es$portfolio.dist), mean(g11_es$portfolio.dist), mean(g12_es$portfolio.dist), mean(g13_es$portfolio.dist), mean(g14_es$portfolio.dist), mean(g15_es$portfolio.dist), mean(g16_es$portfolio.dist), mean(g17_es$portfolio.dist), mean(g18_es$portfolio.dist), mean(g19_es$portfolio.dist), mean(g20_es$portfolio.dist)), type = "l", xlab = "Expected shortfall", ylab = "Expected portfolio return", main = "Efficient frontier for Gaussian Distribution")