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()
Related
I'm currently running a time series script in R-Markdown where I pass the values of Percent Use and another vector of the time values. I am processing the two separate vectors with the regular c() function within r. The percentage vector is able to be passed through when compilining like normal, however I am running into trouble with the date/time vector. The length of both vectors are 749, the percentage vector just has values 0-100 passed within them. THe date/time vector has strings passed into them as such:
dt=c('2022-06-19 14:05:00.0','2022-06-19 14:06:00.0', ....
If I only pass a few arguments into the dt vector, it will compile regularly, however, once I increase the size to around half of what it needs to be I start getting the following error:
Error: Incomplete expression: dt=c('2022-06-19 12:40:00.0','2022-06-19 12:41:00.0','2022-06-19 12:42:00.0','2022-06-19 12:43:00.0','2022-06-19 12:44:00.0','2022-06-19 12:45:00.0','2022-06-19 12:46:00.0','2022-06-19 12:47:00.0','2022-06-19 12:48:00.0','2022-06-19 12:49:00.0','2022-06-19 12:50:00.0','2022-06-19 12:51:00.0','2022-06-19 12:52:00.0','2022-06-19 12:53:00.0','2022-06-19 12:54:00.0','2022-06-19 12:55:00.0','2022-06-19 12:56:00.0','2022-06-19 12:57:00.0','2022-06-19 12:58:00.0','2022-06-19 12:59:00.0','2022-06-19 13:00:00.0','2022-06-19 13:01:00.0','2022-06-19 13:02:00.0','2022-06-19 13:03:00.0','2022-06-19 13:04:00.0','2022-06-19 13:05:00.0','2022-06-19 13:06:00.0','2022-06-19 13:07:00.0','2022-06-19 13:08:00.0','2022-06-19 13:09:00.0','2022-06-19 13:10:00.0','2022-06-19 13:11:00.0','2022-06-19 13:12:00.0','2022-06-19 13:13:00.0','2022-06-19 13:14:00.0','2022-06-19 13:15:00.0','2022-06-19 13:16:00.0','2022-06-19 13:17:00.0','2022-06-19 13:18:00.0','2022-06-19 13:19:00.0','2022
At first I believed it could be a parenthesis in the wrong place, however, there is no mistakes with that. I've looked at other articles with somewhat similar issues and have seen a concept of a maximum size vector allowed, however the percentage vector was able to pass all 700. Is there a way to bypass this error, I feel that it is a memory/storage issue with R.
The full code is a lot but it is:
dt=c('2022-06-19 12:40:00.0','2022-06-19 12:41:00.0','2022-06-19 12:42:00.0','2022-06-19 12:43:00.0','2022-06-19 12:44:00.0','2022-06-19 12:45:00.0','2022-06-19 12:46:00.0','2022-06-19 12:47:00.0','2022-06-19 12:48:00.0','2022-06-19 12:49:00.0','2022-06-19 12:50:00.0','2022-06-19 12:51:00.0','2022-06-19 12:52:00.0','2022-06-19 12:53:00.0','2022-06-19 12:54:00.0','2022-06-19 12:55:00.0','2022-06-19 12:56:00.0','2022-06-19 12:57:00.0','2022-06-19 12:58:00.0','2022-06-19 12:59:00.0','2022-06-19 13:00:00.0','2022-06-19 13:01:00.0','2022-06-19 13:02:00.0','2022-06-19 13:03:00.0','2022-06-19 13:04:00.0','2022-06-19 13:05:00.0','2022-06-19 13:06:00.0','2022-06-19 13:07:00.0','2022-06-19 13:08:00.0','2022-06-19 13:09:00.0','2022-06-19 13:10:00.0','2022-06-19 13:11:00.0','2022-06-19 13:12:00.0','2022-06-19 13:13:00.0','2022-06-19 13:14:00.0','2022-06-19 13:15:00.0','2022-06-19 13:16:00.0','2022-06-19 13:17:00.0','2022-06-19 13:18:00.0','2022-06-19 13:19:00.0','2022-06-19 13:20:00.0','2022-06-19 13:21:00.0','2022-06-19 13:22:00.0','2022-06-19 13:23:00.0','2022-06-19 13:24:00.0','2022-06-19 13:25:00.0','2022-06-19 13:26:00.0','2022-06-19 13:27:00.0','2022-06-19 13:28:00.0','2022-06-19 13:29:00.0','2022-06-19 13:30:00.0','2022-06-19 13:31:00.0','2022-06-19 13:32:00.0','2022-06-19 13:33:00.0','2022-06-19 13:34:00.0','2022-06-19 13:35:00.0','2022-06-19 13:36:00.0','2022-06-19 13:37:00.0','2022-06-19 13:38:00.0','2022-06-19 13:39:00.0','2022-06-19 13:40:00.0','2022-06-19 13:41:00.0','2022-06-19 13:42:00.0','2022-06-19 13:43:00.0','2022-06-19 13:44:00.0','2022-06-19 13:45:00.0','2022-06-19 13:46:00.0','2022-06-19 13:47:00.0','2022-06-19 13:48:00.0','2022-06-19 13:49:00.0','2022-06-19 13:50:00.0','2022-06-19 13:51:00.0','2022-06-19 13:52:00.0','2022-06-19 13:53:00.0','2022-06-19 13:54:00.0','2022-06-19 13:55:00.0','2022-06-19 13:56:00.0','2022-06-19 13:57:00.0','2022-06-19 13:58:00.0','2022-06-19 13:59:00.0','2022-06-19 14:00:00.0','2022-06-19 14:01:00.0','2022-06-19 14:02:00.0','2022-06-19 14:03:00.0','2022-06-19 14:04:00.0','2022-06-19 14:05:00.0','2022-06-19 14:06:00.0','2022-06-19 14:07:00.0','2022-06-19 14:08:00.0','2022-06-19 14:09:00.0','2022-06-19 14:10:00.0','2022-06-19 14:11:00.0','2022-06-19 14:12:00.0','2022-06-19 14:13:00.0','2022-06-19 14:14:00.0','2022-06-19 14:15:00.0','2022-06-19 14:16:00.0','2022-06-19 14:17:00.0','2022-06-19 14:18:00.0','2022-06-19 14:19:00.0','2022-06-19 14:20:00.0','2022-06-19 14:21:00.0','2022-06-19 14:22:00.0','2022-06-19 14:23:00.0','2022-06-19 14:24:00.0','2022-06-19 14:25:00.0','2022-06-19 14:26:00.0','2022-06-19 14:27:00.0','2022-06-19 14:28:00.0','2022-06-19 14:29:00.0','2022-06-19 14:30:00.0','2022-06-19 14:31:00.0','2022-06-19 14:32:00.0','2022-06-19 14:33:00.0','2022-06-19 14:34:00.0','2022-06-19 14:35:00.0','2022-06-19 14:36:00.0','2022-06-19 14:37:00.0','2022-06-19 14:38:00.0','2022-06-19 14:39:00.0','2022-06-19 14:40:00.0','2022-06-19 14:41:00.0','2022-06-19 14:42:00.0','2022-06-19 14:43:00.0','2022-06-19 14:44:00.0','2022-06-19 14:45:00.0','2022-06-19 14:46:00.0','2022-06-19 14:47:00.0','2022-06-19 14:48:00.0','2022-06-19 14:49:00.0','2022-06-19 14:50:00.0','2022-06-19 14:51:00.0','2022-06-19 14:52:00.0','2022-06-19 14:53:00.0','2022-06-19 14:54:00.0','2022-06-19 14:55:00.0','2022-06-19 14:56:00.0','2022-06-19 14:57:00.0','2022-06-19 14:58:00.0','2022-06-19 14:59:00.0','2022-06-19 15:00:00.0','2022-06-19 15:01:00.0','2022-06-19 15:02:00.0','2022-06-19 15:03:00.0','2022-06-19 15:04:00.0','2022-06-19 15:05:00.0','2022-06-19 15:06:00.0','2022-06-19 15:07:00.0','2022-06-19 15:08:00.0','2022-06-19 15:09:00.0','2022-06-19 15:10:00.0','2022-06-19 15:11:00.0','2022-06-19 15:12:00.0','2022-06-19 15:13:00.0','2022-06-19 15:14:00.0','2022-06-19 15:15:00.0','2022-06-19 15:16:00.0','2022-06-19 15:17:00.0','2022-06-19 15:18:00.0','2022-06-19 15:19:00.0','2022-06-19 15:20:00.0','2022-06-19 15:21:00.0','2022-06-19 15:22:00.0','2022-06-19 15:23:00.0','2022-06-19 15:24:00.0','2022-06-19 15:25:00.0','2022-06-19 15:26:00.0','2022-06-19 15:27:00.0','2022-06-19 15:28:00.0','2022-06-19 15:29:00.0','2022-06-19 15:30:00.0','2022-06-19 15:31:00.0','2022-06-19 15:32:00.0','2022-06-19 15:33:00.0','2022-06-19 15:34:00.0','2022-06-19 15:35:00.0','2022-06-19 15:36:00.0','2022-06-19 15:37:00.0','2022-06-19 15:38:00.0','2022-06-19 15:39:00.0','2022-06-19 15:40:00.0','2022-06-19 15:41:00.0','2022-06-19 15:42:00.0','2022-06-19 15:43:00.0','2022-06-19 15:44:00.0','2022-06-19 15:45:00.0','2022-06-19 15:46:00.0','2022-06-19 15:47:00.0','2022-06-19 15:48:00.0','2022-06-19 15:49:00.0','2022-06-19 15:50:00.0','2022-06-19 15:51:00.0','2022-06-19 15:52:00.0','2022-06-19 15:53:00.0','2022-06-19 15:54:00.0','2022-06-19 15:55:00.0','2022-06-19 15:56:00.0','2022-06-19 15:57:00.0','2022-06-19 15:58:00.0','2022-06-19 15:59:00.0','2022-06-19 16:00:00.0','2022-06-19 16:01:00.0','2022-06-19 16:02:00.0','2022-06-19 16:03:00.0','2022-06-19 16:04:00.0','2022-06-19 16:05:00.0','2022-06-19 16:06:00.0','2022-06-19 16:07:00.0','2022-06-19 16:08:00.0','2022-06-19 16:09:00.0','2022-06-19 16:10:00.0','2022-06-19 16:11:00.0','2022-06-19 16:12:00.0','2022-06-19 16:13:00.0','2022-06-19 16:14:00.0','2022-06-19 16:15:00.0','2022-06-19 16:16:00.0','2022-06-19 16:17:00.0','2022-06-19 16:18:00.0','2022-06-19 16:19:00.0','2022-06-19 16:20:00.0','2022-06-19 16:21:00.0','2022-06-19 16:22:00.0','2022-06-19 16:23:00.0','2022-06-19 16:24:00.0','2022-06-19 16:25:00.0','2022-06-19 16:26:00.0','2022-06-19 16:27:00.0','2022-06-19 16:28:00.0','2022-06-19 16:29:00.0','2022-06-19 16:30:00.0','2022-06-19 16:31:00.0','2022-06-19 16:32:00.0','2022-06-19 16:33:00.0','2022-06-19 16:34:00.0','2022-06-19 16:35:00.0','2022-06-19 16:36:00.0','2022-06-19 16:37:00.0','2022-06-19 16:38:00.0','2022-06-19 16:39:00.0','2022-06-19 16:40:00.0','2022-06-19 16:41:00.0','2022-06-19 16:42:00.0','2022-06-19 16:43:00.0','2022-06-19 16:44:00.0','2022-06-19 16:45:00.0','2022-06-19 16:46:00.0','2022-06-19 16:47:00.0','2022-06-19 16:48:00.0','2022-06-19 16:49:00.0','2022-06-19 16:50:00.0','2022-06-19 16:51:00.0','2022-06-19 16:52:00.0','2022-06-19 16:53:00.0','2022-06-19 16:54:00.0','2022-06-19 16:55:00.0','2022-06-19 16:56:00.0','2022-06-19 16:57:00.0','2022-06-19 16:58:00.0','2022-06-19 16:59:00.0','2022-06-19 17:00:00.0','2022-06-19 17:01:00.0','2022-06-19 17:02:00.0','2022-06-19 17:03:00.0','2022-06-19 17:04:00.0','2022-06-19 17:05:00.0','2022-06-19 17:06:00.0','2022-06-19 17:07:00.0','2022-06-19 17:08:00.0','2022-06-19 17:09:00.0','2022-06-19 17:10:00.0','2022-06-19 17:11:00.0','2022-06-19 17:12:00.0','2022-06-19 17:13:00.0','2022-06-19 17:14:00.0','2022-06-19 17:15:00.0','2022-06-19 17:16:00.0','2022-06-19 17:17:00.0','2022-06-19 17:18:00.0','2022-06-19 17:19:00.0','2022-06-19 17:20:00.0','2022-06-19 17:21:00.0','2022-06-19 17:22:00.0','2022-06-19 17:23:00.0','2022-06-19 17:24:00.0','2022-06-19 17:25:00.0','2022-06-19 17:26:00.0','2022-06-19 17:27:00.0','2022-06-19 17:28:00.0','2022-06-19 17:29:00.0','2022-06-19 17:30:00.0','2022-06-19 17:31:00.0','2022-06-19 17:32:00.0','2022-06-19 17:33:00.0','2022-06-19 17:34:00.0','2022-06-19 17:35:00.0','2022-06-19 17:36:00.0','2022-06-19 17:37:00.0','2022-06-19 17:38:00.0','2022-06-19 17:39:00.0','2022-06-19 17:40:00.0','2022-06-19 17:41:00.0','2022-06-19 17:42:00.0','2022-06-19 17:43:00.0','2022-06-19 17:44:00.0','2022-06-19 17:45:00.0','2022-06-19 17:46:00.0','2022-06-19 17:47:00.0','2022-06-19 17:48:00.0','2022-06-19 17:49:00.0','2022-06-19 17:50:00.0','2022-06-19 17:51:00.0','2022-06-19 17:52:00.0','2022-06-19 17:53:00.0','2022-06-19 17:54:00.0','2022-06-19 17:55:00.0','2022-06-19 17:56:00.0','2022-06-19 17:57:00.0','2022-06-19 17:58:00.0','2022-06-19 17:59:00.0','2022-06-19 18:00:00.0','2022-06-19 18:01:00.0','2022-06-19 18:02:00.0','2022-06-19 18:03:00.0','2022-06-19 18:04:00.0','2022-06-19 18:05:00.0','2022-06-19 18:06:00.0','2022-06-19 18:07:00.0','2022-06-19 18:08:00.0','2022-06-19 18:09:00.0','2022-06-19 18:10:00.0','2022-06-19 18:11:00.0','2022-06-19 18:12:00.0','2022-06-19 18:13:00.0','2022-06-19 18:14:00.0','2022-06-19 18:15:00.0','2022-06-19 18:16:00.0','2022-06-19 18:17:00.0','2022-06-19 18:18:00.0','2022-06-19 18:19:00.0','2022-06-19 18:20:00.0','2022-06-19 18:21:00.0','2022-06-19 18:22:00.0','2022-06-19 18:23:00.0','2022-06-19 18:24:00.0','2022-06-19 18:25:00.0','2022-06-19 18:26:00.0','2022-06-19 18:27:00.0','2022-06-19 18:28:00.0','2022-06-19 18:29:00.0','2022-06-19 18:30:00.0','2022-06-19 18:31:00.0','2022-06-19 18:32:00.0','2022-06-19 18:33:00.0','2022-06-19 18:34:00.0','2022-06-19 18:35:00.0','2022-06-19 18:36:00.0','2022-06-19 18:37:00.0','2022-06-19 18:38:00.0','2022-06-19 18:39:00.0','2022-06-19 18:40:00.0','2022-06-19 18:41:00.0','2022-06-19 18:42:00.0','2022-06-19 18:43:00.0','2022-06-19 18:44:00.0','2022-06-19 18:45:00.0','2022-06-19 18:46:00.0','2022-06-19 18:47:00.0','2022-06-19 18:48:00.0','2022-06-19 18:49:00.0','2022-06-19 18:50:00.0','2022-06-19 18:51:00.0','2022-06-19 18:52:00.0','2022-06-19 18:53:00.0','2022-06-19 18:54:00.0','2022-06-19 18:55:00.0','2022-06-19 18:56:00.0','2022-06-19 18:57:00.0','2022-06-19 18:58:00.0','2022-06-19 18:59:00.0','2022-06-19 19:00:00.0','2022-06-19 19:01:00.0','2022-06-19 19:02:00.0','2022-06-19 19:03:00.0','2022-06-19 19:04:00.0','2022-06-19 19:05:00.0','2022-06-19 19:06:00.0','2022-06-19 19:07:00.0','2022-06-19 19:08:00.0','2022-06-19 19:09:00.0','2022-06-19 19:10:00.0','2022-06-19 19:11:00.0','2022-06-19 19:12:00.0','2022-06-19 19:13:00.0','2022-06-19 19:14:00.0','2022-06-19 19:15:00.0','2022-06-19 19:16:00.0','2022-06-19 19:17:00.0','2022-06-19 19:18:00.0','2022-06-19 19:19:00.0','2022-06-19 19:20:00.0','2022-06-19 19:21:00.0','2022-06-19 19:22:00.0','2022-06-19 19:23:00.0','2022-06-19 19:24:00.0','2022-06-19 19:25:00.0','2022-06-19 19:26:00.0','2022-06-19 19:27:00.0','2022-06-19 19:28:00.0','2022-06-19 19:29:00.0','2022-06-19 19:30:00.0','2022-06-19 19:31:00.0','2022-06-19 19:32:00.0','2022-06-19 19:33:00.0','2022-06-19 19:34:00.0','2022-06-19 19:35:00.0','2022-06-19 19:36:00.0','2022-06-19 19:37:00.0','2022-06-19 19:38:00.0','2022-06-19 19:39:00.0','2022-06-19 19:40:00.0','2022-06-19 19:41:00.0','2022-06-19 19:42:00.0','2022-06-19 19:43:00.0','2022-06-19 19:44:00.0','2022-06-19 19:45:00.0','2022-06-19 19:46:00.0','2022-06-19 19:47:00.0','2022-06-19 19:48:00.0','2022-06-19 19:49:00.0','2022-06-19 19:50:00.0','2022-06-19 19:51:00.0','2022-06-19 19:52:00.0','2022-06-19 19:53:00.0','2022-06-19 19:54:00.0','2022-06-19 19:55:00.0','2022-06-19 19:56:00.0','2022-06-19 19:57:00.0','2022-06-19 19:58:00.0','2022-06-19 19:59:00.0','2022-06-19 20:00:00.0','2022-06-19 20:01:00.0','2022-06-19 20:02:00.0','2022-06-19 20:03:00.0','2022-06-19 20:04:00.0','2022-06-19 20:05:00.0','2022-06-19 20:06:00.0','2022-06-19 20:07:00.0','2022-06-19 20:08:00.0','2022-06-19 20:09:00.0','2022-06-19 20:10:00.0','2022-06-19 20:11:00.0','2022-06-19 20:12:00.0','2022-06-19 20:13:00.0','2022-06-19 20:14:00.0','2022-06-19 20:15:00.0','2022-06-19 20:16:00.0','2022-06-19 20:17:00.0','2022-06-19 20:18:00.0','2022-06-19 20:19:00.0','2022-06-19 20:20:00.0','2022-06-19 20:21:00.0','2022-06-19 20:22:00.0','2022-06-19 20:23:00.0','2022-06-19 20:24:00.0','2022-06-19 20:25:00.0','2022-06-19 20:26:00.0','2022-06-19 20:27:00.0','2022-06-19 20:28:00.0','2022-06-19 20:29:00.0','2022-06-19 20:30:00.0','2022-06-19 20:31:00.0','2022-06-19 20:32:00.0','2022-06-19 20:33:00.0','2022-06-19 20:34:00.0','2022-06-19 20:35:00.0','2022-06-19 20:36:00.0','2022-06-19 20:37:00.0','2022-06-19 20:38:00.0','2022-06-19 20:39:00.0','2022-06-19 20:40:00.0','2022-06-19 20:41:00.0','2022-06-19 20:42:00.0','2022-06-19 20:43:00.0','2022-06-19 20:44:00.0','2022-06-19 20:45:00.0','2022-06-19 20:46:00.0','2022-06-19 20:47:00.0','2022-06-19 20:48:00.0','2022-06-19 20:49:00.0','2022-06-19 20:50:00.0','2022-06-19 20:51:00.0','2022-06-19 20:52:00.0','2022-06-19 20:53:00.0','2022-06-19 20:54:00.0','2022-06-19 20:55:00.0','2022-06-19 20:56:00.0','2022-06-19 20:57:00.0','2022-06-19 20:58:00.0','2022-06-19 20:59:00.0','2022-06-19 21:00:00.0','2022-06-19 21:01:00.0','2022-06-19 21:02:00.0','2022-06-19 21:03:00.0','2022-06-19 21:04:00.0','2022-06-19 21:05:00.0','2022-06-19 21:06:00.0','2022-06-19 21:07:00.0','2022-06-19 21:08:00.0','2022-06-19 21:09:00.0','2022-06-19 21:10:00.0','2022-06-19 21:11:00.0','2022-06-19 21:12:00.0','2022-06-19 21:13:00.0','2022-06-19 21:14:00.0','2022-06-19 21:15:00.0','2022-06-19 21:16:00.0','2022-06-19 21:17:00.0','2022-06-19 21:18:00.0','2022-06-19 21:19:00.0','2022-06-19 21:20:00.0','2022-06-19 21:21:00.0','2022-06-19 21:22:00.0','2022-06-19 21:23:00.0','2022-06-19 21:24:00.0','2022-06-19 21:25:00.0','2022-06-19 21:26:00.0','2022-06-19 21:27:00.0','2022-06-19 21:28:00.0','2022-06-19 21:29:00.0','2022-06-19 21:30:00.0','2022-06-19 21:31:00.0','2022-06-19 21:32:00.0','2022-06-19 21:33:00.0','2022-06-19 21:34:00.0','2022-06-19 21:35:00.0','2022-06-19 21:36:00.0','2022-06-19 21:37:00.0','2022-06-19 21:38:00.0','2022-06-19 21:39:00.0','2022-06-19 21:40:00.0','2022-06-19 21:41:00.0','2022-06-19 21:42:00.0','2022-06-19 21:43:00.0','2022-06-19 21:44:00.0','2022-06-19 21:45:00.0','2022-06-19 21:46:00.0','2022-06-19 21:47:00.0','2022-06-19 21:48:00.0','2022-06-19 21:49:00.0','2022-06-19 21:50:00.0','2022-06-19 21:51:00.0','2022-06-19 21:52:00.0','2022-06-19 21:53:00.0','2022-06-19 21:54:00.0','2022-06-19 21:55:00.0','2022-06-19 21:56:00.0','2022-06-19 21:57:00.0','2022-06-19 21:58:00.0','2022-06-19 21:59:00.0','2022-06-19 22:00:00.0','2022-06-19 22:01:00.0','2022-06-19 22:02:00.0','2022-06-19 22:03:00.0','2022-06-19 22:04:00.0','2022-06-19 22:05:00.0','2022-06-19 22:06:00.0','2022-06-19 22:07:00.0','2022-06-19 22:08:00.0','2022-06-19 22:09:00.0','2022-06-19 22:10:00.0','2022-06-19 22:11:00.0','2022-06-19 22:12:00.0','2022-06-19 22:13:00.0','2022-06-19 22:14:00.0','2022-06-19 22:15:00.0','2022-06-19 22:16:00.0','2022-06-19 22:17:00.0','2022-06-19 22:18:00.0','2022-06-19 22:19:00.0','2022-06-19 22:20:00.0','2022-06-19 22:21:00.0','2022-06-19 22:22:00.0','2022-06-19 22:23:00.0','2022-06-19 22:24:00.0','2022-06-19 22:25:00.0','2022-06-19 22:26:00.0','2022-06-19 22:27:00.0','2022-06-19 22:28:00.0','2022-06-19 22:29:00.0','2022-06-19 22:30:00.0','2022-06-19 22:31:00.0','2022-06-19 22:32:00.0','2022-06-19 22:33:00.0','2022-06-19 22:34:00.0','2022-06-19 22:35:00.0','2022-06-19 22:36:00.0','2022-06-19 22:37:00.0','2022-06-19 22:38:00.0','2022-06-19 22:39:00.0','2022-06-19 22:40:00.0','2022-06-19 22:41:00.0','2022-06-19 22:42:00.0','2022-06-19 22:43:00.0','2022-06-19 22:44:00.0','2022-06-19 22:45:00.0','2022-06-19 22:46:00.0','2022-06-19 22:47:00.0','2022-06-19 22:48:00.0','2022-06-19 22:49:00.0','2022-06-19 22:50:00.0','2022-06-19 22:51:00.0','2022-06-19 22:52:00.0','2022-06-19 22:53:00.0','2022-06-19 22:54:00.0','2022-06-19 22:55:00.0','2022-06-19 22:56:00.0','2022-06-19 22:57:00.0','2022-06-19 22:58:00.0','2022-06-19 22:59:00.0','2022-06-19 23:00:00.0','2022-06-19 23:01:00.0','2022-06-19 23:02:00.0','2022-06-19 23:03:00.0','2022-06-19 23:04:00.0','2022-06-19 23:05:00.0','2022-06-19 23:06:00.0','2022-06-19 23:07:00.0','2022-06-19 23:08:00.0','2022-06-19 23:09:00.0','2022-06-19 23:10:00.0','2022-06-19 23:11:00.0','2022-06-19 23:12:00.0','2022-06-19 23:13:00.0','2022-06-19 23:14:00.0','2022-06-19 23:15:00.0','2022-06-19 23:16:00.0','2022-06-19 23:17:00.0','2022-06-19 23:18:00.0','2022-06-19 23:19:00.0','2022-06-19 23:20:00.0','2022-06-19 23:21:00.0','2022-06-19 23:22:00.0','2022-06-19 23:23:00.')
I can understand saying there is a issue with the quotations, but I cant find it anywhere. It seems as there is a maximum allowed in a vector, if I compile only the first quarter of the entries it works, else wise on R studio it gives an error, and in the terminal it just expects more code so it gives you a + to add more instead of >.
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'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.
I am trying to use a R package called GOSemSim, it requires to import a lot of data into variables with a specific format like this:
data1 = c("one", "two", "three")
data2 = c("A", "B", "C")
When the list of data that I try to import into a variable is longer than 293 then I get the following error message:
argument 293 is empty
THere are no error with the "" or comma, I computed it with linux, it does not matter what data it is.
This is really weird basically, I tried on two computers with no luck. I tried to import it as a CSV file but the R package won't allow it.
Anyone knows why you cannot import more than 293 data?
Update:
Here is the code and my data at the same time, it is a one liner in R which has never been a problem for me!
OQ = c("GO:0000003", "GO:0000070", "GO:0000077", "GO:0000079", "GO:0000082", "GO:0000086", "GO:0000122", "GO:0000212", "GO:0000226", "GO:0000278", "GO:0000279", "GO:0000280", "GO:0000724", "GO:0000725", "GO:0000819", "GO:0000910", "GO:0001932", "GO:0002118", "GO:0002121", "GO:0002165", "GO:0003002", "GO:0003006", "GO:0006022", "GO:0006030", "GO:0006040", "GO:0006139", "GO:0006259", "GO:0006260", "GO:0006261", "GO:0006267", "GO:0006270", "GO:0006275", "GO:0006277", "GO:0006281", "GO:0006302", "GO:0006304", "GO:0006305", "GO:0006306", "GO:0006310", "GO:0006323", "GO:0006325", "GO:0006342", "GO:0006351", "GO:0006355", "GO:0006357", "GO:0006366", "GO:0006464", "GO:0006468", "GO:0006479", "GO:0006725", "GO:0006807", "GO:0006928", "GO:0006950", "GO:0006974", "GO:0006996", "GO:0007010", "GO:0007017", "GO:0007018", "GO:0007049", "GO:0007051", "GO:0007059", "GO:0007062", "GO:0007067", "GO:0007076", "GO:0007088", "GO:0007093", "GO:0007095", "GO:0007098", "GO:0007126", "GO:0007127", "GO:0007131", "GO:0007140", "GO:0007141", "GO:0007143", "GO:0007154", "GO:0007155", "GO:0007156", "GO:0007259", "GO:0007266", "GO:0007275", "GO:0007276", "GO:0007281", "GO:0007282", "GO:0007292", "GO:0007304", "GO:0007307", "GO:0007346", "GO:0007350", "GO:0007365", "GO:0007367", "GO:0007379", "GO:0007389", "GO:0007399", "GO:0007400", "GO:0007417", "GO:0007420", "GO:0007423", "GO:0007444", "GO:0007472", "GO:0007476", "GO:0007552", "GO:0007560", "GO:0008104", "GO:0008213", "GO:0008283", "GO:0008284", "GO:0008315", "GO:0008356", "GO:0009059", "GO:0009611", "GO:0009653", "GO:0009790", "GO:0009791", "GO:0009880", "GO:0009886", "GO:0009887", "GO:0009888", "GO:0009889", "GO:0009890", "GO:0009892", "GO:0009893", "GO:0009896", "GO:0009968", "GO:0009987", "GO:0010032", "GO:0010033", "GO:0010092", "GO:0010389", "GO:0010468", "GO:0010498", "GO:0010556", "GO:0010558", "GO:0010564", "GO:0010604", "GO:0010605", "GO:0010608", "GO:0010629", "GO:0010648", "GO:0010948", "GO:0014016", "GO:0014017", "GO:0014070", "GO:0016043", "GO:0016055", "GO:0016070", "GO:0016310", "GO:0016319", "GO:0016321", "GO:0016441", "GO:0016458", "GO:0016568", "GO:0016569", "GO:0016570", "GO:0016571", "GO:0016572", "GO:0017145", "GO:0018130", "GO:0019219", "GO:0019222", "GO:0019438", "GO:0019827", "GO:0019953", "GO:0022402", "GO:0022403", "GO:0022404", "GO:0022412", "GO:0022414", "GO:0022610", "GO:0023052", "GO:0023057", "GO:0030111", "GO:0030154", "GO:0030178", "GO:0030182", "GO:0030261", "GO:0030422", "GO:0030703", "GO:0030727", "GO:0031023", "GO:0031047", "GO:0031050", "GO:0031056", "GO:0031060", "GO:0031123", "GO:0031145", "GO:0031175", "GO:0031323", "GO:0031324", "GO:0031325", "GO:0031326", "GO:0031327", "GO:0031331", "GO:0031398", "GO:0031399", "GO:0031401", "GO:0031570", "GO:0031572", "GO:0031935", "GO:0032268", "GO:0032270", "GO:0032501", "GO:0032502", "GO:0032504", "GO:0032507", "GO:0032774", "GO:0032776", "GO:0032886", "GO:0033043", "GO:0033044", "GO:0033260", "GO:0033301", "GO:0033554", "GO:0034622", "GO:0034641", "GO:0034645", "GO:0034654", "GO:0034754", "GO:0034968", "GO:0035023", "GO:0035107", "GO:0035114", "GO:0035120", "GO:0035186", "GO:0035194", "GO:0035195", "GO:0035220", "GO:0035282", "GO:0035295", "GO:0035825", "GO:0036211", "GO:0036388", "GO:0040029", "GO:0042060", "GO:0042221", "GO:0042445", "GO:0043009", "GO:0043066", "GO:0043069", "GO:0043161", "GO:0043170", "GO:0043331", "GO:0043412", "GO:0043414", "GO:0043549", "GO:0043631", "GO:0043933", "GO:0044237", "GO:0044249", "GO:0044260", "GO:0044271", "GO:0044419", "GO:0044700", "GO:0044702", "GO:0044703", "GO:0044707", "GO:0044728", "GO:0044763", "GO:0044767", "GO:0044770", "GO:0044771", "GO:0044772", "GO:0044773", "GO:0044774", "GO:0044786", "GO:0044818", "GO:0044839", "GO:0044843", "GO:0044848", "GO:0045132", "GO:0045165", "GO:0045168", "GO:0045185", "GO:0045448", "GO:0045455", "GO:0045787", "GO:0045814", "GO:0045859", "GO:0045892", "GO:0045931", "GO:0045934", "GO:0046331", "GO:0046425", "GO:0046483", "GO:0046580", "GO:0046605", "GO:0046777", "GO:0048070", "GO:0048134", "GO:0048135", "GO:0048285", "GO:0048311", "GO:0048468", "GO:0048477", "GO:0048513", "GO:0048518", "GO:0048519", "GO:0048522", "GO:0048523", "GO:0048563", "GO:0048569", "GO:0048583", "GO:0048585", "GO:0048609", "GO:0048646", "GO:0048666", "GO:0048699", "GO:0048704", "GO:0048705", "GO:0048706", "GO:0048707", "GO:0048731", "GO:0048736", "GO:0048737", "GO:0048754", "GO:0048856", "GO:0048863", "GO:0048865", "GO:0048867", "GO:0048869", "GO:0050789", "GO:0050793", "GO:0050794", "GO:0050896", "GO:0051052", "GO:0051058", "GO:0051128", "GO:0051171", "GO:0051172", "GO:0051225", "GO:0051235", "GO:0051246", "GO:0051247", "GO:0051252", "GO:0051253", "GO:0051276", "GO:0051297", "GO:0051299", "GO:0051301", "GO:0051302", "GO:0051321", "GO:0051325", "GO:0051329", "GO:0051338", "GO:0051351", "GO:0051443", "GO:0051445", "GO:0051641", "GO:0051646", "GO:0051651", "GO:0051704", "GO:0051716", "GO:0051726", "GO:0051783", "GO:0051785", "GO:0060255", "GO:0060429", "GO:0060548", "GO:0060688", "GO:0060966", "GO:0060968", "GO:0060993", "GO:0061138", "GO:0065003", "GO:0065004", "GO:0065007", "GO:0070192", "GO:0070507", "GO:0070887", "GO:0070918", "GO:0071103", "GO:0071359", "GO:0071822", "GO:0071824", "GO:0071840", "GO:0071897", "GO:0071900", "GO:0072028", "GO:0072078", "GO:0072079", "GO:0072088", "GO:0080090", "GO:0090068", "GO:0090304", "GO:0090306", "GO:0098609", "GO:1901071", "GO:1901360", "GO:1901362", "GO:1901576", "GO:1901987", "GO:1901988", "GO:1901990", "GO:1901991", "GO:1902275", "GO:1902299", "GO:1902589", "GO:1902679", "GO:1902749", "GO:1903046", "GO:1903047", "GO:1903308", "GO:1903322", "GO:2000026", "GO:2000112", "GO:2000113", "GO:2001141")
The error message in itself is informative. If one tries to make it reproducible, it's best to work with small subsets. It usually helps to have a dead stare at your data before trying to reproduce the behavior. For example,
OQ = c("GO:0000003", "GO:2001141", )
Notice that there are two elements of this character vector. Or are they?
Error in c("GO:0000003", "GO:2001141", ) : argument 3 is empty
Number 3 is the key. R is expecting three elements. Notice the comma after the second element. Once you remove it, you'll be able to create the QQ variable. Scan your real example. I'm sure there's a , , somewhere.
EDIT
I tried copy/pasting your code into a script in Rstudio and it produced the error you describe. If you scroll right, you'll notice that syntax coloring is not working at around position 5000. I have folded the code so that it fits on screen and it runs fine.
This is how I folded the vector and it worked.
OQ = c("GO:0000003", "GO:0000070", "GO:0000077", "GO:0000079", "GO:0000082", "GO:0000086", "GO:0000122",
"GO:0000212", "GO:0000226", "GO:0000278", "GO:0000279", "GO:0000280", "GO:0000724", "GO:0000725",
"GO:0000819", "GO:0000910", "GO:0001932", "GO:0002118", "GO:0002121", "GO:0002165", "GO:0003002",
"GO:0003006", "GO:0006022", "GO:0006030", "GO:0006040", "GO:0006139", "GO:0006259", "GO:0006260",
"GO:0006261", "GO:0006267", "GO:0006270", "GO:0006275", "GO:0006277", "GO:0006281", "GO:0006302",
"GO:0006304", "GO:0006305", "GO:0006306", "GO:0006310", "GO:0006323", "GO:0006325", "GO:0006342",
"GO:0006351", "GO:0006355", "GO:0006357", "GO:0006366", "GO:0006464", "GO:0006468", "GO:0006479",
"GO:0006725", "GO:0006807", "GO:0006928", "GO:0006950", "GO:0006974", "GO:0006996", "GO:0007010",
"GO:0007017", "GO:0007018", "GO:0007049", "GO:0007051", "GO:0007059", "GO:0007062", "GO:0007067",
"GO:0007076", "GO:0007088", "GO:0007093", "GO:0007095", "GO:0007098", "GO:0007126", "GO:0007127",
"GO:0007131", "GO:0007140", "GO:0007141", "GO:0007143", "GO:0007154", "GO:0007155", "GO:0007156",
"GO:0007259", "GO:0007266", "GO:0007275", "GO:0007276", "GO:0007281", "GO:0007282", "GO:0007292",
"GO:0007304", "GO:0007307", "GO:0007346", "GO:0007350", "GO:0007365", "GO:0007367", "GO:0007379",
"GO:0007389", "GO:0007399", "GO:0007400", "GO:0007417", "GO:0007420", "GO:0007423", "GO:0007444",
"GO:0007472", "GO:0007476", "GO:0007552", "GO:0007560", "GO:0008104", "GO:0008213", "GO:0008283",
"GO:0008284", "GO:0008315", "GO:0008356", "GO:0009059", "GO:0009611", "GO:0009653", "GO:0009790",
"GO:0009791", "GO:0009880", "GO:0009886", "GO:0009887", "GO:0009888", "GO:0009889", "GO:0009890",
"GO:0009892", "GO:0009893", "GO:0009896", "GO:0009968", "GO:0009987", "GO:0010032", "GO:0010033",
"GO:0010092", "GO:0010389", "GO:0010468", "GO:0010498", "GO:0010556", "GO:0010558", "GO:0010564",
"GO:0010604", "GO:0010605", "GO:0010608", "GO:0010629", "GO:0010648", "GO:0010948", "GO:0014016",
"GO:0014017", "GO:0014070", "GO:0016043", "GO:0016055", "GO:0016070", "GO:0016310", "GO:0016319",
"GO:0016321", "GO:0016441", "GO:0016458", "GO:0016568", "GO:0016569", "GO:0016570", "GO:0016571",
"GO:0016572", "GO:0017145", "GO:0018130", "GO:0019219", "GO:0019222", "GO:0019438", "GO:0019827",
"GO:0019953", "GO:0022402", "GO:0022403", "GO:0022404", "GO:0022412", "GO:0022414", "GO:0022610",
"GO:0023052", "GO:0023057", "GO:0030111", "GO:0030154", "GO:0030178", "GO:0030182", "GO:0030261",
"GO:0030422", "GO:0030703", "GO:0030727", "GO:0031023", "GO:0031047", "GO:0031050", "GO:0031056",
"GO:0031060", "GO:0031123", "GO:0031145", "GO:0031175", "GO:0031323", "GO:0031324", "GO:0031325",
"GO:0031326", "GO:0031327", "GO:0031331", "GO:0031398", "GO:0031399", "GO:0031401", "GO:0031570",
"GO:0031572", "GO:0031935", "GO:0032268", "GO:0032270", "GO:0032501", "GO:0032502", "GO:0032504",
"GO:0032507", "GO:0032774", "GO:0032776", "GO:0032886", "GO:0033043", "GO:0033044", "GO:0033260",
"GO:0033301", "GO:0033554", "GO:0034622", "GO:0034641", "GO:0034645", "GO:0034654", "GO:0034754",
"GO:0034968", "GO:0035023", "GO:0035107", "GO:0035114", "GO:0035120", "GO:0035186", "GO:0035194",
"GO:0035195", "GO:0035220", "GO:0035282", "GO:0035295", "GO:0035825", "GO:0036211", "GO:0036388",
"GO:0040029", "GO:0042060", "GO:0042221", "GO:0042445", "GO:0043009", "GO:0043066", "GO:0043069",
"GO:0043161", "GO:0043170", "GO:0043331", "GO:0043412", "GO:0043414", "GO:0043549", "GO:0043631",
"GO:0043933", "GO:0044237", "GO:0044249", "GO:0044260", "GO:0044271", "GO:0044419", "GO:0044700",
"GO:0044702", "GO:0044703", "GO:0044707", "GO:0044728", "GO:0044763", "GO:0044767", "GO:0044770",
"GO:0044771", "GO:0044772", "GO:0044773", "GO:0044774", "GO:0044786", "GO:0044818", "GO:0044839",
"GO:0044843", "GO:0044848", "GO:0045132", "GO:0045165", "GO:0045168", "GO:0045185", "GO:0045448",
"GO:0045455", "GO:0045787", "GO:0045814", "GO:0045859", "GO:0045892", "GO:0045931", "GO:0045934",
"GO:0046331", "GO:0046425", "GO:0046483", "GO:0046580", "GO:0046605", "GO:0046777", "GO:0048070",
"GO:0048134", "GO:0048135", "GO:0048285", "GO:0048311", "GO:0048468", "GO:0048477", "GO:0048513",
"GO:0048518", "GO:0048519", "GO:0048522", "GO:0048523", "GO:0048563", "GO:0048569", "GO:0048583",
"GO:0048585", "GO:0048609", "GO:0048646", "GO:0048666", "GO:0048699", "GO:0048704", "GO:0048705",
"GO:0048706", "GO:0048707", "GO:0048731", "GO:0048736", "GO:0048737", "GO:0048754", "GO:0048856",
"GO:0048863", "GO:0048865", "GO:0048867", "GO:0048869", "GO:0050789", "GO:0050793", "GO:0050794",
"GO:0050896", "GO:0051052", "GO:0051058", "GO:0051128", "GO:0051171", "GO:0051172", "GO:0051225",
"GO:0051235", "GO:0051246", "GO:0051247", "GO:0051252", "GO:0051253", "GO:0051276", "GO:0051297",
"GO:0051299", "GO:0051301", "GO:0051302", "GO:0051321", "GO:0051325", "GO:0051329", "GO:0051338",
"GO:0051351", "GO:0051443", "GO:0051445", "GO:0051641", "GO:0051646", "GO:0051651", "GO:0051704",
"GO:0051716", "GO:0051726", "GO:0051783", "GO:0051785", "GO:0060255", "GO:0060429", "GO:0060548",
"GO:0060688", "GO:0060966", "GO:0060968", "GO:0060993", "GO:0061138", "GO:0065003", "GO:0065004",
"GO:0065007", "GO:0070192", "GO:0070507", "GO:0070887", "GO:0070918", "GO:0071103", "GO:0071359",
"GO:0071822", "GO:0071824", "GO:0071840", "GO:0071897", "GO:0071900", "GO:0072028", "GO:0072078",
"GO:0072079", "GO:0072088", "GO:0080090", "GO:0090068", "GO:0090304", "GO:0090306", "GO:0098609",
"GO:1901071", "GO:1901360", "GO:1901362", "GO:1901576", "GO:1901987", "GO:1901988", "GO:1901990",
"GO:1901991", "GO:1902275", "GO:1902299", "GO:1902589", "GO:1902679", "GO:1902749", "GO:1903046",
"GO:1903047", "GO:1903308", "GO:1903322", "GO:2000026", "GO:2000112", "GO:2000113", "GO:2001141")
I have an ARC/INFO generate file whose contents look like:
3594 -124.049541 44.429077
-123.381222 44.530192
-123.479913 44.625517
-123.578917 44.720704
-123.678234 44.815755
-123.777866 44.910669
-123.946044 44.885032
-124.114074 44.858987
-124.281949 44.832529
-124.449663 44.805654
-124.516511 44.684660
-124.583091 44.563597
-124.649404 44.442465
-124.715451 44.321261
-124.615376 44.227772
-124.515601 44.134147
-124.416125 44.040385
-124.316948 43.946486
-124.151513 43.973082
-123.985926 43.999247
-123.820193 44.024987
-123.654322 44.050307
-123.586447 44.170362
-123.518307 44.290360
-123.449899 44.410303
-123.381222 44.530192
END
3595 -123.103772 45.009223
-122.427717 45.101578
-122.525757 45.198252
-122.624122 45.294789
-122.722814 45.391191
-122.821833 45.487459
-122.992014 45.464007
-123.162072 45.440175
-123.332002 45.415959
-123.501798 45.391355
-123.571234 45.271264
-123.640389 45.151121
-123.709266 45.030923
-123.777866 44.910669
-123.678234 44.815755
-123.578917 44.720704
-123.479913 44.625517
-123.381222 44.530192
-123.213811 44.554460
-123.046278 44.578334
-122.878629 44.601816
-122.710869 44.624913
-122.640504 44.744148
-122.569859 44.863337
-122.498931 44.982480
-122.427717 45.101578
END
3676 -122.989567 44.147495
-122.323040 44.238368
-122.419523 44.335217
-122.516322 44.431923
-122.613437 44.528488
-122.710869 44.624913
-122.878629 44.601816
-123.046278 44.578334
-123.213811 44.554460
-123.381222 44.530192
-123.449899 44.410303
-123.518307 44.290360
-123.586447 44.170362
-123.654322 44.050307
-123.556277 43.955264
-123.458534 43.860080
-123.361093 43.764751
-123.263953 43.669279
-123.098838 43.693189
-122.933613 43.716694
-122.768285 43.739802
-122.602857 43.762515
-122.533309 43.881546
-122.463492 44.000532
-122.393403 44.119472
-122.323040 44.238368
END
END
My strategy is to read in the file generating a list of latitude-longitude points and beginning a new unique group id every time I encounter an END. I'll then plot using ggplot" andgeom_polygon".
Alas, I'm not sure how to efficiently accomplish the reading of the file.
Any thoughts?
Read the spatial task view on CRAN and then use readOGR from the rgdal package to read into an sp class object. You'll need a GDAL/OGR install with ARCGEN format support, which despite being listed as 'compiled by default' Link I don't have on my system.
Failing that, open the file as a connection, read each line, build a Polygon, then Polygons and SpatialPolygons.
Here's a fairly sub-optimal but working function:
readUng <- function(f){
require(sp)
stream = file(f,"r")
first = readLines(stream,1)
bits = strsplit(first," ")[[1]]
polys = list();ids=NULL
while(TRUE){
id=bits[1] # label pt = bits[2],bits[3]
ids=c(ids,id)
coords=NULL
while(TRUE){
xy=readLines(stream,1)
if(xy=="END"){
break
}
coords=rbind(coords,strsplit(xy," ")[[1]])
}
polys[[length(polys)+1]] = Polygons(list(Polygon(matrix(as.numeric(coords[,2:3]),ncol=2))),ID=id)
lines = readLines(stream,1)
if(lines == "END"){
break
}
bits = strsplit(lines," ")[[1]]
}
return(SpatialPolygons(polys))
}
Now its a proper spatial data object, you can also give it a coordinate system (looks like lat-long to me, so epsg:4326, but only you know). Now you could modify all this to produce whatever ggplot wants, but if its spatial data then you should keep it as a spatial data class and ggplot should be made capable of dealing with such.