[R] Fractional brownian motion

ED edzabc at gmail.com
Wed Jan 10 14:53:16 CET 2007


Dear All;

I have used fbmSim to simulate a fbm sequence, however, when I tried to
estimate the Hurst effect, none of the nine procedures gave me an answer
close enough to the real value, which is 0.5 (n=1000). So, would you please
advice,

1. which is the best method to estimate the H among the 9 mehods, R/S,
higuchi or Whittle?

2. how to choose the levels (default=50), minnpts, cutoff values or if there
is any other issues to consider? Would you please send your code if you can
get the right estimate for the attached dataset?

3. I also simulated multiple sequences, but some of the estimated results
have H>1, how would you explain this, and how to correct it.

4. if I have a sample size of 30, what are your suggestions for estimating
the hurst effect.

I really appreciate your help, and thanks in advance.

Sincerely;

Qing

list(c(0, 0.0602883614054335, 0.0792551804556156, 0.109563140496959,
0.125726475390261, 0.101877497256596, 0.0571926667502657,
0.00262554012667387,
-0.0651825041744684, -0.0696503406796901, -0.0688337926741136,
-0.0361989921035026, -0.0496524166395568, -0.0718218507093852,
-0.109583486993751, -0.0451891749832022, -0.0260755094163026,
-0.0201228029445018, -0.08818052945358, -0.0285680272562049,
-0.032752629105123, -0.00199553309054604, 0.00508312076775159,
-0.0382257524810562, -0.0398084284989033, -0.0502088451482971,
-0.0502183047870017, -0.0174725900800023, -0.031457396692228,
-0.0497298422601627, -0.0499422794745919, -0.0361646489128442,
0.0266653302426084, -0.00219111075589404, 0.0200305010061950,
-0.054431205284715, -0.00789973517059908, -0.0085425214867059,
-0.0541809440614641, -0.110006586793787, -0.127781398474619,
-0.149159749479771, -0.0888724748560444, -0.0768016585462973,
-0.0427333724364433, -0.0562443902614682, -0.062866380084345,
-0.0883889349587378, -0.08412713827269, -0.124694662643563, -
0.0692092347299578,
-0.062046981143375, -0.0853376756560203, -0.0994827776216755,
-0.0599248702862878, -0.086369926350214, -0.0890587510305954,
-0.0599052704934926, -0.0393644348210234, -0.0364342257334110,
-0.0352152731793845, 0.00857219866032972, 0.0180050537574112,
-0.0356947245542089, -0.0615077884110117, -0.0650043906408892,
-0.0357502728853327, 0.00485292461594026, -0.0458490824792766,
-0.0735087346178546, -0.0820482814570969, -0.124917946967105,
-0.102342757999288, -0.154497291254316, -0.181381845540350, -
0.166583783090691,
-0.130709738638527, -0.124693013403996, -0.120174007595549, -
0.204900741024172,
-0.187508161990262, -0.200032164010717, -0.224311714556569, -
0.207979212824500,
-0.232450923263701, -0.198989498760939, -0.215501785889735, -
0.226880158415004,
-0.203545247216445, -0.182379290331033, -0.197038481353429, -
0.183796217531869,
-0.215299419717715, -0.226095790998872, -0.201080854888686, -
0.261516787028393,
-0.228388047078831, -0.233878235687051, -0.250375640774054, -
0.245746001611628,
-0.232938343176222, -0.204702689747725, -0.259350393091606, -
0.266266298423241,
-0.261596784773769, -0.259875864886204, -0.288719457448324, -
0.318140871065267,
-0.304593775001704, -0.263467947228369, -0.323808258544402, -
0.368149351192211,
-0.347335883864181, -0.400024671714471, -0.371544405921674, -
0.34497157633729,
-0.418657437186672, -0.446912741927435, -0.443280747761798, -
0.384739689692638,
-0.408072099289829, -0.44994167683647, -0.495273953054503, -
0.495799924479583,
-0.43697383254148, -0.440650597876129, -0.458929575025839, -
0.464585181026337,
-0.455843249117235, -0.505571057794026, -0.528316203388037, -
0.58236723395483,
-0.519865988395428, -0.555187599676786, -0.532496699903787, -
0.522718022822901,
-0.559598211031558, -0.547080566663308, -0.554634133071482, -
0.541090495738355,
-0.558915651522986, -0.541917600099341, -0.578140869146858, -
0.53827826166188,
-0.52995892714492, -0.502706986562639, -0.497002923985933, -
0.488096453915811,
-0.51306171001298, -0.536938942207764, -0.515525690037527, -
0.527192086185579,
-0.511195825348352, -0.531104487863626, -0.521038783472886, -
0.519956505596612,
-0.510535242028159, -0.531234776730198, -0.526733559471536, -
0.583555032968273,
-0.551226199377139, -0.489093347568197, -0.505836247919194, -
0.485420121220644,
-0.511258369655571, -0.505044073824655, -0.496646898991666, -
0.49767793016268,
-0.496864531091594, -0.545965735276333, -0.588151447068976, -
0.594481697088302,
-0.633356216248837, -0.649947987535778, -0.698052359479746, -
0.651147537552288,
-0.678918599568793, -0.707817506128165, -0.69962378236799, -
0.664061266573256,
-0.686809968963275, -0.70940931290113, -0.710042179498104, -
0.700294406700751,
-0.69089796594532, -0.687411073827526, -0.660740247977237, -0.61857834066703,

-0.611492191585365, -0.617870447722936, -0.590305790642949, -
0.616713901938279,
-0.619126236382548, -0.618591413756597, -0.604886044395737, -
0.598434945889939,
-0.481524771809731, -0.424388823122359, -0.490492676413456, -
0.506402791875861,
-0.480172973435785, -0.517922095978117, -0.502318428525046, -
0.496731810144271,
-0.531408048623419, -0.546479167267974, -0.583073243018783, -
0.532934946326789,
-0.551949267354811, -0.526047642710959, -0.510085874249115, -
0.516695891848832,
-0.577524134973206, -0.552023065025717, -0.486593593181791, -
0.505244962413324,
-0.50233923782961, -0.520060671266441, -0.528767509770539, -
0.540028769162771,
-0.547477720179665, -0.56531037947022, -0.50193374661204, -0.558262849369318,

-0.494995083625537, -0.497213426057201, -0.500258491337576, -
0.461559035461970,
-0.432242682833388, -0.41117570817118, -0.415407961563806, -
0.372008266805533,
-0.372088740001050, -0.391799048717280, -0.377849107962199, -
0.344184542514291,
-0.348413544886666, -0.380049118743689, -0.39517058021528, -
0.435019093917616,
-0.396825266989612, -0.406175696404068, -0.363939520468838, -
0.276849706077195,
-0.299937554612409, -0.268322605052497, -0.291469417623139, -
0.287282726264279,
-0.264242114241804, -0.276046060096993, -0.217798594356271, -
0.175986521497293,
-0.165844402989144, -0.20246139574587, -0.222148041601722, -
0.204831806500164,
-0.153457268552386, -0.176174002538195, -0.156751478372055, -
0.057878180960031,
-0.0832537960834704, -0.032455515861184, -0.03837970395084,
0.0108416606722306,
0.00435135652516059, -0.0183482961380528, 0.000731922390127772,
0.00999593602758636, 0.0181517073212863, -0.0110713504391765,
-0.00717138666398991, -0.0244655145429809, -0.0435081997698830,
-0.0723796931986123, -0.0812497739510341, -0.167671363055567,
-0.153618438366021, -0.177425107515775, -0.233451742001837, -
0.294328321350601,
-0.262070533304864, -0.304696117238305, -0.308367573404234, -
0.349836635531722,
-0.378702938622887, -0.371478566955775, -0.406474391212437, -
0.427722230759569,
-0.433991757828819, -0.450189851480229, -0.472370053196594, -
0.489099244001516,
-0.489467600993215, -0.485596891281682, -0.49779042627552, -0.46165300963434,

-0.448086803874086, -0.443670308509739, -0.45288256928109, -
0.467547617361111,
-0.411723602022512, -0.454977398955907, -0.444371010673532, -
0.422930806587484,
-0.407809257077638, -0.42665889906324, -0.439836785887343, -0.39462640277382,

-0.456506177242119, -0.417786390221254, -0.3945654618087, -0.380526113175225,

-0.435803862954649, -0.408673249901993, -0.443310483445116, -
0.435285529615121,
-0.391246171534911, -0.423719854507996, -0.408551661568188, -
0.450754029735524,
-0.472602255523911, -0.508979053023968, -0.48467115752141, -
0.483538980395181,
-0.48117952852397, -0.527682963662319, -0.530048158756582, -
0.516321493054729,
-0.505999169443683, -0.481713418746599, -0.509650229070124, -
0.497698126365811,
-0.498898861303871, -0.513440145657546, -0.506018315131867, -
0.513210881839601,
-0.517653436271622, -0.517966483189416, -0.493197050249104, -
0.457616082862051,
-0.445730060732531, -0.445076775112659, -0.418087113066591, -
0.36276831926877,
-0.375690724221468, -0.377881676382754, -0.371887783371946, -
0.425624891596683,
-0.381690010011701, -0.364221504087779, -0.382256197499271, -
0.427591005336707,
-0.400395407899293, -0.384022614564292, -0.418485018168370, -
0.439423306152362,
-0.387746741397525, -0.419803656308417, -0.446329175681457, -
0.485574391008224,
-0.494700574968507, -0.476085998563161, -0.485009731706272, -
0.486608912431453,
-0.446653718022881, -0.449368213900281, -0.448383989719459, -
0.448038045747502,
-0.446146531320324, -0.412949964719027, -0.409053112488588, -
0.381923754371735,
-0.373175213060887, -0.37981480418489, -0.418067908641623, -0.45013561487107,

-0.449711660269272, -0.439274610780474, -0.385770827914678, -
0.372747104530212,
-0.386867708368611, -0.388172052642435, -0.357844117340447, -
0.378325067585380,
-0.367948775430377, -0.356183512307615, -0.368755500006355, -
0.347651383123023,
-0.335184483613088, -0.329242820022007, -0.321073285008462, -
0.340146136324349,
-0.311831319422315, -0.320217397857061, -0.356625496546983, -
0.372687092716155,
-0.304876020028656, -0.313825909214317, -0.258608442215043, -
0.180512658085543,
-0.166410400807192, -0.210524515575469, -0.247866669513052, -
0.263465896656237,
-0.304407304525663, -0.294349034314515, -0.286981852571012, -
0.289152655728849,
-0.222120657232285, -0.201859415849224, -0.224396781073319, -
0.252214854031049,
-0.241481461911827, -0.254155560259281, -0.234920660548415, -
0.266435837149800,
-0.295460049739560, -0.301573757633749, -0.293580836871516, -
0.286929047556988,
-0.277698455775873, -0.273735769220511, -0.300415449827285, -
0.32491594248558,
-0.323226612941381, -0.324478124714544, -0.343944391275079, -
0.332753931098016,
-0.354544735052724, -0.351136814924285, -0.443972485195655, -
0.480378162322601,
-0.456466271168926, -0.490811357473404, -0.495228423726588, -
0.473875567999266,
-0.461202095397426, -0.428236514267841, -0.43316660778517, -
0.441026820386022,
-0.391501536306998, -0.383834931520124, -0.446076138382415, -
0.471862767230103,
-0.484409089668503, -0.512638847721949, -0.57615753476778, -
0.612910121929396,
-0.59409101483892, -0.576073124820514, -0.558917860644171, -0.51741302487145,

-0.444874701848801, -0.452634836362976, -0.501076461133736, -
0.506139148830469,
-0.477096807758099, -0.47903948936219, -0.474304789092253, -
0.519983954579823,
-0.536006287118978, -0.533221006407308, -0.530703365114487, -
0.489301276239221,
-0.482504240833116, -0.482672471692736, -0.496699054247686, -
0.486365502289346,
-0.478302125166766, -0.41699619892767, -0.419811322713434, -
0.461764943597718,
-0.444944325641512, -0.482067207681564, -0.453787868299434, -
0.445862832779393,
-0.500514950852453, -0.54597180181228, -0.59370469930558, -0.622642218811955,

-0.654856796909716, -0.662951487359358, -0.628950615954113, -
0.639558958447642,
-0.68318330922208, -0.696308156898031, -0.672057658680688, -
0.661777110541341,
-0.60490005242343, -0.604182417214347, -0.592238365842442, -0.56106239727129,

-0.604686121506906, -0.56831788795367, -0.575747374999608, -
0.576470254631367,
-0.599655515108991, -0.576602039304853, -0.620742702835376, -
0.623936272446067,
-0.678729593640248, -0.700089982580816, -0.724519962982749, -
0.744258196887065,
-0.77249579598234, -0.767772880386214, -0.830516027924284, -
0.780768413310955,
-0.770728227679768, -0.807497232811467, -0.879388092402374, -
0.877323508209862,
-0.90594146000351, -0.915337533037419, -0.960077494402497, -
0.996453875810982,
-1.00171935025677, -0.971945793840293, -1.02385322331957, -1.03545216489193,

-1.03042547742818, -1.05688598006666, -1.05132715324897, -1.03923643510103,
-1.04441463738159, -1.02606950821931, -1.02207418726124, -1.09384167822686,
-1.10759852188960, -1.03022988278253, -0.979906979321654, -0.98140738601237,

-0.955516498832963, -0.939980354153678, -0.975803371605102, -
0.978220217918253,
-0.91696505099489, -0.94341626578855, -0.89133658495402, -0.855700592049311,

-0.895033078919327, -0.89627932177841, -0.938753796944238, -
0.950927820337492,
-0.96380837671684, -0.974038158542758, -1.01277050185216, -1.02419954622913,

-1.08846210248524, -1.10659736422949, -1.09389877696098, -1.12153568516753,
-1.09295618980445, -1.07529181899901, -1.03767165603343, -1.02442281832222,
-1.07439802122888, -1.05377483017375, -1.0485017136285, -1.05988777343173,
-1.05902132289274, -1.00325932676755, -0.97680173673455, -0.974825258532387,

-1.01655342958140, -0.977698033964966, -0.92760046211231, -0.955327438302301,

-0.99713716544251, -0.971561439917154, -0.980895949333255, -
0.968004870171057,
-0.97614953144124, -1.03366215332186, -1.04409457003187, -1.05158861021004,
-1.09196993450594, -1.11264011244533, -1.06717708434083, -1.06480370833963,
-1.08575759875714, -1.09341227228900, -1.09177567504520, -1.06910800772099,
-1.08712926865485, -1.09686409752487, -1.09774736862532, -1.11213772549334,
-1.15403915449551, -1.22373850734154, -1.21648676898269, -1.22018337968052,
-1.24952942783237, -1.27656503710055, -1.25972114225839, -1.24703179692361,
-1.20919050916248, -1.1842207616518, -1.19287391961853, -1.19721171380420,
-1.19501363104455, -1.18939215272863, -1.14678652882893, -1.15465661360344,
-1.11813358574639, -1.10040843911149, -1.16032096536330, -1.14919807942398,
-1.19048875998298, -1.19989772675020, -1.20398620784167, -1.26586491882295,
-1.25851869647003, -1.26544165917516, -1.22693569033306, -1.25725659511552,
-1.21010427615498, -1.34232359165001, -1.32660610047776, -1.29610987865021,
-1.26013011718281, -1.21766261392912, -1.17439826445927, -1.17622836607061,
-1.19410480592852, -1.20980753722721, -1.23210362227815, -1.22356037273062,
-1.28685023169650, -1.31147384975096, -1.31325943087577, -1.32294429181585,
-1.27764956444795, -1.28185065227558, -1.30324581243296, -1.35720257804154,
-1.31485149250653, -1.31249005177055, -1.33192918666462, -1.32949449292688,
-1.34651906053735, -1.37881888127112, -1.36863148386580, -1.35635968160029,
-1.31362991720627, -1.33018045547167, -1.31016961130448, -1.28971757113766,
-1.27303302564893, -1.31791694333118, -1.30060097705289, -1.31787074800653,
-1.29868914777550, -1.30857532844143, -1.28573219864486, -1.32410102371601,
-1.33639969073421, -1.34007415835075, -1.38507735502386, -1.34107837659624,
-1.32257301638253, -1.31159626761432, -1.30609936880576, -1.29975894236571,
-1.27624841277619, -1.25405147690197, -1.22910085941265, -1.20795472697643,
-1.24853851088728, -1.27438582338932, -1.26400606495224, -1.23302999402678,
-1.23178880710196, -1.24518199542743, -1.23331612167119, -1.22903127206189,
-1.30490842217184, -1.28576361752545, -1.31074792460780, -1.30894341383401,
-1.25965302705916, -1.24232528469753, -1.19631320277742, -1.18300183600512,
-1.17795022721718, -1.17403776779038, -1.12859939113535, -1.17956095926659,
-1.14352875660976, -1.15024124231263, -1.12826553775509, -1.13025738690559,
-1.13484753813885, -1.13056767812354, -1.15714576109888, -1.14068343396429,
-1.12261477644592, -1.09392813920472, -1.13016278463302, -1.12804293615061,
-1.15706655799472, -1.19894078640923, -1.13702840262820, -1.13987336894122,
-1.08890542138072, -1.07906342544466, -1.09459214157365, -1.10581794169538,
-1.05652075609808, -1.09629145098239, -1.06052971723101, -1.03095393415127,
-1.03757475536182, -1.04433286723235, -1.04144782691406, -1.05020397226132,
-1.11335366423297, -1.13308242788617, -1.11881260434659, -1.10274785767592,
-1.10409688950570, -1.07848493782634, -1.05652152591989, -1.10303694008154,
-1.13047872199259, -1.11449882444427, -1.10209602713833, -1.08728702613152,
-1.08016151023848, -1.14207794356016, -1.14108708993475, -1.11437528006918,
-1.13074737179464, -1.10080001358879, -1.13380131056501, -1.02974340003599,
-1.03118429349646, -1.01159974250007, -1.01706186526495, -1.01352921118355,
-0.990699389037849, -1.03145411246970, -1.08279225328073, -1.08031351179812,

-1.12183163205751, -1.13694519045933, -1.10873154083012, -1.09048086125836,
-1.10496312564202, -1.12271481183688, -1.14404655640216, -1.10727842881992,
-1.09827438761092, -1.08984538074143, -1.02245675642832, -0.982400792644428,

-1.00758031375945, -0.997985602590168, -1.01785304453019, -1.01973792693547,

-1.05359645335894, -1.04117220263969, -1.00443264885732, -0.94462766276219,
-0.93052780914025, -0.901150831889696, -0.876218739372558, -
0.889781755283235,
-0.922689983868312, -0.986812563893075, -1.00577437687499, -1.05528836194302,

-1.06361757010428, -1.05509423969802, -1.07764076165386, -1.05393419503949,
-1.03141748838034, -1.07666158378067, -1.08826736188342, -1.03825851781967,
-1.05838978627658, -1.05513301624640, -1.10147358321021, -1.09601673752277,
-1.10974357792115, -1.10308596910112, -1.13992099599508, -1.15495903878899,
-1.12672669661503, -1.15364353419803, -1.14456118003669, -1.13725384673776,
-1.12789761787739, -1.08972321956154, -1.06817521360676, -1.00612585643724,
-1.01721188607368, -1.00086920668929, -1.01165701637624, -0.996640056212344,

-0.98313986002211, -0.984299091792497, -0.96583710325524, -1.00232690455114,

-0.982922825052307, -1.02459712149803, -1.07881579816258, -1.03163543147277,

-1.07729789184894, -1.01180574553278, -1.06293062885976, -1.04345761615412,
-1.04507118212109, -1.06662413792108, -1.10636689376820, -1.12222320840262,
-1.16318956376199, -1.22192850138826, -1.24348755430912, -1.25220663040648,
-1.22215708457057, -1.25236192402018, -1.22111129025221, -1.24724680332852,
-1.22611890659054, -1.23082404646155, -1.21893054743993, -1.22794164100034,
-1.24178410222557, -1.23491411880606, -1.21518417771487, -1.17966072752994,
-1.15118328785644, -1.11210377352095, -1.09023084325471, -1.11849942190869,
-1.08675131543657, -1.09597339581439, -1.10809871281075, -1.10221229540946,
-1.07320636516498, -1.12163406259843, -1.12479846098100, -1.11007095163810,
-1.18712453125076, -1.17536343414050, -1.18780217004682, -1.18605983953425,
-1.18095143598242, -1.18140299327442, -1.20254462890750, -1.18702354604552,
-1.18338989324559, -1.25961943500943, -1.26207336479591, -1.27358337177702,
-1.28871107136099, -1.23882265409051, -1.25505293124848, -1.18747335150791,
-1.22188479870869, -1.22512921075359, -1.27267669879780, -1.29897687823820,
-1.30983284713907, -1.31092491248901, -1.30162429691371, -1.33581710215692,
-1.37449534029653, -1.3624275416685, -1.38856284059204, -1.37690066885045,
-1.35812270817097, -1.26439412795243, -1.29012098316903, -1.30021268606416,
-1.29777924276391, -1.3330372548974, -1.29014528078931, -1.27029202224966,
-1.28153308247510, -1.27170307834150, -1.29079391277124, -1.27810834543040,
-1.29612763984711, -1.29674584868698, -1.26923823763037, -1.27642570501084,
-1.23329549212125, -1.24558112749617, -1.21620176421750, -1.23394969740689,
-1.21143465389589, -1.19397690475935, -1.20538829239449, -1.22385912544701,
-1.22335448659759, -1.18764026930546, -1.18460739308550, -1.18544657910073,
-1.17488653166367, -1.21217413413652, -1.24804440395677, -1.27705906222957,
-1.24670841596142, -1.22433660324650, -1.14970333463809, -1.16375279138769,
-1.21568964613097, -1.21945917869375, -1.22821658985229, -1.19934594306181,
-1.14643724079028, -1.12643944455629, -1.12233497987804, -1.07034913586819,
-1.05097392060227, -1.1717978905894, -1.16695008848828, -1.20329793233844,
-1.22449878968466, -1.23620777156718, -1.29283095807763, -1.30456841492716,
-1.29956251426570, -1.35321339995369, -1.30781507790772, -1.25820955708785,
-1.20144277085079, -1.20233755352439, -1.1957463134694, -1.24548177598815,
-1.23362870938309, -1.23751123003283, -1.16881036271405, -1.17867785852131,
-1.13299980838684, -1.13575190845236, -1.15265011862258, -1.13951844565801,
-1.20445908347423, -1.19738791438949, -1.17891493280400, -1.14633612729150,
-1.18945230462923, -1.25911000157446, -1.21092734531082, -1.19892430224433,
-1.20094139684876, -1.20934910682816, -1.19981222813604, -1.20380034427012,
-1.24216556039083, -1.27832948787905, -1.29714241228392, -1.24544589067096,
-1.22418764631655, -1.1588736470015, -1.13559976408933, -1.14546273748605,
-1.07796573614627, -1.08111480654649, -1.05358267953188, -1.08838275274950,
-1.06160134081697, -1.03782319420854, -1.01718992134545, -0.98004212412867,
-0.995178761325842, -1.03019373420257, -1.05128326974742, -1.05013618518556,

-1.09826884471200, -1.09809384521198, -1.08519850594646, -1.01064238573864,
-1.03473773757524, -1.04378027348155, -1.06594802906565, -1.05967833110358,
-1.09239793955592, -1.10372033897668, -1.11398074918965, -1.07142615682281,
-1.04147669240454, -1.02369590981874, -1.00377810316095, -0.992689744428071,

-1.0585698310738, -1.08344112379608, -1.10562560362484, -1.09246302856515,
-1.08639186584488, -1.10794139071165, -1.08098657419634, -1.03139971026396,
-1.03491525139187, -1.00675304495887, -1.0124392044228, -0.992602075541911
))



More information about the R-help mailing list