I'm trying to do a stacked area plot in ggplot2. I seem to have followed the code correctly but the output is different. The code seems to work but the areas aren't exactly filled. Anything wrong with my code? or is it the ggplot2 package? Code below:
ggplot(stacks, aes(x=year, y=remit, fill=polcat3, colour=polcat3) ) + geom_area(aes(fill=polcat3), position='fill')
my plot
My data
year,polcat3,remit
2000,Auto,57131999.97
2000,Auto,139199996.9
2000,Auto,4821840600
2000,Auto,10352872.85
2000,Auto,2852000000
2000,Auto,2161018799
2000,Auto,103622840.9
2000,Auto,39000000
2000,Auto,1075000000
2000,Auto,56885414.12
2000,Auto,180000000
2000,Auto,1340000000
2000,Demo,86343658.45
2000,Demo,87064697.27
2000,Demo,1967528809
2000,Demo,58235450.74
2000,Demo,126870002.7
2000,Demo,1649354980
2000,Demo,26245557.79
2000,Demo,1610122925
2000,Demo,135975250.2
2000,Demo,1838800049
2000,Demo,1322301025
2000,Demo,596200012.2
2000,Demo,641289794.9
2000,Demo,280652252.2
2000,Demo,1190199951
2000,Demo,12883465820
2000,Demo,891700012.2
2000,Demo,49942501.07
2000,Demo,72300003.05
2000,Demo,178600006.1
2000,Demo,7524542969
2000,Demo,73163955.69
2000,Demo,12000000
2000,Demo,177000000
2000,Demo,746886.1938
2000,Demo,9450584.412
2000,Demo,111498191.8
2000,Demo,16399999.62
2000,Demo,6961000000
2000,Demo,278000000
2000,Demo,96000000
2000,Demo,1275199951
2000,Demo,233485305.8
2000,Demo,1765400024
2000,Demo,1696785400
2000,Demo,38100002.29
2000,Demo,4560000000
2000,Demo,33000000
2000,Demo,17000000
2000,Demo,343705993.7
2000,Demo,13300000.19
2000,Semi,597799987.8
2000,Semi,67335105.9
2000,Semi,119317176.8
2000,Semi,790000000
2000,Semi,53158882.14
2000,Semi,273500000
2000,Semi,32400001.53
2000,Semi,536000000
2000,Semi,1845134033
2000,Semi,121802436.8
2000,Semi,537900024.4
2000,Semi,8844514.847
2000,Semi,1165826782
2000,Semi,478087554.9
2000,Semi,342368408.2
2000,Semi,14424637.79
2000,Semi,1391800049
2000,Semi,717663818.4
2000,Semi,7226507.187
2000,Semi,34234287.26
2000,Semi,795953369.1
2000,Semi,8000015.259
2000,Semi,238100006.1
2000,Semi,1288000000
2001,Auto,104125000
2001,Auto,149199996.9
2001,Auto,6538500000
2001,Auto,12119418.14
2001,Auto,2911399902
2001,Auto,3260920166
2001,Auto,117284111
2001,Auto,39011703.49
2001,Auto,1461000000
2001,Auto,52876461.03
2001,Auto,170000000
2001,Auto,1100000000
2001,Demo,189600006.1
2001,Demo,83602699.28
2001,Demo,2104551514
2001,Demo,826200317.4
2001,Demo,135279998.8
2001,Demo,1774780029
2001,Demo,25919340.13
2001,Demo,2056358887
2001,Demo,198443420.4
2001,Demo,1981800049
2001,Demo,1420558105
2001,Demo,45900001.53
2001,Demo,633753601.1
2001,Demo,746951049.8
2001,Demo,296019500.7
2001,Demo,1046000000
2001,Demo,14273018550
2001,Demo,1058300049
2001,Demo,1184989990
2001,Demo,401579681.4
2001,Demo,79150001.53
2001,Demo,112400001.5
2001,Demo,243320007.3
2001,Demo,10146262700
2001,Demo,88167266.85
2001,Demo,25000000
2001,Demo,215000000
2001,Demo,8815530.777
2001,Demo,146985153.2
2001,Demo,73099998.47
2001,Demo,753211425.8
2001,Demo,8769000000
2001,Demo,263799987.8
2001,Demo,116000000
2001,Demo,1402890015
2001,Demo,304678405.8
2001,Demo,1926000000
2001,Demo,1252447510
2001,Demo,40900001.53
2001,Demo,2786000000
2001,Demo,141000000
2001,Demo,6038.0301
2001,Demo,19000000
2001,Demo,297390441.9
2001,Demo,11899999.62
2001,Semi,699299987.8
2001,Semi,49915523.53
2001,Semi,116004821.8
2001,Semi,670000000
2001,Semi,18313718.8
2001,Semi,181300003.1
2001,Semi,682000000
2001,Semi,2011334106
2001,Semi,171272735.6
2001,Semi,550000000
2001,Semi,11113728.52
2001,Semi,746885.9553
2001,Semi,367105255.1
2001,Semi,21972925.19
2001,Semi,1166599976
2001,Semi,5924500.942
2001,Semi,68502540.59
2001,Semi,927148376.5
2001,Semi,15253000.26
2001,Semi,348574798.6
2001,Semi,1294612549
2002,Auto,181707000.7
2002,Auto,140300003.1
2002,Auto,10292581917.288
2002,Auto,2893100098
2002,Auto,204929046.6
2002,Auto,2877152588
2002,Auto,106488113.4
2002,Auto,678485900.9
2002,Auto,39011703.49
2002,Auto,3554000000
2002,Auto,44669322.97
2002,Auto,135000000
2002,Auto,1770000000
2002,Demo,733570007.3
2002,Demo,206630004.9
2002,Demo,76230842.59
2002,Demo,2858057861
2002,Demo,1176951538
2002,Demo,112900001.5
2002,Demo,2448964111
2002,Demo,27150941.85
2002,Demo,2479942627
2002,Demo,250646484.4
2002,Demo,2194300049
2002,Demo,1438022827
2002,Demo,43520000.46
2002,Demo,1600192017
2002,Demo,438744873
2002,Demo,278989471.4
2002,Demo,1258917603
2002,Demo,15735736330
2002,Demo,1260300049
2002,Demo,433000000
2002,Demo,1309079956
2002,Demo,389541656.5
2002,Demo,109248046.9
2002,Demo,138399993.9
2002,Demo,323739990.2
2002,Demo,11029448240
2002,Demo,137653686.5
2002,Demo,56312149.05
2002,Demo,215000000
2002,Demo,7959585.667
2002,Demo,84599998.47
2002,Demo,705382507.3
2002,Demo,9735000000
2002,Demo,201500000
2002,Demo,143000000
2002,Demo,1359354980
2002,Demo,344496490.5
2002,Demo,1953699951
2002,Demo,1379995483
2002,Demo,79099998.47
2002,Demo,1936000000
2002,Demo,209000000
2002,Demo,36053001.4
2002,Demo,19000000
2002,Demo,288378570.6
2002,Demo,12399999.62
2002,Semi,58734134.67
2002,Semi,120092422.5
2002,Semi,1176490.545
2002,Semi,1070000000
2002,Semi,32986892.7
2002,Semi,230453414.9
2002,Semi,851000000
2002,Semi,2143018799
2002,Semi,36716598.51
2002,Semi,843693.614
2002,Semi,435000000
2002,Semi,18965887.07
2002,Semi,1208900024
2002,Semi,5622274.399
2002,Semi,103212654.1
2002,Semi,1070313232
2002,Semi,12000000
2002,Semi,422582519.5
2002,Semi,1294000000
2003,Auto,170957000.7
2003,Auto,222399993.9
2003,Auto,14542468993.9088
2003,Auto,2960899902
2003,Auto,147501953.1
2003,Auto,3613916504
2003,Auto,85281295.78
2003,Auto,771072143.6
2003,Auto,39011703.49
2003,Auto,3964000000
2003,Auto,65225296.02
2003,Auto,889000000
2003,Auto,2100000000
2003,Demo,888748596.2
2003,Demo,273420013.4
2003,Demo,55357849.12
2003,Demo,3191664795
2003,Demo,1718485474
2003,Demo,158166412.4
2003,Demo,2821291016
2003,Demo,38872653.96
2003,Demo,3076018555
2003,Demo,320879547.1
2003,Demo,2325399902
2003,Demo,1633427002
2003,Demo,65100158.69
2003,Demo,2146959961
2003,Demo,516599182.1
2003,Demo,294864593.5
2003,Demo,1488709351
2003,Demo,20999150390
2003,Demo,1398400024
2003,Demo,538000000
2003,Demo,556885986.3
2003,Demo,114833999.6
2003,Demo,173100006.1
2003,Demo,486630004.9
2003,Demo,16653685550
2003,Demo,153798950.2
2003,Demo,128600006.1
2003,Demo,215000000
2003,Demo,12216717.72
2003,Demo,106800003.1
2003,Demo,868500000
2003,Demo,10243000000
2003,Demo,222100006.1
2003,Demo,124000000
2003,Demo,1452994995
2003,Demo,510932281.5
2003,Demo,2122399902
2003,Demo,1607350342
2003,Demo,86800003.05
2003,Demo,729000000
2003,Demo,330000000
2003,Demo,61753292.08
2003,Demo,208000000
2003,Demo,434367492.7
2003,Demo,11800000.19
2003,Semi,53747390.75
2003,Semi,141603500.4
2003,Semi,12388155.94
2003,Semi,1750000000
2003,Semi,46470191.96
2003,Semi,235979049.7
2003,Semi,1178000000
2003,Semi,2201128418
2003,Semi,78157081.6
2003,Semi,1437750000
2003,Semi,571052612.3
2003,Semi,25481748.58
2003,Semi,1062800049
2003,Semi,6931479.454
2003,Semi,148687118.5
2003,Semi,1250253906
2003,Semi,8899999.619
2003,Semi,298785186.8
2003,Semi,1269900024
2003,Semi,36299999.24
2004,Auto,227578002.9
2004,Auto,256700012.2
2004,Auto,19578498944.7746
2004,Auto,3340699951
2004,Auto,1032000000
2004,Auto,165835662.8
2004,Auto,4220835938
2004,Auto,117792945.9
2004,Auto,822608947.8
2004,Auto,39011703.49
2004,Auto,3945000000
2004,Auto,82797576.9
2004,Auto,855000000
2004,Auto,2310000000
2004,Demo,1160672119
2004,Demo,311779998.8
2004,Demo,63013359.07
2004,Demo,3583817139
2004,Demo,1722769531
2004,Demo,210558670
2004,Demo,3575069092
2004,Demo,91521476.75
2004,Demo,3189709961
2004,Demo,319546356.2
2004,Demo,244929977.4
2004,Demo,2501199951
2004,Demo,1838041748
2004,Demo,303235168.5
2004,Demo,82371200.56
2004,Demo,2627500000
2004,Demo,664671203.6
2004,Demo,1716792236
2004,Demo,1866314941
2004,Demo,18750376950
2004,Demo,1623198120
2004,Demo,620000000
2004,Demo,627204895
2004,Demo,324496734.6
2004,Demo,229399993.9
2004,Demo,705239990.2
2004,Demo,19861748050
2004,Demo,155463485.7
2004,Demo,202500000
2004,Demo,215000000
2004,Demo,15236835.48
2004,Demo,59712108.61
2004,Demo,108900001.5
2004,Demo,1132713623
2004,Demo,11471000000
2004,Demo,238199996.9
2004,Demo,132000000
2004,Demo,2495106934
2004,Demo,632895324.7
2004,Demo,2564100098
2004,Demo,1621875977
2004,Demo,86900001.53
2004,Demo,804000000
2004,Demo,411000000
2004,Demo,69896163.94
2004,Demo,143000000
2004,Demo,522837402.3
2004,Demo,11600000.38
2004,Semi,44417327.88
2004,Semi,159198211.7
2004,Semi,14954054.83
2004,Semi,2460000000
2004,Semi,133742813.1
2004,Semi,2330324463
2004,Semi,188672546.4
2004,Semi,58446876.53
2004,Semi,1589569946
2004,Semi,801868408.2
2004,Semi,2272699951
2004,Semi,9526646.614
2004,Semi,179017074.6
2004,Semi,1431351074
2004,Semi,13798418.05
2004,Semi,310532470.7
2004,Semi,1282598999
2004,Semi,48400001.53
2005,Auto,623235000
2005,Auto,198600000
2005,Auto,8831739256.9678
2005,Auto,1032000000
2005,Auto,62021921.3283798
2005,Auto,4589226819.79732
2005,Auto,129488927.894728
2005,Auto,1211823252.50637
2005,Auto,39011703.5110533
2005,Auto,4280000000
2005,Auto,95421869.7624513
2005,Auto,823000000
2005,Auto,3150000000
2005,Demo,1289704315.92613
2005,Demo,432090000
2005,Demo,146880908.452325
2005,Demo,4642385283.80997
2005,Demo,1612912427.00584
2005,Demo,337040056.699114
2005,Demo,2805358000
2005,Demo,117891862.0425
2005,Demo,3345610321.40797
2005,Demo,420346410
2005,Demo,104817181.330673
2005,Demo,2719200000
2005,Demo,2460000000
2005,Demo,446005398.63
2005,Demo,99184576
2005,Demo,3066600000
2005,Demo,692807908.925609
2005,Demo,1912955033.65206
2005,Demo,5419620728.75305
2005,Demo,22125089480.4982
2005,Demo,1761628767.01559
2005,Demo,424991045.929029
2005,Demo,4924310615.23333
2005,Demo,599357676.667103
2005,Demo,534305076.983413
2005,Demo,381000000
2005,Demo,915080000
2005,Demo,22741840918
2005,Demo,177186056.198184
2005,Demo,180360000
2005,Demo,609194.496202642
2005,Demo,22531928.6466736
2005,Demo,17561579.0633268
2005,Demo,66369885.0351597
2005,Demo,129600000
2005,Demo,1440067536.07357
2005,Demo,13732557721.8325
2005,Demo,161300000
2005,Demo,4707910000
2005,Demo,3436552689.38484
2005,Demo,788824890.465872
2005,Demo,3028600000
2005,Demo,1187095814.94597
2005,Demo,92400000
2005,Demo,887000000
2005,Demo,2408000000
2005,Demo,76740390
2005,Demo,148000000
2005,Demo,613577432.378763
2005,Demo,13000000
2005,Semi,56642283.3408211
2005,Semi,163173835.554478
2005,Semi,11375094.5558432
2005,Semi,169999999.999993
2005,Semi,5017300000
2005,Semi,173519405.683879
2005,Semi,2420874471.08604
2005,Semi,313250258.3113
2005,Semi,31854955.6480974
2005,Semi,1975536000
2005,Semi,1116972560.32095
2005,Semi,14640084346.8768
2005,Semi,6866648.84058121
2005,Semi,192500458.801737
2005,Semi,1392672712.79192
2005,Semi,19366300
2005,Semi,321808492.95616
2005,Semi,1282599000
2005,Semi,8900000
2005,Semi,52870000
The geom_area charts are just geom_ribbon charts with the range fixed to 0. Since you have multiple overlapping points in your ribbon chart, I think you need to modify your position tag to 'stack' instead of 'fill'.
You also might need to tell ggplot to provide an order to your raw y variable by using:
ggplot(stacks, aes(x=year, y=order(remit), fill=polcat3, colour=polcat3)
Without seeing your data finding the root cause is difficult
Your problem is you have multiple values for each year, and ggplot doesnt know what to do with them. To fix this we should aggregate our data.
Here I'll use dplyr:
library(dplyr)
stacks2 <-stacks %>% group_by(polcat3, year) %>%
summarise(remit = sum(remit))
Now with our data summed by year, the plot works nicely, using your code (modified to run on stacks2, our summarised data). You can use position = "fill" if you'd like it all to end at 1.:
ggplot(stacks2, aes(x=year, y=remit, fill=polcat3, colour=polcat3) ) +
geom_area(aes(fill=polcat3), position='stack')
Now you have to provide almost every value, which was set by default previously:
ggplot(stacks, aes(x=year, y=order(remit), group=polcat3, fill=polcat3, colour=polcat3)+
geom_area(position='stack')
Pay attention to group variable
Related
Good afternoon. I have a vector 'a' containing 16000 values. I get the descriptive statistics with the help of the following:
library(pastecs)
library(timeDate)
stat.desc(a)
skewness(a)
kurtosis(a)
Especially skewness=-0.5012, kurtosis=420.8073 (1)
Then I build a histogram of my empirical data:
hist(a, col="lightblue", breaks = 140, border="white", main="",
xlab="Value",xlim=c(-0.001,0.001))
After this I try to fit a theoretical distribution to my empirical data. I choose Variance-Gamma distribution and try to get its parameter estimates on my data:
library(VarianceGammma)
a_VG<-vgFit(a)
The parameter estimates are the following:
vgC=-11.7485, sigma=0.4446, theta=11.7193, nu=0.1186 (2)
Further, I create a sample from the Variance-Gamma distribution with the parameters from (2)
and build a histogram of created theoretical values:
VG<-rvg(length(a),vgC=-11.7485,sigma=0.4446,theta=11.7193,nu=0.1186)
hist(VG,breaks=140,col="orange",main="",xlab="Value")
Bu the second histogram differs absolutely from the first (empirical) histogram. Moreover, it is built on the basis of the parameters (2), which I got on the empirical data.
What's wrong with my code? How can I fix it?
P.S. When I type dput(a[abs(a) > 5e-4]) I get:
c(0.000801110480004752, 0.000588162271316861, 0.000555169128569233,
0.000502563410256229, 0.000854633994686438, 0.00593622112246628,
-0.000506168123513007, -0.000502909585836875, 0.000720924373137422,
0.00119141739181039, 0.000548159382141478, -0.000516511318695123,
-0.000744590777740584, 0.000595213912401249, 0.000514055190913965,
-0.000589061375421807, -0.00175392114572581, 0.000745548313668465,
-0.00075910234096277, -0.00059987613053103, 0.000583568488865538,
0.00426484136013094, 0.000610760059768012, 0.000575522836335551,
0.000823785810599276, 0.00181936036509178, -0.00073316272551871,
-0.00184238143420679, -0.000519146793923397, -0.00120324664043103,
-0.000882469414168696, -0.00148118339830283, 0.000929612782487155,
0.000565364610238817, 0.000578158613453894, 0.00060479145432879,
-0.00520576206828594, 0.000708404040882016, 0.00105224485893451,
0.000636486872540587, -0.00359655507585543, 0.000769164650506582,
0.000635701125126786, 0.000570489501935612, -0.000641260260277221,
0.000735092947873994, 0.000757195823062773, 0.000556002742616357,
-0.00207489740356159, -0.000553386431560554, 0.000511326871983186,
0.000504591469525195, -0.000749886905655472, -0.0013939718643865,
-0.000513742626250036, -0.00105021597423516, -0.00156667292147716,
0.000864563166150134, 0.00433724128055069, 0.00053855648931922,
-0.00150732363190365, 0.00052621785349416, 0.000987781100809215,
0.000560725818171903, 0.00176012436713435, -0.000594895431092368,
-0.000686229580335151, 0.00138682284509528, -0.000531964338888358,
-0.00179959148771403, 0.000574543871314503, -0.000686996216439084,
-0.000559043343629995, 0.00055881173674166, -0.000636332688477736,
-0.000623778186703561, -0.00173834148094443, -0.000567224129968125,
-0.00122578683434504, 0.00130960156515414, -0.000548203197176633,
-0.000522749285863711, -0.000820371086264871, 0.000756014225812507,
-0.000714081490558627, -0.000617600335221624, 0.000523639760748651,
-0.000578502663833191, 0.00107478825239227, 0.000612725356974764,
-0.00065509337422931, 0.000505887803587513, -0.000566716376848575,
0.000511727090058756, 0.000572807738912218, -0.000756026937699161,
0.000547948751494332, 0.000628323894238392, -0.000541350489317693,
-0.00133529454372372, -0.000590618859845904, -0.000700581963648972,
0.000735987224462775, 0.000528958898682319, 0.000838250041022448,
-0.000519084424130511, -0.00052258402856431, -0.000538130765869838,
-0.000631819887885854, 0.00054800880764283, 0.00266115500510899,
-0.000839092093771754, 0.000559253571783103, -0.000801028189803432,
-0.000608879021022801, -0.000538018076854385, -0.000689859734395171,
0.00329650346269972, 0.000765494493951024, -0.000689450477848297,
-0.000560199139975737, 0.00159082699266122, -0.00208548663121455,
-0.000598493596793759, 0.000563544422691464, 0.000626996183768824,
-0.000653166846808162, -0.000851350174739807, -0.00140687473245116,
-0.000887003220306326, -0.000765614651347946, -0.00100676206277761,
0.000724714394852555, 0.00108872127644233, -0.000678558537305918,
-0.000705087556212902, 0.000544828152248655, -0.000791700964308362,
0.000606125736727137, -0.00119335967326073, 0.00075413211796338,
0.000526038939010931, 0.00086543737231537, -0.000817788712950573,
-0.000584070926663571, 0.000619657281937691, 0.000680783312420274,
-0.000513831718574664, -0.00050972403875349, -0.00114542220685365,
-0.00070564389723593, -0.01057964950882, -0.000610357922434801,
0.000818264221596365, 0.000940825400308043, -0.000726555639413817,
-0.000591089505560305, 0.000564738888193972, -0.00068515060569041,
0.000668920238348747, -0.00110103375121717, -0.0015480433031172,
0.000663030855223568, 0.000500097431997304, -0.000600730311271391,
-0.000672397772962796, -0.000607852365856587, 0.000536711920570809,
0.000595055206488837, 0.000523123873687581, 0.000977280737528119,
0.000616410821629998, 0.000788593666889881, -0.000671642905915704,
0.000717328711735021, -0.000551853104219902, -0.000565153434708421,
-0.000802585212152707, 0.000536342062561701, 0.000682048510343591,
-0.000541902545439399, 0.000779676683974273, 0.000698841439971787,
0.000559313965908359, -0.00064986819016255, 0.000795421518319017,
0.00364973919549527, 0.000669658692276087, 0.00109045476974678,
0.000514411572742901, 0.000503832507211754, -0.000507376233564116,
0.001232871590787, 0.000561820312542594, -0.000501190337518054,
-0.000769036505996468, -0.000695537959007453, -0.000572065848166048,
-0.00167929926328192, 0.000597078186826749, 0.00710238430870014,
0.000745192112519888, -0.00116091022028009, -0.000791139281769659,
-0.00148898466632552, 0.000565144038962018, -0.000514019821833855,
-0.00148427996685285, -0.000822717245339888, -0.00062922111212238,
-0.000636011367371125, 0.00119640327632808, 0.000548455410294579,
0.000652678152560426, 0.000509244387833618, 0.000961872348987924,
0.000662064072514568, -0.00068116858054168, -0.000569930302445343,
0.00188358126928101, 0.00130560555273895, 0.000593470885775105,
0.00160093110088155, 0.000785262438315115, -0.000912313442922752,
0.000609996052359563, 0.000720137994393966, 0.000568163899000496,
0.00128685533068307, -0.000756787473447318, 0.000765932134255465,
0.00064884753100003, 0.000687571386270847, -0.000582094290400903,
-0.000693177295971736, -0.000601776208094762, 0.000503616387996786,
-0.000615095866544735, -0.000799593899689199, 0.000773750859128342,
-0.000522576090260074, 0.000503578107212022, -0.00104492224837571,
0.000547928732299141, 0.00310304337507183, 0.000893382870797765,
-0.000577792878910799, -0.000647710366578735, -0.00061992948706191,
0.000825702487162516, 0.000606579510524341, 0.000552792484727505,
0.000688600840895504, 0.000505093563534231, -0.000728420573667066,
-0.00157924525963438, -0.000603846616019865, -0.000521941317177976,
0.00150498158245682, -0.000584572670337735, 0.000713757870583365,
0.000524287801789924, 0.00107217649464886, 0.00213147531822244,
0.000566012832157625, -0.00069828890607937, 0.000641567963736378,
-0.000509531713644762, -0.000547564140049417, -0.00115275240244728,
0.000560465768010943, -0.000651807371497171, -0.00096487058986483,
0.000753687665266511, -0.000665599418910645, -0.000691278087025182,
-0.000578010050725553, -0.000685833148198256, 0.000698470819832764,
0.00102943368139208, -0.000725840586788706, 0.00125882415960632,
-0.000630791474954151, -0.000764813558678412, -0.000638539347184164,
0.000654486496518558, 0.000547453642294471, 0.000572020020495501,
-0.000605791001705214, 0.00660211658324172, 0.00114928683282756,
0.000985676480677711, -0.000694668292547718, -0.000528955637964401,
0.000647975568638159, 0.00116454536417443, 0.000506748841724303,
-0.000500925156604382, -0.000567015088082101, 0.00128711230206946,
0.000533633762033858, 0.00505991432758357, 0.000518058378462527,
-0.000592822519784875, 0.00177414999018666, 0.00059845426944527,
-0.000511614433724716, 0.0016614697907098, 0.000852196464322219,
0.00241689725305427, -0.000614317948913978, -0.000729717143318709,
-0.000612900648802039, -0.000727983564232204, -0.000694965869158182,
-0.000527752006066251, -0.000584233784708843, 0.000522097476268968,
0.000543092880677776, 0.000947121210698398, -0.00241810275096377,
0.00181893137435019, 0.000931873879297385, 0.000512116215015013,
0.000724985702444059, -0.000566713495050664, 0.000603953591362227
)
After fitting the data look like the following (empirical histogram-blue, theoretical histogram-orange):
The same when include freq=FALSE in hist
This will all be due to anomalous values in a not represented by the histogram you've shown. This could be the cause of both the very high kurtotsis, and the vgFit() algorithm failing to find a good fit.
Type dput(a[abs(a) > 5e-4]) in the console and copy the output into your question. People then may be able to recreate aomething like the vector a without having to get all 16000 values and debug the vgFit issue.
Thanks for the extra data. There are some extreme values in there, but I don;t think those are what is causing the problem in vgFit. Fitting 4 parameters which can be almost any value is difficult, but you can help it along by rescaling your data to something typical. Try this:
b <- (a-mean(a))/sd(a)
vgf <- vgFit(b)
vgf$param
VG <- rvg(16000, param = vgf$param)
VG_rescaled <- VG*sd(a)+mean(a)
hist(VG_rescaled, breaks=140, col="orange", main="", xlab="Value")
and see if the two histograms are close enough now.
I am using ggplot command to graph my data. I get this picture:
The graph is ugly and discontinus; I look at another data and it looks like this:
(Their code is in Step 2 on link: https://www.datascience.com/blog/introduction-to-forecasting-with-arima-in-r-learn-data-science-tutorials). How do I make my data look like that?
My code is:
daily_data = read.csv2('satni3.csv', header=TRUE, stringsAsFactors=FALSE)
daily_data$Date = as.Date(daily_data$datum, format = "%Y.%m.%d")
ggplot(daily_data, aes(Date, total)) +
geom_line() +
scale_x_date('day') +
ylab("opterecenje") +
xlab("") # total is the column with numbers that represent some data for a given date.
Here is half of my data (there is 15 days missing), pasted, first are names of columns and then there is data:
datum;sat;a11;a21;a31;b11;b21;b31;b41;b51;c11;c21;xx1;a12;a22;a32;b12;b22;b32;b42;b52;c12;c22;xx2;sum1;sum2;total
01.01.2017;00:00;0;16;0;0;0;0;0;0;0;0;0;0;62;4;1;0;0;0;0;1;0;0;16;68;84
01.01.2017;01:00;0;76;3;0;0;0;0;0;0;0;0;0;3;0;0;0;0;0;0;0;0;0;79;3;82
01.01.2017;02:00;0;151;3;2;0;0;0;0;0;0;0;0;15;0;1;0;0;0;0;0;0;0;156;16;172
01.01.2017;03:00;0;72;1;0;0;0;0;0;0;0;0;0;19;1;0;0;0;0;0;0;0;0;73;20;93
01.01.2017;04:00;0;32;1;0;0;0;0;0;0;0;0;0;18;0;0;0;0;0;0;0;0;0;33;18;51
01.01.2017;05:00;0;22;3;0;0;0;0;0;0;0;0;0;15;0;0;0;0;0;0;0;0;0;25;15;40
01.01.2017;06:00;0;7;0;0;0;0;0;0;0;0;0;0;12;0;0;0;0;0;0;1;0;0;7;13;20
01.01.2017;07:00;0;10;0;0;0;0;0;0;0;0;0;0;14;0;0;0;0;0;0;1;0;0;10;15;25
01.01.2017;08:00;0;9;0;0;0;0;0;0;0;0;0;0;11;0;0;0;0;0;0;0;0;0;9;11;20
01.01.2017;09:00;0;12;0;0;0;0;0;0;0;0;0;0;15;0;1;0;0;0;0;1;0;0;12;17;29
01.01.2017;10:00;0;15;2;0;0;0;0;0;1;0;0;0;21;1;1;0;0;0;0;0;0;0;18;23;41
01.01.2017;11:00;0;37;1;0;0;0;0;0;0;0;0;0;28;3;0;0;0;0;0;0;0;0;38;31;69
01.01.2017;12:00;0;59;2;0;0;0;0;0;0;0;0;0;78;2;0;0;0;0;0;1;0;0;61;81;142
01.01.2017;13:00;0;70;0;1;0;0;0;0;2;0;0;3;88;4;0;0;0;0;0;1;0;1;73;97;170
01.01.2017;14:00;1;123;5;1;0;0;0;0;0;0;0;0;110;3;0;0;0;0;1;0;0;0;130;114;244
01.01.2017;15:00;3;101;0;0;0;0;0;0;6;0;0;1;124;4;0;0;0;0;0;1;0;0;110;130;240
01.01.2017;16:00;0;128;1;0;0;0;0;0;2;0;0;0;117;3;0;0;0;0;0;0;0;0;131;120;251
01.01.2017;17:00;1;126;2;0;0;0;0;1;0;0;0;0;99;2;0;0;0;0;0;1;0;0;130;102;232
01.01.2017;18:00;0;91;1;0;0;0;0;0;1;0;0;0;84;1;1;0;0;0;0;0;0;0;93;86;179
01.01.2017;19:00;0;69;2;1;0;0;0;0;1;0;0;0;53;4;1;0;0;0;0;0;0;0;73;58;131
01.01.2017;20:00;0;54;0;0;0;0;0;0;0;0;0;0;42;1;0;0;0;0;0;2;0;0;54;45;99
01.01.2017;21:00;0;36;0;0;1;0;0;0;1;0;0;0;22;1;0;0;0;0;0;2;0;0;38;25;63
01.01.2017;22:00;0;11;0;1;0;0;0;0;1;0;0;0;19;2;1;0;0;0;0;0;0;0;13;22;35
01.01.2017;23:00;0;15;0;0;0;0;0;0;2;0;0;0;12;2;0;0;0;0;0;0;0;0;17;14;31
02.01.2017;00:00;0;13;0;0;0;0;0;0;0;0;0;0;10;0;0;0;0;0;0;1;0;0;13;11;24
02.01.2017;01:00;0;11;1;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;12;2;14
02.01.2017;02:00;0;2;0;0;0;0;0;0;0;0;0;0;5;0;0;0;0;0;1;0;0;0;2;6;8
02.01.2017;03:00;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0
02.01.2017;04:00;0;5;0;0;0;0;0;1;1;0;0;0;1;1;0;0;0;0;0;0;0;0;7;2;9
02.01.2017;05:00;0;3;2;0;0;0;0;1;2;0;0;0;3;1;0;0;0;0;0;2;0;0;8;6;14
02.01.2017;06:00;0;8;1;0;0;0;0;0;0;0;0;0;7;1;0;0;0;0;0;2;0;1;9;11;20
02.01.2017;07:00;0;29;0;0;0;0;0;1;1;0;0;1;37;1;0;0;1;0;0;3;0;0;31;43;74
02.01.2017;08:00;0;60;1;0;0;1;0;0;0;0;0;0;75;5;3;1;2;0;0;1;0;0;62;87;149
02.01.2017;09:00;0;54;1;2;2;4;0;0;1;0;0;0;62;3;5;1;2;1;1;0;0;0;64;75;139
02.01.2017;10:00;0;52;7;1;1;9;1;0;1;0;0;0;63;4;2;1;7;0;0;1;0;0;72;78;150
02.01.2017;11:00;0;76;5;5;0;6;0;1;2;0;0;0;81;5;2;2;7;0;1;0;0;0;95;98;193
02.01.2017;12:00;0;86;6;6;0;10;0;1;0;0;0;0;79;9;1;1;9;0;3;2;0;0;109;104;213
02.01.2017;13:00;0;92;2;2;2;4;0;0;2;0;0;0;69;3;0;3;7;0;1;0;0;0;104;83;187
02.01.2017;14:00;0;104;2;3;2;9;0;1;0;0;0;0;99;7;3;4;5;1;1;1;0;1;121;122;243
02.01.2017;15:00;0;117;8;0;1;11;0;0;1;0;0;0;85;2;0;1;11;0;0;2;0;0;138;101;239
02.01.2017;16:00;0;113;4;0;0;6;0;0;4;0;0;0;95;1;2;0;8;0;0;3;0;0;127;109;236
02.01.2017;17:00;0;115;3;2;0;3;0;0;2;0;0;0;114;3;0;0;5;0;0;1;0;0;125;123;248
02.01.2017;18:00;0;61;1;1;1;1;0;1;0;0;0;0;94;4;1;0;0;0;0;2;0;0;66;101;167
02.01.2017;19:00;0;72;4;0;0;0;0;0;1;0;0;0;57;3;0;0;0;0;0;1;0;0;77;61;138
02.01.2017;20:00;0;33;1;0;0;0;0;0;2;0;0;0;45;1;0;0;0;0;0;1;0;0;36;47;83
02.01.2017;21:00;0;21;1;1;0;0;0;1;0;0;0;0;32;4;0;0;0;0;0;4;0;0;24;40;64
02.01.2017;22:00;0;22;1;0;0;0;0;0;0;0;0;0;24;1;0;0;0;0;0;1;0;0;23;26;49
02.01.2017;23:00;0;23;0;0;0;0;0;1;1;0;0;0;20;0;0;0;0;0;0;0;0;0;25;20;45
03.01.2017;00:00;0;12;0;0;0;0;0;0;1;0;0;0;8;1;0;0;0;0;0;0;0;0;13;9;22
03.01.2017;01:00;0;7;0;0;0;0;0;0;0;0;0;0;3;0;0;0;0;0;0;1;0;0;7;4;11
03.01.2017;02:00;0;4;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;4;1;5
03.01.2017;03:00;0;2;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;2;1;3
03.01.2017;04:00;0;2;0;0;0;0;0;0;0;0;0;0;1;1;0;0;0;0;0;0;0;0;2;2;4
03.01.2017;05:00;0;3;1;0;0;0;0;0;3;0;0;0;3;0;0;0;0;0;0;1;0;0;7;4;11
03.01.2017;06:00;0;7;0;0;0;0;0;0;1;0;0;0;3;2;0;0;0;0;0;2;0;0;8;7;15
03.01.2017;07:00;0;25;2;0;0;0;0;0;2;0;0;0;51;5;0;1;0;0;1;2;0;0;29;60;89
03.01.2017;08:00;0;64;4;1;1;5;0;0;0;0;0;1;86;4;4;1;2;1;0;2;0;0;75;101;176
03.01.2017;09:00;0;65;4;2;3;6;0;0;1;0;0;0;57;8;5;0;7;0;0;1;0;1;81;79;160
03.01.2017;10:00;0;84;13;3;3;8;1;2;1;0;0;0;75;9;1;0;5;0;0;0;0;1;115;91;206
03.01.2017;11:00;0;74;3;3;1;6;0;2;1;0;0;0;75;7;1;3;4;0;0;0;0;0;90;90;180
03.01.2017;12:00;1;96;10;3;3;9;0;0;0;0;1;1;85;6;5;1;7;0;0;0;0;0;123;105;228
03.01.2017;13:00;1;115;9;5;0;1;0;0;1;0;1;0;103;7;2;2;6;0;0;1;1;1;133;123;256
03.01.2017;14:00;2;120;8;3;0;6;0;0;1;0;3;0;114;6;2;1;5;0;0;1;0;0;143;129;272
03.01.2017;15:00;0;116;6;1;0;7;0;0;1;0;0;1;107;5;1;0;6;0;0;3;0;0;131;123;254
03.01.2017;16:00;0;139;4;0;0;5;0;0;4;0;0;0;115;5;1;0;4;0;0;2;0;0;152;127;279
03.01.2017;17:00;0;117;1;2;0;4;0;0;1;0;0;0;126;5;0;2;5;0;0;1;0;0;125;139;264
03.01.2017;18:00;0;86;4;1;0;1;0;0;0;0;0;0;86;4;0;0;0;0;0;2;0;0;92;92;184
03.01.2017;19:00;0;57;2;0;0;0;0;0;1;0;0;0;70;3;1;0;0;0;1;0;0;1;60;76;136
03.01.2017;20:00;0;40;4;0;0;0;0;1;1;0;0;0;48;3;2;1;0;0;0;0;0;0;46;54;100
03.01.2017;21:00;0;37;0;1;0;0;0;0;1;0;0;0;37;3;0;0;0;0;0;3;0;0;39;43;82
03.01.2017;22:00;0;25;1;1;0;0;0;0;0;0;0;0;20;0;0;0;0;0;0;1;0;0;27;21;48
03.01.2017;23:00;0;23;0;0;0;1;0;0;2;0;0;0;14;0;0;0;0;0;0;0;0;0;26;14;40
04.01.2017;00:00;0;20;0;0;0;0;0;0;0;0;0;0;10;0;1;0;0;0;0;0;0;0;20;11;31
04.01.2017;01:00;0;4;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;4;1;5
04.01.2017;02:00;0;4;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;4;0;4
04.01.2017;03:00;0;1;0;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;1;2;3
04.01.2017;04:00;0;1;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;2;0;2
04.01.2017;05:00;0;4;1;0;0;0;0;0;2;0;0;0;5;1;0;0;0;0;0;1;0;0;7;7;14
04.01.2017;06:00;0;10;0;0;0;0;0;0;0;0;0;0;5;2;0;0;0;0;0;1;0;0;10;8;18
04.01.2017;07:00;0;31;3;0;0;0;0;0;2;0;0;1;43;6;0;1;1;0;0;3;0;0;36;55;91
04.01.2017;08:00;0;74;4;1;0;4;0;0;0;0;1;0;91;2;6;1;3;0;0;2;0;0;84;105;189
04.01.2017;09:00;0;65;10;2;2;7;0;0;1;0;0;0;73;8;3;1;8;0;0;1;0;0;87;94;181
04.01.2017;10:00;0;76;4;4;2;9;0;0;1;0;0;0;66;9;2;1;8;0;0;0;0;0;96;86;182
04.01.2017;11:00;1;94;1;3;2;7;0;1;1;0;0;0;89;4;2;2;8;0;0;0;0;0;110;105;215
04.01.2017;12:00;0;96;4;1;0;8;0;0;0;0;0;0;92;5;2;2;7;0;0;0;0;1;109;109;218
04.01.2017;13:00;1;98;4;3;1;3;0;0;1;0;0;1;100;10;4;5;5;0;0;2;0;0;111;127;238
04.01.2017;14:00;1;120;3;0;0;8;0;1;1;0;0;0;84;4;3;0;5;0;0;2;0;1;134;99;233
04.01.2017;15:00;0;114;3;1;1;4;0;0;1;0;0;0;109;5;0;1;7;0;1;2;0;1;124;126;250
04.01.2017;16:00;1;139;6;1;0;6;0;0;4;0;0;0;109;4;4;1;7;0;0;2;0;0;157;127;284
04.01.2017;17:00;1;115;5;3;0;7;0;0;1;0;0;1;124;5;0;0;8;0;0;1;0;0;132;139;271
04.01.2017;18:00;0;77;2;1;0;0;0;0;1;0;0;0;97;3;1;0;0;0;0;2;0;0;81;103;184
04.01.2017;19:00;0;51;4;0;0;0;0;0;2;0;0;0;57;1;0;0;0;0;0;0;0;0;57;58;115
04.01.2017;20:00;0;41;3;1;0;0;0;0;1;0;0;0;34;2;0;0;0;0;1;2;0;0;46;39;85
04.01.2017;21:00;0;24;2;0;0;0;0;0;1;0;0;0;38;0;0;0;0;0;0;4;0;0;27;42;69
04.01.2017;22:00;0;23;1;0;0;0;0;0;1;0;0;0;30;2;0;0;0;0;0;0;0;0;25;32;57
04.01.2017;23:00;0;16;0;0;0;0;0;0;2;0;0;0;14;0;0;0;0;0;0;0;0;0;18;14;32
05.01.2017;00:00;0;19;0;0;0;0;0;0;0;0;0;0;10;1;1;0;0;0;0;1;0;0;19;13;32
05.01.2017;01:00;0;9;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;9;1;10
05.01.2017;02:00;0;2;0;0;0;0;0;0;0;0;0;0;3;0;0;0;0;0;0;0;0;0;2;3;5
05.01.2017;03:00;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;1
05.01.2017;04:00;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;1
05.01.2017;05:00;0;4;1;0;0;1;0;0;2;0;0;0;4;1;0;0;0;0;0;1;0;0;8;6;14
05.01.2017;06:00;0;7;0;0;0;0;0;0;1;0;0;0;5;1;0;0;0;0;0;2;0;0;8;8;16
05.01.2017;07:00;0;30;1;0;0;0;0;0;2;0;0;1;49;5;0;0;0;0;0;4;0;0;33;59;92
05.01.2017;08:00;0;57;1;3;2;8;0;0;0;0;0;0;76;4;3;3;9;0;0;2;0;0;71;97;168
05.01.2017;09:00;0;64;4;4;1;12;0;0;1;0;0;0;59;6;2;1;9;0;0;1;0;0;86;78;164
05.01.2017;10:00;2;85;8;2;2;10;0;0;1;0;0;1;40;8;3;2;11;0;0;0;0;0;110;65;175
05.01.2017;11:00;0;86;9;6;2;8;0;0;2;0;0;0;91;6;5;3;11;0;0;0;0;0;113;116;229
05.01.2017;12:00;0;127;9;2;6;10;0;0;0;0;0;1;93;7;1;3;11;0;0;2;0;1;154;119;273
05.01.2017;13:00;0;105;5;1;2;4;0;0;3;0;0;0;107;9;3;3;5;0;0;0;1;1;120;129;249
05.01.2017;14:00;0;104;5;3;0;14;0;0;1;0;0;0;117;7;2;0;10;0;0;1;0;1;127;138;265
05.01.2017;15:00;0;117;5;1;0;10;0;0;1;0;0;0;107;6;3;0;10;0;1;3;0;1;134;131;265
05.01.2017;16:00;0;126;4;1;0;11;0;0;4;0;0;0;110;3;3;2;10;0;1;1;0;0;146;130;276
05.01.2017;17:00;0;132;4;0;1;10;0;0;1;0;0;0;118;1;2;0;13;0;0;0;0;0;148;134;282
05.01.2017;18:00;0;71;2;2;0;1;0;0;1;0;0;0;89;3;0;1;1;0;0;1;0;0;77;95;172
05.01.2017;19:00;0;57;5;0;0;0;0;0;2;0;0;0;57;2;0;0;0;0;0;1;0;0;64;60;124
05.01.2017;20:00;0;43;4;1;0;0;0;0;1;0;0;0;58;3;0;0;0;0;1;2;0;0;49;64;113
05.01.2017;21:00;0;36;0;0;0;0;0;1;1;0;0;0;50;1;0;0;0;0;0;4;0;0;38;55;93
05.01.2017;22:00;0;30;2;0;0;0;0;0;1;0;0;0;25;1;0;0;0;0;0;1;0;0;33;27;60
05.01.2017;23:00;0;17;1;0;0;0;0;0;2;0;0;0;22;0;1;0;0;0;0;0;0;0;20;23;43
06.01.2017;00:00;0;22;0;0;0;0;0;0;0;0;0;0;13;1;0;0;0;0;0;1;0;0;22;15;37
06.01.2017;01:00;0;5;0;0;0;0;0;0;0;0;0;0;6;1;0;0;0;0;0;0;0;0;5;7;12
06.01.2017;02:00;0;4;0;0;0;0;0;0;0;0;0;0;5;0;0;0;0;0;0;0;0;0;4;5;9
06.01.2017;03:00;0;1;1;0;0;0;0;0;0;0;0;0;5;0;0;0;0;0;0;0;0;0;2;5;7
06.01.2017;04:00;0;2;0;0;0;0;0;0;0;0;0;0;1;1;0;0;0;0;0;0;0;0;2;2;4
06.01.2017;05:00;0;5;1;0;0;0;0;0;2;0;0;0;7;0;0;0;0;0;0;0;0;0;8;7;15
06.01.2017;06:00;0;6;0;0;0;0;0;0;1;0;0;0;7;3;0;0;0;0;0;0;0;0;7;10;17
06.01.2017;07:00;0;7;1;0;0;0;0;0;0;0;0;0;21;4;0;1;0;0;0;1;0;0;8;27;35
06.01.2017;08:00;0;10;3;0;0;0;0;0;0;0;0;0;27;1;1;0;0;0;0;1;0;0;13;30;43
06.01.2017;09:00;0;15;3;1;0;1;0;0;1;0;0;0;24;1;2;0;0;0;0;1;0;0;21;28;49
06.01.2017;10:00;0;41;3;0;3;2;0;0;0;0;0;0;39;1;2;1;2;0;0;0;0;0;49;45;94
06.01.2017;11:00;0;61;3;1;1;1;0;0;1;0;0;0;52;1;0;0;1;0;0;0;0;0;68;54;122
06.01.2017;12:00;0;61;5;1;0;0;0;0;0;0;0;0;58;1;1;0;1;0;0;2;0;1;67;64;131
06.01.2017;13:00;0;73;4;0;0;0;0;0;1;0;0;0;59;1;1;0;1;0;0;1;0;0;78;63;141
06.01.2017;14:00;0;67;2;2;0;0;0;0;1;0;0;0;91;3;1;0;0;0;0;2;0;0;72;97;169
06.01.2017;15:00;0;93;4;1;0;0;0;0;0;0;0;0;78;2;2;0;0;0;0;0;0;0;98;82;180
06.01.2017;16:00;0;67;3;0;0;0;0;0;1;0;0;0;72;1;0;1;0;0;0;1;0;0;71;75;146
06.01.2017;17:00;0;73;1;1;0;0;0;0;1;0;0;0;59;2;1;0;0;0;0;0;0;0;76;62;138
06.01.2017;18:00;0;59;1;0;0;0;0;0;0;0;0;0;61;3;0;0;0;0;0;0;0;0;60;64;124
06.01.2017;19:00;0;47;2;0;0;0;0;0;0;0;0;0;42;0;0;0;0;0;0;1;0;0;49;43;92
06.01.2017;20:00;0;31;0;0;0;0;0;0;1;0;0;0;28;1;0;0;0;0;1;1;0;0;32;31;63
06.01.2017;21:00;0;51;0;0;0;0;0;0;1;0;0;0;26;0;0;0;0;0;0;2;0;0;52;28;80
06.01.2017;22:00;0;16;0;0;0;0;0;0;1;0;0;0;19;0;0;0;0;0;0;0;0;0;17;19;36
06.01.2017;23:00;0;24;0;0;0;0;0;0;1;1;0;0;31;0;1;0;0;0;0;0;0;0;26;32;58
07.01.2017;00:00;0;12;1;0;0;0;0;0;0;0;0;0;9;2;0;0;0;0;0;0;0;0;13;11;24
07.01.2017;01:00;0;9;0;0;0;0;0;0;0;0;0;0;5;1;0;0;0;0;0;1;0;0;9;7;16
07.01.2017;02:00;0;2;1;0;0;0;0;0;0;0;0;0;8;0;0;0;0;0;0;0;0;0;3;8;11
07.01.2017;03:00;0;2;0;0;0;0;0;0;0;0;0;0;4;0;0;0;0;0;0;0;0;0;2;4;6
07.01.2017;04:00;0;1;1;0;0;0;0;0;0;0;0;0;4;0;0;0;0;0;0;0;0;0;2;4;6
07.01.2017;05:00;0;2;0;0;0;0;0;0;2;0;0;0;7;1;0;0;0;0;0;0;0;0;4;8;12
07.01.2017;06:00;0;4;0;0;0;0;0;0;1;0;0;0;5;1;0;0;0;0;0;2;0;0;5;8;13
07.01.2017;07:00;0;7;3;0;0;0;0;0;0;0;0;0;28;5;1;0;0;0;0;2;0;0;10;36;46
07.01.2017;08:00;0;24;2;0;1;8;0;0;0;0;1;0;31;2;2;0;7;0;1;2;0;0;36;45;81
07.01.2017;09:00;0;32;5;1;3;11;0;0;1;0;0;0;30;2;3;1;16;0;0;1;0;0;53;53;106
07.01.2017;10:00;0;47;4;3;3;14;0;0;0;0;0;0;41;3;2;3;12;0;0;1;0;1;71;63;134
07.01.2017;11:00;0;83;5;3;2;12;0;0;1;0;0;0;54;2;2;1;9;0;0;0;0;0;106;68;174
07.01.2017;12:00;0;90;2;1;1;7;0;0;0;0;0;0;78;2;1;1;10;0;0;1;0;2;101;95;196
07.01.2017;13:00;2;84;3;2;2;6;0;0;1;0;0;0;81;5;3;0;7;0;0;1;0;1;100;98;198
07.01.2017;14:00;0;78;6;2;2;8;1;0;1;0;0;0;86;3;1;1;10;0;0;2;0;1;98;104;202
07.01.2017;15:00;0;92;3;2;4;8;0;0;0;0;0;0;81;3;1;1;10;0;0;1;0;2;109;99;208
07.01.2017;16:00;0;83;2;1;0;0;0;1;2;0;0;0;92;2;1;0;0;0;1;0;0;0;89;96;185
07.01.2017;17:00;0;69;2;1;1;3;0;0;1;0;0;0;83;0;1;0;0;0;0;0;0;0;77;84;161
07.01.2017;18:00;0;59;2;0;0;0;0;0;0;0;0;0;70;1;1;0;1;0;0;0;0;0;61;73;134
07.01.2017;19:00;0;38;2;1;0;1;0;0;2;0;0;0;42;1;3;0;1;0;0;1;0;0;44;48;92
07.01.2017;20:00;0;41;0;0;0;0;0;0;0;0;0;0;51;2;0;0;0;0;0;0;0;0;41;53;94
07.01.2017;21:00;0;22;1;0;0;0;0;0;2;0;0;0;31;2;1;0;0;0;0;1;0;0;25;35;60
07.01.2017;22:00;0;26;0;0;0;0;0;0;0;0;0;0;18;0;0;0;0;0;0;1;0;0;26;19;45
07.01.2017;23:00;0;26;0;0;0;0;0;0;2;0;0;0;16;0;1;0;0;0;0;0;0;0;28;17;45
08.01.2017;00:00;0;19;0;0;0;0;0;0;0;0;0;0;13;1;1;0;0;0;0;1;0;0;19;16;35
08.01.2017;01:00;0;9;0;0;0;0;0;0;0;0;0;0;8;0;0;0;0;0;0;0;0;0;9;8;17
08.01.2017;02:00;0;4;0;0;0;0;0;1;0;0;0;0;4;0;0;0;0;0;0;0;0;0;5;4;9
08.01.2017;03:00;0;1;0;0;0;0;0;0;0;0;0;0;4;0;0;0;0;0;0;0;0;0;1;4;5
08.01.2017;04:00;0;1;0;0;0;0;0;0;1;0;0;0;1;0;0;0;0;0;0;0;0;0;2;1;3
08.01.2017;05:00;0;2;2;0;0;0;0;1;0;0;0;0;5;1;0;0;0;0;0;1;0;0;5;7;12
08.01.2017;06:00;0;0;0;0;0;0;0;0;0;0;0;0;6;1;0;0;0;0;0;1;0;0;0;8;8
08.01.2017;07:00;0;9;2;1;0;0;0;0;0;0;0;0;15;5;0;0;0;0;1;2;0;0;12;23;35
08.01.2017;08:00;0;16;1;0;0;5;0;0;0;0;0;0;16;0;1;0;3;0;0;0;0;0;22;20;42
08.01.2017;09:00;0;18;1;2;1;4;0;0;1;0;0;0;19;1;0;0;4;0;0;1;0;0;27;25;52
08.01.2017;10:00;0;49;1;1;1;5;0;1;1;0;0;0;30;1;0;0;5;0;1;0;0;0;59;37;96
08.01.2017;11:00;0;47;0;1;0;4;0;2;0;0;0;0;35;1;2;1;5;0;1;0;0;0;54;45;99
08.01.2017;12:00;0;79;2;1;0;7;0;0;0;0;0;0;52;2;0;0;6;0;0;1;0;0;89;61;150
08.01.2017;13:00;0;67;3;2;0;4;0;0;1;0;0;0;63;1;1;0;6;0;0;1;0;0;77;72;149
08.01.2017;14:00;1;74;3;2;0;0;0;1;0;0;0;0;85;1;0;0;0;0;0;2;0;0;81;88;169
08.01.2017;15:00;2;74;5;1;1;0;0;0;0;0;0;0;61;8;1;0;0;0;0;1;0;0;83;71;154
08.01.2017;16:00;0;65;0;0;0;0;0;0;4;0;0;0;81;1;0;0;0;0;0;0;0;0;69;82;151
08.01.2017;17:00;0;71;2;0;0;0;0;0;1;0;0;0;62;2;0;0;1;0;0;0;0;0;74;65;139
08.01.2017;18:00;0;55;0;0;0;0;0;0;0;0;0;0;61;0;0;2;0;0;0;2;0;0;55;65;120
08.01.2017;19:00;0;33;1;0;0;0;0;0;0;0;0;0;36;0;0;0;0;0;0;1;0;0;34;37;71
08.01.2017;20:00;0;42;0;0;0;0;0;0;1;0;0;0;38;1;0;0;0;0;0;1;0;0;43;40;83
08.01.2017;21:00;0;23;0;0;0;0;0;0;2;0;0;1;28;1;0;0;0;0;0;3;0;0;25;33;58
08.01.2017;22:00;0;22;0;0;0;0;0;0;0;0;0;0;20;0;0;0;0;0;0;2;0;0;22;22;44
08.01.2017;23:00;0;9;1;0;0;0;0;0;1;1;0;0;9;2;0;0;0;0;0;0;0;0;12;11;23
09.01.2017;00:00;0;4;0;0;0;0;0;0;0;0;0;0;5;0;0;0;0;0;0;0;0;0;4;5;9
09.01.2017;01:00;0;1;0;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;1;2;3
09.01.2017;02:00;0;5;0;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;5;2;7
09.01.2017;03:00;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;1
09.01.2017;04:00;0;3;0;0;0;0;0;0;1;0;0;0;3;1;0;0;0;0;0;0;0;0;4;4;8
09.01.2017;05:00;0;2;1;0;0;0;0;0;3;0;0;0;2;0;0;0;0;0;0;1;0;0;6;3;9
09.01.2017;06:00;0;6;0;0;0;0;0;0;0;0;0;0;3;2;0;0;1;0;1;1;0;0;6;8;14
09.01.2017;07:00;0;42;0;1;0;1;0;1;2;0;0;0;50;6;1;1;0;0;0;3;0;0;47;61;108
09.01.2017;08:00;0;84;2;5;2;5;0;0;0;0;0;0;111;6;5;1;6;0;0;2;0;0;98;131;229
09.01.2017;09:00;1;66;5;0;0;9;0;0;1;0;0;0;74;4;8;2;6;0;0;0;0;0;82;94;176
09.01.2017;10:00;0;80;8;2;4;8;0;0;1;0;0;0;64;5;2;3;11;0;0;0;0;0;103;85;188
09.01.2017;11:00;1;80;3;8;2;6;0;2;2;0;0;0;76;7;2;4;7;1;1;0;0;0;104;98;202
09.01.2017;12:00;0;107;7;5;3;9;0;1;0;0;0;0;80;9;3;2;11;0;1;1;0;0;132;107;239
09.01.2017;13:00;0;103;8;3;1;5;0;1;1;0;1;0;76;6;3;2;8;0;0;1;0;0;123;96;219
09.01.2017;14:00;0;106;5;2;2;5;0;0;1;0;0;0;92;8;4;1;9;1;0;2;0;0;121;117;238
09.01.2017;15:00;0;97;8;4;1;8;0;0;1;0;0;0;104;2;5;4;7;1;0;2;0;0;119;125;244
09.01.2017;16:00;0;112;6;5;4;8;1;0;3;0;0;0;106;4;4;3;9;0;1;1;0;0;139;128;267
09.01.2017;17:00;0;109;6;3;3;7;0;0;2;0;0;0;129;6;1;0;8;0;0;1;0;1;130;146;276
09.01.2017;18:00;0;86;5;1;0;0;1;0;0;0;0;0;92;4;1;0;0;0;1;2;0;0;93;100;193
09.01.2017;19:00;0;47;2;0;0;0;0;0;2;0;0;0;60;3;0;0;0;0;0;1;0;0;51;64;115
09.01.2017;20:00;0;42;2;0;0;0;0;0;1;0;0;0;45;2;0;0;0;0;0;1;0;0;45;48;93
09.01.2017;21:00;0;38;2;0;0;1;0;1;0;0;0;0;34;0;0;0;0;0;0;2;0;0;42;36;78
09.01.2017;22:00;0;29;0;0;0;1;0;0;1;0;0;0;22;0;0;0;0;0;0;1;0;0;31;23;54
09.01.2017;23:00;0;17;2;0;0;0;0;0;2;0;0;0;25;0;0;0;0;0;0;0;0;0;21;25;46
10.01.2017;00:00;0;14;0;0;0;0;0;0;0;0;0;0;6;1;0;0;0;0;0;1;0;0;14;8;22
10.01.2017;01:00;0;8;0;0;0;0;0;0;0;0;0;0;4;0;0;0;0;0;0;0;0;0;8;4;12
10.01.2017;02:00;0;2;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;2;1;3
10.01.2017;03:00;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;0;0;0;1;1;2
10.01.2017;04:00;0;1;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;2;0;2
10.01.2017;05:00;0;7;0;0;0;0;0;1;2;0;0;0;5;1;1;0;0;0;0;1;0;0;10;8;18
10.01.2017;06:00;0;6;0;1;0;0;0;0;0;0;0;0;9;1;0;0;0;0;0;2;0;0;7;12;19
10.01.2017;07:00;0;38;3;0;0;0;0;0;2;0;0;0;49;7;0;1;1;0;0;3;0;0;43;61;104
10.01.2017;08:00;0;78;2;1;2;7;0;0;0;0;0;0;102;0;3;0;9;2;0;2;0;1;90;119;209
10.01.2017;09:00;0;85;9;2;3;8;0;0;1;0;0;0;91;6;3;2;11;0;0;0;0;0;108;113;221
10.01.2017;10:00;0;62;8;2;4;10;1;0;1;0;0;0;68;9;3;1;9;0;0;0;0;0;88;90;178
10.01.2017;11:00;1;100;11;4;4;9;0;0;2;0;0;0;80;10;6;1;9;0;0;0;0;0;131;106;237
10.01.2017;12:00;0;96;5;6;5;9;0;0;0;0;0;0;101;5;1;2;14;1;0;1;0;0;121;125;246
10.01.2017;13:00;1;115;7;3;2;4;0;0;2;0;0;1;91;9;2;2;7;0;0;1;0;0;134;113;247
10.01.2017;14:00;0;93;7;2;4;6;2;0;1;0;0;1;107;5;3;3;4;0;0;1;0;3;115;127;242
10.01.2017;15:00;0;126;7;1;2;7;1;0;1;0;0;0;110;2;1;3;9;2;0;3;0;0;145;130;275
10.01.2017;16:00;0;133;4;1;2;10;0;1;3;0;0;0;114;6;4;0;7;0;0;2;0;0;154;133;287
10.01.2017;17:00;0;111;7;6;4;5;1;0;1;0;0;0;118;7;1;2;8;0;1;0;0;0;135;137;272
10.01.2017;18:00;0;78;6;1;0;2;0;0;0;0;0;0;84;2;2;0;0;0;1;2;0;0;87;91;178
10.01.2017;19:00;0;51;5;0;1;0;0;0;1;0;0;0;60;1;3;0;0;0;0;0;0;0;58;64;122
10.01.2017;20:00;0;39;1;2;1;0;0;0;1;0;0;0;54;2;2;0;0;0;1;0;0;0;44;59;103
10.01.2017;21:00;0;32;0;0;0;0;0;0;1;0;0;0;34;2;0;0;0;0;0;3;0;0;33;39;72
10.01.2017;22:00;0;22;0;0;0;0;0;0;1;0;0;0;25;0;0;0;0;0;0;1;0;0;23;26;49
10.01.2017;23:00;0;22;2;0;0;0;0;0;2;0;0;0;21;0;1;0;0;0;0;0;0;0;26;22;48
11.01.2017;00:00;0;9;0;0;0;0;0;0;0;0;0;0;5;1;0;0;0;0;0;0;0;0;9;6;15
11.01.2017;01:00;0;6;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;7;0;7
11.01.2017;02:00;0;1;0;0;0;0;0;0;0;0;0;0;2;2;0;0;0;0;0;0;0;0;1;4;5
11.01.2017;03:00;0;1;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;1;1;2
11.01.2017;04:00;0;0;1;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;2;0;2
11.01.2017;05:00;0;5;0;0;0;0;0;1;1;0;0;0;4;0;0;0;0;0;0;1;0;0;7;5;12
11.01.2017;06:00;0;8;0;0;0;0;0;0;0;0;0;0;9;1;0;0;0;0;0;2;0;0;8;12;20
11.01.2017;07:00;0;35;3;0;0;1;0;0;2;0;0;0;52;9;0;1;1;0;0;3;0;0;41;66;107
11.01.2017;08:00;0;80;2;4;2;7;0;0;0;0;0;0;99;3;4;1;9;0;1;2;0;0;95;119;214
11.01.2017;09:00;0;78;3;2;2;8;0;0;1;0;0;0;84;4;4;2;7;0;0;0;0;3;94;104;198
11.01.2017;10:00;0;61;6;6;4;9;0;0;1;0;0;1;80;10;3;3;7;0;0;0;0;0;87;104;191
11.01.2017;11:00;0;83;6;3;1;5;0;0;2;0;0;0;63;6;4;4;5;0;0;0;0;0;100;82;182
11.01.2017;12:00;1;96;7;3;3;13;0;1;0;0;0;0;91;8;4;2;13;0;0;0;0;1;124;119;243
11.01.2017;13:00;0;86;11;4;1;3;0;0;2;0;0;0;80;7;3;2;5;0;0;1;0;0;107;98;205
11.01.2017;14:00;0;100;6;3;2;6;0;0;1;0;0;0;92;4;2;4;8;0;0;1;0;0;118;111;229
11.01.2017;15:00;0;109;8;1;2;4;0;0;1;0;0;1;100;4;2;2;4;0;0;3;0;0;125;116;241
11.01.2017;16:00;0;118;4;2;1;4;0;0;4;0;0;0;108;7;1;2;6;0;1;3;0;0;133;128;261
11.01.2017;17:00;0;117;5;4;2;6;0;1;1;0;0;0;116;5;2;2;4;0;0;0;0;0;136;129;265
11.01.2017;18:00;0;93;2;2;0;0;0;1;0;0;0;0;97;4;1;0;0;0;0;2;0;0;98;104;202
11.01.2017;19:00;0;51;3;0;0;0;0;0;1;0;0;0;56;3;0;0;0;0;1;0;0;0;55;60;115
11.01.2017;20:00;0;36;3;0;0;0;0;0;1;0;0;0;54;0;1;0;0;0;0;0;0;0;40;55;95
11.01.2017;21:00;0;38;0;0;0;0;0;0;1;0;0;0;35;1;0;0;0;0;0;3;0;1;39;40;79
11.01.2017;22:00;0;21;1;0;0;0;0;0;1;0;0;0;19;0;0;0;0;0;0;1;0;0;23;20;43
11.01.2017;23:00;0;16;2;0;0;0;0;1;0;0;0;0;13;0;1;0;0;0;0;0;0;0;19;14;33
12.01.2017;00:00;0;11;0;0;0;0;0;0;0;0;0;0;11;0;0;0;0;0;0;1;0;0;11;12;23
12.01.2017;01:00;0;7;0;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;7;2;9
12.01.2017;02:00;0;1;0;0;0;0;0;1;0;0;0;0;1;0;0;0;0;0;0;0;0;0;2;1;3
12.01.2017;03:00;0;0;0;0;0;0;0;0;0;0;0;0;3;0;0;0;0;0;0;0;0;0;0;3;3
12.01.2017;04:00;0;0;1;0;0;0;0;0;1;0;0;0;1;0;0;0;0;0;0;0;0;0;2;1;3
12.01.2017;05:00;0;5;0;0;0;0;0;0;3;0;0;0;3;0;0;0;0;0;0;1;0;0;8;4;12
12.01.2017;06:00;0;6;0;0;0;0;0;0;0;0;0;0;8;2;0;0;0;0;1;1;0;0;6;12;18
12.01.2017;07:00;0;43;2;0;0;0;0;0;2;0;0;0;57;7;1;0;2;0;0;4;0;0;47;71;118
12.01.2017;08:00;0;74;5;3;1;6;0;0;0;0;0;0;103;6;4;1;6;0;0;2;0;1;89;123;212
12.01.2017;09:00;0;66;7;6;3;4;0;0;1;0;0;0;68;7;3;1;6;0;0;1;0;2;87;88;175
12.01.2017;10:00;0;68;6;7;4;7;0;1;0;0;0;0;76;10;4;3;5;0;0;0;0;0;93;98;191
12.01.2017;11:00;0;86;8;4;2;4;0;0;2;0;1;0;70;6;6;2;5;0;0;0;0;1;107;90;197
12.01.2017;12:00;0;97;6;4;3;9;0;0;0;0;0;0;66;4;3;4;10;0;0;1;0;1;119;89;208
12.01.2017;13:00;0;76;12;7;3;6;0;0;1;0;0;0;83;5;5;7;6;0;1;1;0;0;105;108;213
12.01.2017;14:00;0;97;3;4;5;6;0;0;1;0;1;0;100;12;2;4;7;1;0;1;0;1;117;128;245
12.01.2017;15:00;1;105;8;4;1;4;0;0;1;0;0;0;94;2;2;1;5;0;0;3;0;0;124;107;231
12.01.2017;16:00;0;131;7;1;2;4;1;0;4;0;0;0;110;8;1;2;5;0;0;2;0;1;150;129;279
12.01.2017;17:00;0;118;5;3;2;5;0;2;1;0;0;0;100;4;2;1;5;0;1;1;0;0;136;114;250
12.01.2017;18:00;0;73;5;3;0;0;0;2;0;0;0;0;87;6;0;0;1;0;0;2;0;0;83;96;179
12.01.2017;19:00;0;43;4;2;1;0;0;0;1;0;0;0;63;0;0;0;0;0;0;1;0;0;51;64;115
12.01.2017;20:00;0;48;2;0;0;0;0;0;1;0;0;0;51;1;1;0;1;1;0;1;0;0;51;56;107
12.01.2017;21:00;0;37;0;1;0;0;0;0;2;0;0;0;33;1;0;0;0;0;0;4;0;0;40;38;78
12.01.2017;22:00;0;22;0;0;0;0;0;0;1;0;0;0;28;1;0;0;0;0;0;0;0;0;23;29;52
12.01.2017;23:00;0;25;0;0;0;0;0;0;2;0;0;0;13;0;1;0;0;0;0;0;0;0;27;14;41
13.01.2017;00:00;0;9;1;0;0;0;0;0;0;0;0;0;9;0;0;0;0;0;0;0;0;0;10;9;19
13.01.2017;01:00;0;9;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;9;1;10
13.01.2017;02:00;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;1
13.01.2017;03:00;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;1;1
13.01.2017;04:00;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;1
13.01.2017;05:00;0;4;0;1;0;0;0;0;3;0;0;0;6;1;0;0;0;0;0;1;0;0;8;8;16
13.01.2017;06:00;0;6;0;0;0;0;0;0;1;0;0;0;6;1;0;1;0;0;0;1;0;0;7;9;16
13.01.2017;07:00;0;38;1;0;0;0;0;0;2;0;0;0;60;10;0;0;1;0;0;3;0;0;41;74;115
13.01.2017;08:00;0;60;2;0;1;2;0;0;0;0;1;0;94;6;4;2;8;1;0;2;0;1;66;118;184
13.01.2017;09:00;0;65;10;1;3;6;0;0;1;0;1;0;64;3;6;0;7;0;0;0;0;0;87;80;167
13.01.2017;10:00;0;79;10;3;3;6;0;0;1;0;0;0;78;12;5;1;7;0;0;0;0;0;102;103;205
13.01.2017;11:00;0;97;9;3;8;7;0;0;2;0;0;0;77;13;4;1;9;0;0;0;0;0;126;104;230
13.01.2017;12:00;0;80;9;5;6;9;0;0;0;0;0;0;86;5;4;3;10;0;0;2;0;0;109;110;219
13.01.2017;13:00;0;90;4;1;1;8;0;0;2;0;1;0;83;5;6;3;3;0;0;2;0;0;107;102;209
13.01.2017;14:00;0;99;8;2;2;5;0;0;1;0;0;0;98;7;0;1;7;0;0;1;0;1;117;115;232
13.01.2017;15:00;0;90;9;3;4;5;0;0;1;0;0;0;83;12;2;0;6;0;0;3;0;1;112;107;219
13.01.2017;16:00;0;118;9;5;0;3;0;0;4;0;0;0;97;3;2;1;4;0;0;2;0;0;139;109;248
13.01.2017;17:00;0;126;10;2;1;4;0;0;1;0;0;0;91;4;2;0;2;0;0;0;0;0;144;99;243
13.01.2017;18:00;0;72;1;2;0;0;0;0;0;0;0;0;88;4;0;0;0;0;0;1;0;0;75;93;168
13.01.2017;19:00;0;43;4;2;0;1;0;1;1;0;0;0;58;0;0;0;0;0;1;0;0;0;52;59;111
13.01.2017;20:00;0;44;3;0;0;0;0;0;2;0;0;0;52;1;0;0;0;0;0;1;0;0;49;54;103
13.01.2017;21:00;0;40;2;0;0;1;0;0;0;0;0;0;35;3;0;0;0;0;0;2;0;0;43;40;83
13.01.2017;22:00;0;17;0;1;0;0;0;0;1;0;0;0;24;2;0;0;0;0;0;1;0;0;19;27;46
13.01.2017;23:00;0;20;0;0;0;0;0;0;0;0;0;0;23;0;0;0;0;0;0;1;0;1;20;25;45
14.01.2017;00:00;0;15;1;1;0;0;0;0;2;0;0;0;11;0;1;0;0;0;0;1;0;0;19;13;32
14.01.2017;01:00;0;10;0;0;0;0;0;0;1;0;0;0;3;1;0;0;0;0;0;0;0;0;11;4;15
14.01.2017;02:00;0;2;1;0;0;0;0;0;0;0;0;0;5;0;0;0;0;0;0;0;0;0;3;5;8
14.01.2017;03:00;0;3;0;0;0;0;0;0;0;0;0;0;2;0;0;0;0;0;0;0;0;0;3;2;5
14.01.2017;04:00;0;1;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;2;0;2
14.01.2017;05:00;0;2;1;0;0;0;0;0;2;0;0;0;3;1;0;0;0;0;0;0;0;0;5;4;9
14.01.2017;06:00;0;3;0;0;0;0;0;0;1;0;0;0;8;2;0;0;0;0;0;2;0;0;4;12;16
14.01.2017;07:00;0;14;0;1;0;0;0;0;0;0;0;0;38;2;1;0;0;0;0;2;0;0;15;43;58
14.01.2017;08:00;0;30;1;4;0;2;0;1;0;0;0;0;46;3;3;0;5;0;0;2;0;0;38;59;97
14.01.2017;09:00;0;47;9;2;1;8;0;0;1;0;0;0;40;6;4;0;6;0;0;1;0;0;68;57;125
14.01.2017;10:00;0;68;6;3;3;9;0;1;0;0;1;0;55;2;1;2;7;0;1;0;0;1;91;69;160
14.01.2017;11:00;0;81;5;3;2;5;0;0;1;0;0;0;50;3;5;2;6;0;0;0;0;0;97;66;163
14.01.2017;12:00;0;79;5;2;4;3;0;0;0;0;0;0;72;3;2;0;10;1;1;0;0;0;93;89;182
14.01.2017;13:00;0;94;5;3;1;3;1;0;1;0;0;0;74;3;1;0;3;0;1;0;0;0;108;82;190
14.01.2017;14:00;0;79;3;2;2;9;0;0;1;0;0;0;99;2;0;2;11;0;0;2;0;0;96;116;212
14.01.2017;15:00;0;96;3;1;3;8;0;0;0;0;1;0;82;1;1;1;6;0;0;1;0;0;112;92;204
14.01.2017;16:00;2;99;3;1;0;0;0;0;2;0;0;2;90;7;0;1;1;1;0;0;0;0;107;102;209
14.01.2017;17:00;0;92;2;1;0;2;0;0;2;0;0;0;99;5;1;0;0;0;0;0;0;0;99;105;204
14.01.2017;18:00;0;70;2;2;0;0;0;0;0;0;0;0;86;1;2;0;0;0;0;0;0;0;74;89;163
14.01.2017;19:00;0;63;1;0;0;0;0;0;1;0;0;0;70;1;0;0;0;0;0;1;0;0;65;72;137
14.01.2017;20:00;0;40;0;0;1;1;0;0;0;0;0;0;45;0;0;0;0;0;0;0;0;0;42;45;87
14.01.2017;21:00;0;40;0;0;0;0;0;0;1;0;0;0;28;3;1;0;0;0;0;2;0;0;41;34;75
14.01.2017;22:00;0;27;0;0;0;0;0;1;1;0;0;0;31;0;0;0;0;0;0;1;0;0;29;32;61
14.01.2017;23:00;0;24;0;0;0;0;0;0;2;0;0;0;18;0;0;0;0;0;0;0;0;0;26;18;44
The issue you have, if you look at your data closely, is actually that you have several values per day (i.e., hourly data), but as.Date() produces only a date with no time. Therefore, the figure you're plotting is correct because you're plotting all values on the same day together.
To fix this, you need to produce a column of the hourly (date and time) data to use in your plots. I am particularly fond of the lubridate package, so I've used that here to help:
library(lubridate)
library(ggplot2)
daily_data$datetime <- paste(daily_data$datum, daily_data$sat, sep=" ")
daily_data$Date <- dmy_hm(daily_data$datetime)
ggplot(daily_data, aes(Date, total)) +
geom_line() +
ylab("opterecenje") +
xlab("") # total is the column with numbers that represent some data for a given date.
I'm trying to create a map with the spatstat package of R so that the sidebar shows the values of the third (preferably) or fourth column of my data frame and that the colors are also reflective of that third (or fourth column) chosen.
My script:
x=c(6.839887, 6.671494, 6.651083, 6.655289, 6.591903, 6.653641, 6.661709, 6.671664, 6.660044, 6.624659, 6.648162, 6.536877, 6.654134, 6.674678,6.618935, 6.677705, 6.643918, 6.644119, 6.670517, 6.583619, 6.649991, 6.647649, 6.656308, 6.645772, 6.648740, 6.643103, 6.652199, 6.666641,6.633400, 6.621282, 6.635427, 6.646127, 6.630862, 6.657919, 6.671616, 6.622935, 6.648225, 6.676911, 6.640234, 6.719334, 6.653202, 6.656747,6.724692, 6.639747, 6.630575, 6.657916, 6.618957, 6.640006, 6.645280, 6.614058, 6.576136, 6.631994, 6.617391, 6.782351, 6.620072, 6.661061,6.597216, 6.648755, 6.618436, 6.659507, 6.653993, 6.663255, 6.630893, 6.656322, 6.617265, 6.649022, 6.629346, 6.595224, 6.540263, 6.623435,6.652709, 6.608565, 6.618335, 6.645100, 6.790914, 6.643620, 6.462808, 6.680115, 6.716004, 6.668781, 6.765199, 6.674251, 6.647542, 6.724564,6.724556)
y=c(17.16749, 17.16727, 17.16678, 17.16673, 17.16813, 17.16663, 17.16652, 17.16636, 17.16629, 17.16856, 17.16521, 17.16519, 17.17002, 17.16465,17.17015, 17.16407, 17.16356, 17.17122, 17.16334, 17.17152, 17.16282, 17.16278, 17.16272, 17.17257, 17.16198, 17.17279, 17.16169, 17.16161,17.16146, 17.17352, 17.17389, 17.16076, 17.17420, 17.16046, 17.15917, 17.17571, 17.15895, 17.15881, 17.15860, 17.15827, 17.15797, 17.15776,17.17761, 17.15664, 17.15622, 17.15610, 17.15571, 17.15561, 17.15527,17.15514, 17.15494, 17.15447, 17.15438, 17.18041, 17.18053, 17.15402,17.18090, 17.15384, 17.18121, 17.15355, 17.15352, 17.15349, 17.18213,17.15242, 17.15201, 17.14978, 17.18591, 17.18688, 17.18707, 17.18761,17.14712, 17.18788, 17.18794, 17.14619, 17.18868, 17.14588, 17.14511,17.14471, 17.14440, 17.14430, 17.19116, 17.19140, 17.14222, 17.14123,17.33627)
z=c(32.23228,526.46061, -1300.03539, -376.04329, 139.67322,-913.24800, -526.46061, 354.55511, 483.48424, 161.16141, 182.64960, 419.0196, 75.20866, -225.62598, -1536.40546, -397.53148, -1106.64169, -440.50786, 118.18504,-290.09054, -1471.94089, 440.50786,-848.78343, -1385.98814, -676.87793, -1622.35821, -1450.45271,75.20866, -1557.89365, 161.16141, 376.04329, 354.55511, -32.23228,-1171.10626,-75.20866, 547.94880, -805.80706, 870.27162, -698.36612,-32.23228, -2331.46842, -182.64960, 75.20866, -719.85431,-1837.24009,913.24800, -1106.64169, 698.36612, 483.48424, -676.87793, -3019.09045, 891.75981, 1106.64169, 333.06692, -913.24800,333.06692, 934.73619, 354.55511, 75.20866, -891.75981, -247.11416, -1966.16922, 139.67322, -784.31887, -569.43699, -118.18504,-440.50786, 397.53148, -655.38974, 139.67322, 53.72047, -633.90155,-633.90155, 419.01967, -547.94880, 75.20866, 569.43699, 290.09054, -376.04329, 547.94880, 75.20866, -10.74409, 182.64960,-397.53148, -479.53833 )
w=c(96326.91, 96769.46, 95127.94, 95960.41, 96423.22, 95476.93, 95825.18,96615.67, 96731.03, 96442.47, 96461.73, 96673.36, 96365.44, 96095.53,94914.31, 95941.10, 95302.53, 95902.47, 96403.96, 96037.64, 94972.60,96692.58, 95535.03, 95050.29, 95689.84, 94836.56, 94992.03, 96365.44,94894.87, 96442.47, 96634.90, 96615.67, 96269.09, 95244.36, 96230.54,96788.68, 95573.74, 97076.62, 95670.50, 96269.09, 94193.69, 96134.12,96365.44, 95651.15, 94642.01, 97114.98, 95302.53, 96923.12, 96731.03,95689.84, 93567.91, 97095.80, 97287.46, 96596.43, 95476.93, 96596.43,97134.15, 96615.67, 96365.44, 95496.30, 96076.24, 94525.17, 96423.22,95593.10, 95786.52, 96191.98, 95902.47, 96654.13, 95709.18, 96423.22,96346.17, 95728.52, 95728.52, 96673.36, 95805.85, 96365.44, 96807.89,96557.96, 95960.41, 96788.68, 96365.44, 96288.37, 96461.73,95941.10, 99451.20)
shap.lo=data.frame(x,y,z,w)
library(spatstat)
shap.lo.win <- owin(range(shap.lo[,1]), range(shap.lo[,2]))
centroid.owin(shap.lo.win) ; area.owin(shap.lo.win)
shap.lo.ppp <- as.ppp(shap.lo[,c(1,2,3)], shap.lo.win) # making a ppp object
plot(density(shap.lo.ppp,0.02), col=topo.colors(25), main='', xlab='x',
ylab='y')
points(x, y)
the result is shown below
I would like to know why the sidebar shows different values than the ones shown in the third column of my data frame, that is, in addition to displaying no negative values, shows values much larger than those contained in the third column.
Is it possible to do this, that is, make the colors and the sidebar represent the third or fourth column of the data frame?
I thank the help of all you!
Let me quote directly from the help file for density.ppp:
This function is often misunderstood.
The result of density.ppp is not a spatial smoothing of the
marks or weights attached to the point pattern. To perform
spatial interpolation of values that were observed at the points
of a point pattern, use Smooth.ppp.
The result of density.ppp is not a probability density. It is
an estimate of the intensity function of the point process that
generated the point pattern data. Intensity is the expected
number of random points per unit area. The units of intensity are
“points per unit area”. Intensity is usually a function of
spatial location, and it is this function which is estimated by
density.ppp. The integral of the intensity function over a
spatial region gives the expected number of points falling in this
region.
So try Smooth.ppp (note the upper case S), and see if you can make that produce the results you expected.
Can someone explain why my plot apparently overwrites the first line with the second? The values are in no way the same.
per.acc<-c(0.508407842930696, 0.508407842930696, 0.508776550216615, 0.508895063272804,
0.509619309727288, 0.51039622865119, 0.510712263467692, 0.511871057794867,
0.512252933198141, 0.513280046351773, 0.514122805862446, 0.514478345031011,
0.514768043612805, 0.515676643710249, 0.516084855348231, 0.51682226992007,
0.517757206252222, 0.519008177400877, 0.519890441263613, 0.520680528304868,
0.521286261703164, 0.522181693683254, 0.523024453193927, 0.524222751873165,
0.524420273633479, 0.524986502679712, 0.525631740430071, 0.52609262453747,
0.526777366639891, 0.52733042756877, 0.52788348849765, 0.528278532018277,
0.528673575538905, 0.529661184340475, 0.530569784437919, 0.531109677249444,
0.530833146785004, 0.531425712065946, 0.531280862775049, 0.53063562502469,
0.531465216418009, 0.531215022188278, 0.531188685953569, 0.53153105700478,
0.53180758746922, 0.532110454168368, 0.532005109229533, 0.532070949816305,
0.53233431216339, 0.532189462872493, 0.532307975928682, 0.532321144046036,
0.531741746882448, 0.531504720770071, 0.530214245269354, 0.529661184340475,
0.529002778472762, 0.529845537983434, 0.53054344820321, 0.529108123411596,
0.528186355196797, 0.527896656615004, 0.526830039109309, 0.526263810063075,
0.525091847618546, 0.524354433046707, 0.523538009770743, 0.522655745908008,
0.522036844392357, 0.520812209478411, 0.519640247033882, 0.518520957058769,
0.517335826496886, 0.515531794419352, 0.514715371143388, 0.51314836517823,
0.511831553442804, 0.510975625814777, 0.510040689482625, 0.508118144348902,
0.507104199312624, 0.505813723811907, 0.505168486061548, 0.504312558433521,
0.503311781514597, 0.502284668360964, 0.501297059559395, 0.500546476870202,
0.500072424645448, 0.499018975257107, 0.497886517164641, 0.496964748949843,
0.495516256040874, 0.494120435601322, 0.493383021029483, 0.492803623865896,
0.492421748462622, 0.492118881763474, 0.492026704941994, 0.491684333890784
)
num.pred<-c(38609L, 38609L, 38637L, 38646L, 38701L, 38760L, 38784L, 38872L,
38901L, 38979L, 39043L, 39070L, 39092L, 39161L, 39192L, 39248L,
39319L, 39414L, 39481L, 39541L, 39587L, 39655L, 39719L, 39810L,
39825L, 39868L, 39917L, 39952L, 40004L, 40046L, 40088L, 40118L,
40148L, 40223L, 40292L, 40333L, 40312L, 40357L, 40346L, 40297L,
40360L, 40341L, 40339L, 40365L, 40386L, 40409L, 40401L, 40406L,
40426L, 40415L, 40424L, 40425L, 40381L, 40363L, 40265L, 40223L,
40173L, 40237L, 40290L, 40181L, 40111L, 40089L, 40008L, 39965L,
39876L, 39820L, 39758L, 39691L, 39644L, 39551L, 39462L, 39377L,
39287L, 39150L, 39088L, 38969L, 38869L, 38804L, 38733L, 38587L,
38510L, 38412L, 38363L, 38298L, 38222L, 38144L, 38069L, 38012L,
37976L, 37896L, 37810L, 37740L, 37630L, 37524L, 37468L, 37424L,
37395L, 37372L, 37365L, 37339L)
plot(x=1:100,y=per.acc, main="RSI predicting trend",ylab=NA,xlab="RSI cutoff value",col="blue",type="p")
par(new=T)
plot(x=1:100,y=num.pred,axes=F,xlab=NA,main=NA,ylab=NA,col="red",type="p")
axis(side=4)
Because your per.acc and num.pred has perfect co-linearity:
plot(per.acc, num.pred)
Therefore, your second plot and first plot just differ by a vertical shift. When you do par(new = TRUE), they coincide.
On the other hand, if you do some random shuffling:
set.seed(0); num.pred <- sample(num.pred)
the issue goes away:
I have my data in a csv file in the given format.
csv file data
Latitude,Longitude,Temperature
20,84.01,15.93913043
28.48,77.13,16.62857143
28.68,77.2,17.81333333
31.32,78.16,2.472222222
31.531,77.112,5.228
28.11,77,21.85
31.5,77.09,7.910526316
31.43,76.57,11.444
28.7,77.15,17.708
30.55,77.35,15.30526316
26.95,78.96,16.46818182
27.44,79.39,15.74090909
26.58,81.59,15.90952381
25.33,80.43,18.465
29.45,77.34,14.15238095
20.42,86.47,19.83181818
29.52,75.5,14.135
24.17,72.43,20.065
25.1,76.19,18.59444444
30.975,76.517,14.88421053
28.8,76.13,16.05
29.54,75.04,15.295
24.32,72.3,18.84782609
23.86,72.13,20.49444444
30.19,74.95,13.996
22.36,73.16,22.365
30.78,75.84,13.75652174
21.86,73.52,24.725
21.5,70.44,22.812
21.36,69.75,22.33125
30.32,78.05,15.35
20.92,72.89,21.3
17.69,74.02,23.45
28.3278,77.2467,17.87857143
20.17,79.98,21.11428571
24,76.73,23.67857143
22.76,74.59,19.97619048
22.03,74.97,20.85
29.57,80.23,12.70869565
29.95,79.9,13.425
25.92,83.56,16.67
27.6,81.58,15.128
24.68,78.4,19.1
18.11,76.02,22.84285714
20.39,78.11,24.98571429
24.57,80.82,23.35714286
21.283,76.198,22.98235294
21.81,80.18,20.16666667
24.5,81.3,16.22857143
22.09,82.17,18.93636364
30.35,76.87,14.77
19.7,81.7,19.98823529
18.9,81.35,19.16956522
28.9917,77.701,15.43636364
28.39,77.83,15.745
27.58,77.98,16.52631579
27.03,79.92,17.40526316
26.57,80.48,12.67
25.17,80.91,24.75
26.55,79.55,17.12727273
22.443,77.03,18.825
30.19,78.18,15.72857143
29.87,77.88,14.75454545
26.54,81.83,15.75416667
29.93,77.97,12.96666667
26.127,81.94,17.54666667
26.43,82.57,16.63684211
29.34,80.09,11.82631579
27.14,83.53,15.82
21.1,86.5,19.82
20.25,85.82,21.31
21.13,86.57,19.352
20.23,86.18,19.52173913
20.46,85.9,20.74
27.17,78.03,17.292
20.83,84.33,19.224
21.89,84.03,19.47142857
20.01,83.01,20.295
19,83.01,21.24285714
18.53,73.83,22.47142857
18.81,82.71,19.04545455
18.01,82.01,20.73076923
25,84.01,16.952
25.03,85.6,18.48125
19.92,83.16,19.975
26.21,84.35,16.15454545
26.58,84.38,16.1952381
25.73,85.23,15.9375
26.38,85.73,16.39
25.98,85.66,16.176
30.48,78.05,8.985
18.35,81.88,21.736
26.54,85.72,16.53043478
26,85.01,16.104
26.9,75.8,18.97272727
25.92,86.8,16.01904762
26.42,74.62,21.04545455
24.87,85.53,14.9
24.8,85.04,16.236
25.91,86.55,16.17
25.3,83,20.13333333
25.52,87.57,17.13181818
25.37,86.48,17.56190476
25.01,86,16.42727273
26.9,76.35,18.836
25.23,87.03,15.3875
26.07,76.37,17.324
25.41,75.65,18.03684211
25.12,75.93,17.93333333
19.65,78.52,21.51052632
24.58,76.13,19.8
26.16,75.78,19.24285714
24.9,74.61,18.74583333
18.56,77.88,23.68571429
18.83,78.93,22.75238095
17.01,78.97,21.63684211
18.43,79.08,22.78181818
17.21,77.58,22.36363636
26.51,85.28,15.988
18.01,79.6,24.295
16.5,78.23,22.368
28.62,77.27,17.135
15.45,78.46,24.09444444
16.85,79.47,23.58571429
16.35,80.43,23.988
17.62,74.07,22.97777778
16.36,80.84,22.976
14.68,77.67,23.62
17.87,82.34,19.12
17.7,83,21.62083333
14.11,78.15,23.52352941
14.43,78.8,23.92857143
17.24,81.1,22.6
23.55,74.44,26.72
23.51,74.38,19.105
28.28,75.12,18.33333333
12.38,78.21,23.22857143
28.35,75.58,16.25714286
11.33,76.8,15.97391304
13.15,80.23,23.65263158
13.63,79.73,22.688
12.71,77.81,21.03809524
11.62,79.53,24.01818182
13.61,79.34,22.35454545
11.01,76.94,24.765
11.37,76.63,14.824
13.25,80.31,26.384
11.13,79.07,25.87
11.22,78.87,24.86
26.25,82,16.25217391
10.77,76.7,26.32
11.15,78.15,26.03888889
10.53,76.28,27.46666667
24.57,73.7,9.8
9.62,76.42,28.135
12.95,78.25,22.26111111
23.86,73.86,20.24
12.97,77.18,23.35652174
9.55,77.934,26.18636364
10.117,77.583,24.93157895
13.663,76.916,22.21363636
9.2,77.88,26.98
8.43,76.99,26.33809524
8.73,77.7,27.02380952
11.26,75.77,28.36470588
8.43,77.3,26.47619048
11.0617,76.21,21.59
12.15,76.8,22.73809524
13.34,75.8,20.49047619
11.91,76.93,23.37619048
25.79,73.3,18.55652174
13.72,75.72,22.75263158
12.42,75.73,19.85
26.8,75.8,18.80416667
25.75,71.4,24.05
25.38,72.5,22.38
27.9,78.07,17.15789474
26.12,91.82,18.00555556
14.611,74.846,21.89
12.65,74.96,25.74736842
16.18,77.32,24.73888889
15.48,74.98,22.80909091
27.58,75.13,17.195
27.15,73.78,21.04166667
17.35,76.8,23.78888889
27.93,74.98,17.32105263
26.351,92.633,17.57142857
25.87,93.43,15.88
14.78,75.33,23.655
15.9,75.55,23.05454545
16.98,75.75,23.73333333
27.53,76.6,17.01052632
26.44,89.94,15.684
26.35,90.65,18.98
26.48,90.9,15.688
29.49,73.5,13.752
14.03,77.27,23.605
29.93,73.88,14.30588235
26.65,91.33,16.59565217
16.15,75.6,19.57692308
13.96,75.56,22.95714286
26.71,93.13,15.95714286
26.72,77.89,19.03846154
25.11,76.47,19.04444444
26.16,90.63,16.86818182
26.51,93.96,16.56956522
26.467,90.283,16.615
27.2,77.46,15.84545455
26.72,94.18,16.61
26.55,77,18.58695652
27.467,95.02,16.452
10.08,77.05,16.04
25.36,74.63,19.18571429
14.283,74.45,25.65238095
24.6,92.55,18.965
25.56,91.86,11.03478261
24.85,92.3,23.05
25.454,92.2069,12.044
9.834,76.95,22.456
28.333,79.417,15.99090909
27.71,81.91,14.88571429
25.87,84.13,12.81666667
26.78,78.98,20.1
26.33,79.97,17.7125
25.68,91.92,14.98888889
25.2,90.63,19.11578947
25.52,91.27,10.316
26.25,81.37,17.04210526
24.58,72.7,14.35
26.95,94.62,17.43333333
27.48,95.31,16.108
23.53,84.81,16.07619048
26.97,93.87,15.712
25.5,90.58,16.09473684
23.78,85.88,17.42631579
24.83,87.2,17.484
25.23,86.65,22.075
27.47,94.55,17.31052632
24.5,87.81,17.99130435
26.52,90.48,16.925
24.49,86.66,21.56875
17.43,78.47,23.636
24.2,86.3,16.83
24.16,83.8,16.635
22.47,88.09,18.27777778
24.8,92.78,18.168
22.69,86.43,17.15238095
22.27,87.92,19.71666667
29.46,79.65,9.25
22.3,87.3,18.97
22.36,88.43,18.58571429
22.662,88.873,19.12
22.97,88.48,21.24545455
23.4,88.49,17.99166667
23.53,91.46,17.92
32.37,75.55,14.44583333
23.21,87.88,19.11666667
32.65,74.8,13.45454545
23.24,87.04,18.892
23.91,87.52,19.42272727
31.18,77.13,4.424
24.092,88.27,21.4
25.62,88.14,17.50555556
31.33,75.33,13.91428571
23.8,91.27,19.10555556
26.52,88.72,13.65625
24.4,89.39,15.97826087
34.05,74.8,3.419047619
34.14,74.87,2.961904762
25.31,88.76,17.65263158
34.06,74.71,3.036363636
27.07,88.46,13.75652174
27.15,88.4,11.696
24.07,91.61,18.56842105
34.35,74.4,1.118181818
25.18,93.01,18.8
33.888,74.9808,5.125
33.64,74.94,1.635
23.43,85.3,17.4
24.18,91.82,22.1
23.48,93.3,14.836
24.208,92.675,18.988
23.95,92.49,18.305
33.22,75.26,7.57
23.33,92.83,17.485
22.5055,92.895,18.21578947
32.9,75.167,11.42
22.95,92.93,16.78
22.48,92.98,19.40454545
26.72,92.08,16.804
18.22,74.48,23.43333333
17.26,74.42,17.07647059
30.15,75.41,17.36666667
31.63,74.83,14.29444444
31.46,77.22,5.661904762
30.73,76.8,15.395
31.32,75.3,13.80769231
30.91,74.66,13.28947368
20.76,73.69,23.83333333
31.37,74.97,12.724
31.38,76.024,14.865
20.4,72.833,20.672
23.303,70.36,20.628
31.39,75.36,14.52142857
21.144,72.75,24.92222222
23.37,68.54,19.37222222
20.8,70.7,22.295
20.143,74.794,26.94444444
18.97,73.04,24.73157895
18.33,73.79,22.23684211
21.02,75.53,23.15833333
28.47,77.03,21.1
17.67,75.92,24.94545455
19.03,73.17,24.74615385
22.94,79.22,17.4
25.688,78.4419,18.9
22.94,81.08,15.76818182
22.6,82.13,18.308
21.1,81.034,21.95
21.97,81.24,20.3375
22.37,82.733,17.77083333
21.23,81.7,14.96428571
25.35,81.3,16.92
26.79,82.72,18.1125
26.48,83.77,16.86666667
29.36,79.46,9.733333333
25.92,84.13,16.616
28.98,79.4,16.068
25.72,82.68,16.572
25.17,82.6,16.175
29.62,79.62,13.245
30.7,76.9,16.096
30.36,75.53,13.78888889
22.8,73.22,20.07083333
22.58,72.28,20.46
22.73,72.88,23.5875
23.27,72.65,21.43888889
30.36,76.45,13.97
22.88,74.23,21.42
17.69,74.02,20.988
18.53,73.8,21.98095238
28.1,77,14.965
28.08,76.59,17.8
31.1,77.17,6.628571429
28.28,76.15,17.94
30.1,77.28,14.07777778
29.68,77,13.628
29.7,77.03,14.21764706
29.38,77,12.85384615
28.611,76.651,15.08947368
28.8836,76.619,15.695
29.78,76.4,13.67
31.35,76.7525,13.14347826
31.45,76.26,13.85238095
28.66,77.31,17.71578947
29.58,76.51,14.625
29.3,75.45,15.28421053
29.068,75.476,16.05263158
28.63,77.15,16.575
28.8,76.97,15.12
30.90444,76.96,8.48
30.89916,76.96,12.392
31.5858,77.0748,6.884
32.3521,76.0608,7.844
32.1994,76.3249,11.67777778
28.58,76.9,15.6
29.13,75.7,15.608
30.2,75.88,14.72857143
31.09,76.38,15.6
31.13,76.12,14.7
30.6,76.1,13.49090909
30.67,74.74,13.825
21.7,73.01,14.45
21.73,72.55,23.42352941
32.05,75.42,13.04736842
30.62,76.41,13.49090909
31.6,75.98,16.1
20.54,73.16,22.70833333
22.73,71.62,24.21538462
28.63,77,17.245
19.18,72.98,27.00666667
20.9,74.8,23.6
23.83,73,21.71333333
21.16,73.79,23.51578947
28.63,77.36,16.836
28.63,77.37,17.252
28.23,77.18,15.85454545
22.34,70.91,21.075
21.37,74.23,24.24210526
19.09,74.75,21.21818182
19.36,75.95,23.3047619
16,73.7,25.14166667
16.67,74.16,22.89545455
18.4,76.56,23.275
15.48,73.91,26.81764706
19.85,75.87,21
20.94,77.78,23.285
18.72,76.36,23.89375
19.73,77.15,21.85263158
20.13,77.13,25.08125
20.7,77.65,20.05
21.16,79.65,21.32
21.49,80.15,18.30416667
22.3,79.67,23.19090909
26.4772,77.9866,17.328
26.5727,78.81,16.856
22.0683,79.5488,16.036
24.5972,77.75,18.8625
24.42,74.8708,18.93809524
24.06,75.1,18.93333333
22.88,84.56,16.75
23.444,75.044,18.97916667
22.7,75.89,20.97
20.53,76.19,23.49
21.82,75.62,22.50952381
23.83,78.69,26.7875
23.18,77.06,20.295
25.03,79.5,18.79
23.82,79.4,17.095
25.67,76.69,21.26
22.314,74.353,19.064
22.056,78.94,16.83333333
21.86,77.937,19.58333333
23.21,79.96,17.56190476
22.6,80.3,16.14761905
23.28,81.33,16.37142857
24.73,80.21,19.6625
23.525,80.84,17.63888889
24.41,81.86,18.38461538
23.09,83.14,15.60666667
24.07,82.61,16.95555556
23.12,81.7,17.71111111
21.88,83.4,19.548
23.515,77.81,20.57142857
22.58,81.13,19.17368421
21.103,82.083,20.92380952
20.66,81.53,20.27647059
20.25,81.5,20.64444444
30.756,75.163,12.795
18.8,80.8,21.0952381
28.01,79.11,16.25
27.55,78.65,15.1
27.22,79.03,17.38947368
16.98,73.32,24.70526316
30.05,79,8.456
30.74,76.73,16.32727273
25.21,82.27,22.54
19.8,85.83,21.5952381
20.84,86.32,20.155
22.03,84.05,19.05217391
21.33,83.62,19.915
20.12,85.08,19.904
21.61,85.56,18.728
20.83,85.1,19.656
21.36,83.88,18.896
20.52,86.43,19.76190476
19.3,84.8,22.65333333
19.25,82.54,19.03333333
24.95,84,18.05
23.55,83.96,16.288
25.0286,73.89,18.44117647
11.28,76.23,25.91666667
26.26,72.99,18.52
15.42,75.63,23.7
22.9,88.37,16.58095238
27.28,88.23,10.93333333
26.02,94.53,10.45
I tried drawing a contour using contour(), contourplot() etc, I am not sure how to make my data regular to get a contour . I tried akima interp() but getting error given below. not sure what that means and how to correct my data to get a contour.
*Error in is.finite(x) : default method not implemented for type 'list'*
someone please help me out. Thanks.
This works for me (by the way, dumping this size data file, even though it's only moderate-sized, is kind of a pain for responders: if you can post it somewhere that's easier. I had to save the HTML page source and extract the info that way, as it was too big for me to cut & paste from the clipboard):
X <- read.csv("tempcontour.txt")
library(akima)
## needed to specify the 'duplicate' argument (see ?interp)
interp1 <- with(X,interp(Latitude,Longitude,Temperature,
duplicate="mean"))
names(interp1) <- colnames(X) ## assumes colnames are in Lat/Long/Temp order
png("interp1.png")
with(interp1,contour(Longitude,Latitude,Temperature,
ylab="Latitude",xlab="Longitude"))
dev.off()