Skip to content

Instantly share code, notes, and snippets.

@jebyrnes

jebyrnes/dcov.test

Created October 23, 2014 19:09
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jebyrnes/b50e92e72f0ee1852cc5 to your computer and use it in GitHub Desktop.
Save jebyrnes/b50e92e72f0ee1852cc5 to your computer and use it in GitHub Desktop.
Download ZIP
misbehaving_dcov.Rnw
Raw
dcov.test
\documentclass{article}
\begin{document}
<<echo=FALSE>>=
activityDF <- structure(list(julian_day = 1:365, food = c(92.6182247978374,
199.361820535986, 128.619308542186, 131.048504694701, 44.6344539529983,
139.02542881751, 118.815437972886, 85.5963235153311, 73.9459098960665,
85.6671795660157, 195.954639346317, 192.942941724749, 193.341526021202,
142.419581610027, 130.23002915561, 232.843857188757, 105.14817987285,
131.58751463704, 123.30056211349, 150.447279408636, 138.40775613945,
218.648860534883, 194.421101663946, 184.27966957954, 45.5636832866875,
194.751288703184, 72.3476195748857, 91.8024365506414, 117.912307182212,
152.356152920377, 123.774505920761, 210.431311364469, 92.4850481923936,
83.4929918472923, 140.89057603969, 212.185010689255, 110.89077098351,
158.915635073238, 243.602462842277, 113.85151636109, 187.461965135917,
123.39807077405, 109.366896172138, 62.750055759989, 158.27572933655,
183.643527951499, 118.062641750957, 163.809693408412, 161.470046830902,
44.9620261936814, 75.8956856397685, 145.374340857118, 84.6335256124607,
96.2343370636656, 96.4908079226898, 157.672272981851, 103.053906280901,
208.941405091109, 138.16950904558, 203.165328487543, 208.005235171877,
176.090096993895, 124.939457574059, 123.578722676302, 210.158020974461,
123.701611792461, 88.1279901849868, 178.172167248393, 115.485200534093,
48.1753644223312, 114.437008847433, 216.890600885956, 160.577452694601,
187.636720159085, 35.4104915656917, 93.1159396769545, 184.552912326746,
155.851153652582, 53.9315864058179, 67.4922411996844, 177.731419240215,
161.022137995797, 177.751646725712, 66.5464744721085, 114.611506843323,
182.11848029471, 167.665754880109, 184.261466089455, 98.1500535371342,
141.557492252183, 146.407447374401, 206.002324389895, 108.687087822477,
156.134825622264, 150.818575370124, 176.64496458315, 66.8279889121797,
111.351363515011, 29.7490965154032, 140.984820272781, 107.266407825426,
129.711537146045, 138.585158027869, 106.218682958644, 102.364474480861,
171.707832406573, 145.469685071349, 63.8891191439936, 206.378175234226,
219.947568084034, 148.530738273619, 128.978722777434, 85.1779466864904,
91.5788554671927, 131.809777122641, 133.972433951106, 104.742488588141,
94.8160253943713, 94.8061837420664, 47.2254266613728, 159.966268353496,
158.680285268277, 91.6868871824241, 165.087107860318, 100.371924681926,
131.861695167093, 107.377570698527, 173.176540687664, 119.916807158549,
175.84740313946, 148.889936186203, 118.940375372028, 125.850615964101,
92.0543617414466, 133.508587619841, 144.456864225918, 105.51642023199,
272.998243984191, 56.6593421796523, 135.863931666786, 104.783127844088,
141.198630680474, 101.573381155845, 79.6706587566052, 30.1688376061025,
86.9480235828328, 39.3584310443109, 127.954860731941, 174.80679302295,
63.1647831008937, 55.5521916103169, 109.26569417767, 95.1007059548783,
113.728428017884, 98.5464837770939, 130.209634320759, 120.539447047208,
50.3800756706825, 120.790433070926, 52.2447785778057, 95.9734917332465,
108.113289822878, 120.905723138996, 127.210004315811, 114.10274979741,
163.036756058898, 106.899123936635, 146.251101437674, 39.5692889581971,
165.081200732599, 67.9882618671799, 76.8936966368098, 165.180352624766,
83.5460525341158, 126.859662498911, 116.001986704702, 226.654829230487,
116.805185799183, 150.158616578881, 109.60194484322, 138.22405923649,
36.4872550696884, 145.714902046142, 188.889564444867, 124.077543440467,
133.50037197129, 168.342067008401, 136.293888202937, 43.2729896598689,
115.615200466435, 121.459532145185, 79.490813746232, 60.2907034543131,
49.2959436668144, 88.9400335439815, 176.241030989141, 104.258356012765,
111.891459341663, 80.9619026337839, 84.9492019336249, 210.432170822072,
0.0505059285401686, 59.9886859678684, 135.769230954418, 161.686164685214,
53.2285882381232, 131.564427951866, 145.69854380207, 67.3821122354483,
136.713945420718, 154.34070211432, 119.024284113464, 108.536388202401,
87.7333333717414, 81.1017571536505, 122.454993384285, 81.6826716057264,
90.842087903626, 137.217049025757, 77.7069052204681, 131.275702287324,
134.571308044248, 138.837453965803, 115.037018184555, 19.6187559069408,
117.344248072864, 107.247561554318, 56.5311440599193, 60.1230860733006,
88.6946864266355, 140.455748933063, 121.076124021638, 13.0042748654361,
109.554059373576, 203.966997251439, 105.868501831117, 165.112278909329,
162.687576585181, 62.6109733141616, 82.0072428216217, 126.436842863547,
128.749345279836, 26.0735930142839, 65.3027640679725, 106.653544529857,
76.1501211329605, 70.3916992631395, 141.745428973418, 104.418747241088,
135.767604732908, 166.160771286207, 91.2317356723762, 207.908308330388,
73.4356928257744, 92.3686480336286, 74.8653656581363, 132.926730915638,
79.1355036712226, 140.742358185361, 98.3589025888085, 25.0322947731106,
177.057842724031, 137.611119437581, 219.468077255829, 179.320624766792,
87.7064039806991, 74.8696951970104, 76.9263674398227, 140.320366002433,
120.503377133764, 101.658338269288, 120.672012014565, 118.133786664406,
135.454259196793, 181.757496925907, 109.255179399929, 74.9745432858777,
104.158812435609, 62.16128549595, 170.872467830164, 135.697225606721,
152.996255038669, 110.330340853093, 113.692766555864, 148.337081198471,
90.8588359172082, 125.331430999266, 149.017166417168, 118.498082040037,
162.4879453457, 137.829938308553, 20.9188316277781, 82.1302578490659,
102.942717568631, 188.409558657804, 109.087968360525, 170.3730016912,
69.6434022655983, 124.960082493789, 77.4454998751145, 137.098599142723,
183.508794923258, 139.112977180323, 155.272469076326, 130.259092844694,
140.079415888445, 175.818619302811, 15.672605522213, 186.430407139864,
241.291900090408, 168.616695783813, 145.854112461188, 223.599087704853,
121.795764613106, 181.868222898767, 102.834379921068, 288.786910812126,
59.76730875369, 127.839451469262, 131.29992454163, 181.465832158727,
78.9908288037463, 61.7818907288452, 66.2868999265771, 122.218941481859,
122.779973563298, 102.422052395125, 189.12726909292, 175.181018874536,
109.527435358322, 85.30231949349, 130.154540869946, 126.714087863695,
102.600730020963, 106.134725094268, 231.583754568073, 182.441499208991,
150.800919209349, 159.947295055341, 187.172018629861, 116.412725194469,
72.3613765338995, 274.011822953084, 108.371416182906, 128.796904715756,
167.667985005819, 324.649464017018, 111.966809961498, 106.261888481444,
196.257954588932, 225.254470194564, 127.60267701444, 176.052138905108,
188.416932766027, 86.907331167983, 123.355324982167, 202.013292949984,
115.584069953517, 283.862724847073, 168.241748242238, 148.073617394392,
193.650363685957, 150.541826260407, 136.292118070661, 103.718550991133
), activity = c(11.920097538528, 12.0195990013862, 11.1807927813684,
13.6595810456339, 8.79081163944262, 13.5467691746706, 16.2972410245348,
9.64664120641838, 9.41749570906209, 10.2533466620528, 9.7946806929984,
13.8676095942115, 12.6687121034226, 13.819277584983, 12.7079897846093,
25.8694817855055, 6.70457421134495, 17.0507544561441, 10.6684396335407,
19.2286971138972, 8.90868495997642, 22.7114071168936, 10.5134538562663,
12.5746159259715, 14.2857589351014, 16.2780722812274, 11.2171662900959,
16.3620201812271, 19.0007011619652, 5.57892250169615, 12.3289771053315,
27.667008083528, 17.6917895504001, 10.3714141381093, 12.6055277285178,
22.4063094853805, 23.4152391089365, 14.5232812211249, 23.6905028022803,
13.885034076517, 14.1034536681749, 17.5392662747931, 20.607723435228,
16.8038524533458, 18.6172601112022, 17.7581539602919, 21.6803546324845,
16.5356535179442, 18.8760228507647, 11.5369787073068, 22.6866958509062,
23.324724164803, 11.3249270716243, 14.8103597281481, 17.8516176684395,
20.8920078506492, 17.7380258297224, 17.8564594141057, 11.8812179200784,
20.1145086883129, 30.308887383542, 21.5315628321876, 16.8524044606184,
15.63015656438, 27.4799963849614, 14.4371539348239, 15.2722882892596,
16.4018684809308, 20.2887401218192, 16.3891978201985, 20.5586975217359,
21.8190728726519, 21.8623855905704, 16.0797197742758, 5.35293082904277,
12.2438869334952, 20.3238799931006, 23.6223673369945, 15.2702035177931,
11.0843896306172, 30.6145007847849, 19.1356587453439, 16.322235075022,
18.8896127526115, 15.8849303368809, 36.625333825011, 27.0174271539914,
13.740261760329, 17.5743382748244, 16.8911276348754, 15.4122667059764,
20.1777482137377, 13.7293803111715, 14.7156009844141, 20.0049765785071,
26.5393039062544, 4.19761479499594, 12.5683712144758, 8.95986239538307,
17.0676780283871, 14.0352600425057, 18.9037869896361, 21.3236987630101,
12.1305933568469, 16.8559814397427, 14.470045898631, 20.8790125553639,
16.481901190716, 23.4000242318959, 30.9566264496972, 18.9602032653072,
15.2307637412031, 9.64994793226168, 13.8789860629882, 13.3163585986214,
8.12332360968242, 7.27384584588638, 16.0143963854747, 14.6589266027771,
2.68660199656501, 22.2801893225556, 19.7048486886616, 7.98037601903306,
11.4583247237043, 16.7757371254774, 17.0721865184138, 12.380913520529,
11.4916630564651, 11.2737209219238, 14.5988040390794, 10.9404700644643,
12.3890612985543, 19.4704466200111, 15.4399988603335, 11.7492796452933,
13.8313412133427, 6.96146405877394, 14.4140770202713, 10.8283408026079,
15.1829482524231, 16.6560133699902, 16.8392587590229, 16.0494303065846,
13.5317640751209, 14.2200874242571, 10.8177541706431, 5.84400350024724,
12.5912668364878, 17.0044879012595, 11.7496520498878, 8.47256554481267,
8.01523672437145, 12.1850711617907, 13.0410983954508, 13.2854523839741,
9.44664116749571, 6.96125278661691, 9.84211804738037, 11.0537270610413,
10.2057999988903, 11.1003508546259, 15.3126973957463, 10.2343320088724,
11.0034655540186, 5.39010009999515, 9.11054909920172, 12.924843639867,
9.91913344974142, 7.66593668857548, 10.8076156323269, 10.6107197581744,
7.20438776908952, 14.1081954150818, 5.63830826588118, 12.043364292193,
12.9924377281289, 12.8822950316483, 16.9980014489642, 11.6488607321872,
9.81284624648331, 12.6407949406939, 10.1482591671917, 15.8590916956152,
12.1539587032845, 11.0057305120311, 10.4456096296013, 15.7003838452032,
16.8815627387071, 7.07407145108536, 7.15717261098184, 13.8044480483327,
15.2572517828595, 3.40253150095381, 6.19080719599411, 8.95091128718426,
12.3036294245981, 9.78182170752043, 14.3313083053085, 8.68371044161526,
9.71718780182442, 13.4940279017968, 6.52429862120466, 10.2171629284111,
9.06677022310845, 10.3176238052634, 10.0799514043202, 13.2812663163393,
12.2000716283678, 10.3729171995664, 14.3963918619489, 13.088300206672,
15.6944450162824, 13.2895114566821, 10.8487790214636, 15.488046268643,
13.3912247770387, 13.0115202141572, 9.09935503704645, 11.632013251895,
8.91454313086708, 11.0549500809369, 8.93601778679673, 12.120908506501,
12.7578506417873, 10.3275340621882, 9.07141686490589, 12.9656604200678,
11.3823805313684, 4.31110374341891, 7.90579204208663, 11.9482355174396,
8.17313612678552, 5.79651170525914, 15.1355794090919, 10.9429423642176,
10.0323390651674, 15.1179222228796, 12.3236364005908, 5.61311551278944,
10.1193594255859, 10.0805354547339, 5.38167967875483, 11.5115140121868,
7.2524630056032, 11.6862740981679, 11.791498166506, 12.0092392215093,
8.44761736648053, 11.513655245768, 7.13642926891478, 16.8037097294941,
8.61484713185853, 13.7063664053723, 9.77298212613597, 13.5960648552787,
9.25125795985577, 7.71493241790108, 10.0027694294737, 14.0180710443908,
8.86127381633078, 6.24693471524571, 9.91682634642715, 18.267824647319,
15.8071087311249, 14.2839553107402, 9.83206277565418, 14.0359619293604,
16.1920115100821, 12.52463518731, 10.9821801853018, 4.60633433318818,
9.46902047448621, 10.8754860160782, 13.3666637274135, 10.7216876463376,
9.21223370097406, 11.943629129986, 10.5235734548516, 6.06410370911916,
15.577080910686, 14.8836467584441, 15.6688166039721, 9.54169072089836,
11.8674238756015, 16.310578424976, 15.326289426417, 7.01938169078748,
12.4913861569567, 19.7719745502327, 12.2114568633856, 12.6263437023503,
9.33273170240612, 19.1605100881563, 15.4857092382258, 18.9976165805992,
19.1913725411898, 17.341384400225, 11.7901119471321, 12.8504018687935,
7.71710163568595, 10.5522591456909, 16.9590875131269, 12.6196306545294,
13.4013150232804, 7.84670944772032, 11.9408665270102, 11.3328199123383,
11.0554922652021, 10.6737510486332, 26.1390112233044, 13.7207289932775,
19.7972250012392, 23.905684165791, 13.184129614166, 17.3582604252191,
15.5391525226031, 25.7654217864843, 13.1676984677302, 12.7481869619464,
14.9580845779175, 29.476317251101, 13.6114536356502, 16.9024974804217,
13.5649459465052, 11.6107711845281, 9.79654793422858, 13.9479485694661,
15.9751845741874, 32.9831119347267, 13.9782623776588, 18.4147198534295,
19.5441687289734, 15.7174847344547, 13.9166119288116, 15.0624394935599,
29.6996950305479, 24.0778927794494, 11.6990163009617, 17.9522072735396,
24.7603354872186, 11.7556555582864, 13.1089013028834, 38.0504659705217,
9.68137043403555, 23.9984042062278, 24.7607279183477, 53.7052183947374,
16.9048945666354, 15.2099412517356, 19.8423406648329, 20.0459505834903,
20.0794759371227, 16.1826511824307, 18.0905753872859, 20.1929410763538,
16.4990521573372, 24.8321991237985, 10.6706217483247, 34.4648984041675,
23.200796261837, 14.418704908135, 12.7128176352372, 19.834403356468,
30.6737126793986, 18.1085932426884), temperature = c(6.35507454696068,
11.7560310188448, 11.2495939621507, 3.73346103334718, 10.3661182111891,
10.3085866824344, 12.4775058121735, 7.66655437752012, 6.52970733076538,
-1.87605477385517, 5.10449558666635, -3.58753342645542, 4.11697274513642,
7.55052878453628, 9.71727905677056, 18.3816666805458, 3.40597222774437,
12.3284594849146, 5.85670971678146, 13.9990541297326, 2.61658736772251,
17.2757499465976, 8.20670474820072, 10.6436211167783, 3.8105542944817,
9.88203582957166, -1.01891147435467, 17.1171969067754, 18.8112766130188,
5.51855073235419, 10.0823185080486, 19.6336430299285, 14.6519704513102,
8.76312278805939, 9.08112418367295, 18.7597283043102, 14.4754080807138,
11.9579813415786, 18.769827644851, 7.82361835423932, 10.5921686559314,
14.815221838897, 15.2389453501096, 7.0346160592125, 13.8826118276844,
13.1989011556973, 15.848460083646, 15.2229349493246, 16.745582953479,
10.077652053366, 11.701056319025, 18.2001760021992, -1.02299811407325,
7.53487637489243, 10.4604602337602, 16.4495919711491, 5.74980334489716,
15.657071950732, 6.86834725920662, 17.2712233929732, 21.2567043876971,
15.5905097246646, 12.2061104174593, 10.7613208974803, 21.0105558821797,
13.1128164423221, 12.724058975623, 11.3453520127134, 16.6644948888331,
11.2046012710649, 15.5368316514543, 19.1606968100224, 18.8000303633099,
15.3822905559617, 1.69447335492661, 12.8052277664745, 16.6961136895892,
20.4878052795757, 11.0084263348347, 4.07921618364478, 19.7649414364788,
16.5319391572095, 11.4427983074454, 15.3044843072756, 11.5916056989608,
21.123439052776, 18.5902156372037, 11.1387496052556, 12.2963507043683,
13.8434033244688, 6.02042902035332, 18.0550526952859, 11.4393155728915,
8.93540712364661, 17.0318537356545, 16.7245382753348, 7.83824206075789,
7.28725778160235, 7.80902414538936, 11.4239236297597, 3.69670162323183,
13.9593548128048, 15.2145272068447, 11.8144549043748, 13.9397749203779,
10.6619643369176, 10.652365588167, 13.1135015334356, 20.1695164895092,
22.0007190659263, 13.0790836329263, 10.5128803520964, 4.75360345709931,
13.9802442798788, 10.57415191661, 0.993505161328855, 9.21837225516134,
11.7117762834927, 7.53439048471802, 6.57477166198307, 17.4507990248835,
10.8297962144046, -1.60288173509049, 6.11950415303803, 12.7243528351451,
6.12965368564584, 10.8622984779187, 3.62810563247419, 10.9106692753153,
1.25912752789464, 3.98359173610398, 0.649749324726187, 13.9272789456756,
-0.360229254947622, 5.29188144929541, 6.53636655266501, 1.77561532429037,
7.9863341711694, -1.29254103768951, 7.88079021317134, 9.3737811680975,
-0.270891641997307, 11.4874215919484, 0.842744662673051, -1.59025763579771,
-0.139966659092431, 0.288418493257311, 5.57911704530901, 2.77937323848888,
-4.05792659341121, -9.93315610826954, 4.02977064498619, 0.454413904605345,
-4.12016451807492, -0.985672791844681, 4.01817532061821, -1.26771707532113,
0.348129202921482, 1.65048137080458, -0.946973786554503, 9.34211973703921,
0.943745404457393, 1.90774718934274, 2.4201373254997, -3.20922284112533,
2.15445000340291, 1.05831235602602, -0.152013396463132, 0.920469669145042,
-3.22195034761052, 7.3617554008298, -3.0596864152631, -0.787914564987367,
1.11195217714569, -8.18408834881589, -1.81254507405036, -1.31700994538615,
3.05937611049721, 3.21866262012074, -6.69956721304698, -1.39128827049852,
-5.00495856084543, 1.47330485654575, -11.3124750889325, 3.40930531848876,
-8.44061397748895, -1.19771984624843, -0.553653502647532, -8.56097641265537,
-7.93704939343875, -1.14928573189326, 5.237883501835, -1.37788580850199,
-7.70892950512333, -1.78252413118413, -9.05582996077794, -5.55904373470637,
-12.2155064810008, -0.43372593388759, -0.476701034084234, -2.35341342128998,
-7.29360453685343, -5.15938007593694, -0.566890048838369, -4.24066811182081,
-5.92966943925736, -10.9281260724478, -5.23924201579057, -5.99726548635786,
-9.97692267701527, -0.54089655262869, 2.15866446254868, -4.16089168134964,
-7.45150781795663, -8.96541202442287, -4.30006227539221, -4.6896939342599,
-4.37330505687518, -23.6115909160958, 2.09774858294872, -3.93099526052606,
-3.69285010873705, 1.77266754566051, -19.650558869622, -4.5818510536868,
-8.91735662953401, 5.08143232575424, -6.21266432657428, -2.58822558957274,
-6.27268784073668, -4.91746235819558, -3.83779917252961, -5.18660839054893,
0.287836438126709, 2.96372593666194, 1.61015225709053, -2.37791976479566,
-0.935519553271928, -8.38828377992071, -9.14731306134678, -2.54711435936451,
-4.41428863071396, -1.12918751152106, -3.71855168347905, -1.18196815380449,
-3.06993243667706, -14.6724808578292, 4.03410254430953, -1.092990356009,
-1.9865532121938, 2.48925447920705, -7.81840920346058, -1.67158419013435,
-7.24778769773717, -0.692394718336579, -1.24364051651762, 3.52558627612833,
0.247800814906269, -3.09610273845469, 5.65979772727581, -1.52105673876609,
5.43562817357103, 2.1677641758387, -0.749757018743299, 5.72306662255618,
3.63721765437149, -0.721087476822461, 9.53814896918926, -1.66001547274739,
3.39420212417525, -6.15352721569824, 4.28594135848916, 4.67071280557993,
6.79935332577831, 4.75448093530858, 2.59805920875356, 6.03392408354104,
-5.21392889607793, 0.313240763969512, 6.43304232276564, 4.15642907672921,
7.28364751186963, -5.76127895579716, 8.05526050602846, 0.849580369627819,
3.96205223486873, 10.1445164629725, 4.14968385758892, 10.1970257589183,
7.94635605534818, 8.61720263455791, -1.43918392210539, 16.2649586852542,
5.76104229437924, 15.6250225956021, 12.533339810259, 12.6990949801351,
6.86694036030026, 9.80525801882815, -1.41015769459808, 12.7844368472141,
12.8245647184544, 4.31285702299309, 4.23875861117211, 7.12165547418165,
7.82930950483037, 7.50369175994354, 8.3693924845056, 4.08454724874286,
19.7837683972037, 11.7884155873084, 17.0462496549376, 18.4714071102731,
4.9885993929714, 13.0661900600271, 11.4031423185502, 21.1149326339772,
9.85678178451254, 12.3878460891814, 8.12126888342424, 17.0356118294437,
11.2456128864316, 12.2962779007745, 12.6366078524336, 7.51315401911126,
6.29964523499395, 12.0784462894698, 15.4992624606988, 22.1155286521687,
10.2500732463639, 10.3588973994376, 18.518525479015, 10.7389621618097,
5.01454219361799, 13.9064643666595, 19.7671191453412, 19.0371132917916,
12.0811696500256, 15.6026283986953, 18.3312031263019, 10.7946273170363,
13.6926232563993, 23.6863126410597, 3.18107395473114, 18.4823669552533,
17.0652908005299, 25.9227271916762, 8.10212413203572, 13.7537822802554,
13.3925712357859, 14.4677764631929, 15.0091045915852, 14.9654683662397,
11.5097228593827, 12.9612547798403, 12.4749236313945, 17.2461154476705,
6.32470350143276, 23.1142704245122, 18.1315648167586, 11.5675048469081,
4.76959787482954, 18.6582520500873, 20.046567968312, 12.8374454402959
)), .Names = c("julian_day", "food", "activity", "temperature"
), class = "data.frame", row.names = c(NA, -365L))
@
In the consolte, the p value of dcov.test is 0.005
<<test>>=
library(energy)
with(activityDF, {
print(dcov.test(julian_day, temperature))
print(dcor(julian_day, temperature))
})
@
\end{document}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment