[R] Fractional brownian motion

RABINOVITCH Peter Peter.Rabinovitch at alcatel-lucent.com
Wed Jan 10 15:42:00 CET 2007


You have applied the Hurst parameter estimating function to the fBm, not
the fGn as described in the help page for Long Range Dependence
Modelling in fSeries.

"The function aggvarFit computes the Hurst exponent from the variance of
an aggregated FGN or FARIMA time series process"

Trying perFit(diff(x)) (where x is your fBm series) gives an estimate of
the Hurst Exponent as 0.53527922 

-----Original Message-----
From: r-help-bounces at stat.math.ethz.ch
[mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of ED
Sent: Wednesday, January 10, 2007 8:53 AM
To: r-help at stat.math.ethz.ch
Subject: [R] Fractional brownian motion

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
))

______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list