Hello,
I am working with the ETAS pachkage. I created a ppx file from my data
(attached file example.csv) and now i am trying to run the etas model but an
error occurs. Could you help me please on how to fix this problem?I am
copy-pasting the procedure i followed:> library(spatstat)Loading required package: mgcvLoading required package:
nlmeThis is mgcv 1.7-26. For overview type
'help("mgcv-package")'.Loading required package: deldirdeldir
0.1-1Loading required package: abindLoading required package: tensorLoading
required package: polyclippolyclip 1.2-0
spatstat 1.35-0 (nickname: ‘Multiple Personality’) For an introduction to
spatstat, type ‘beginner’> library(SAPP)> library(ETAS)>
data<-read.csv("example.csv")> edit(data) X No longitude
latitude magnitude time depth year month day1 1 1 23.33 35.47
6.2 0.0000 42.46 2013 10 122 2 2 23.26 35.29 3.6
0.0035 32.07 2013 10 123 3 3 23.34 35.41 2.5 0.0056
53.99 2013 10 124 4 4 23.24 35.19 2.9 0.0069 28.85 2013
10 125 5 5 23.33 35.34 3.2 0.0069 30.38 2013 10 126
6 6 23.28 35.32 2.7 0.0083 36.29 2013 10 127 7 7
23.34 35.48 2.4 0.0090 38.38 2013 10 128 8 8 23.39
35.41 2.3 0.0104 44.65 2013 10 129 9 9 23.75 35.13
1.8 0.0111 102.00 2013 10 1210 10 10 23.30 35.39 2.4
0.0118 28.33 2013 10 1211 11 11 23.29 35.38 2.8 0.0125
35.81 2013 10 1212 12 12 23.32 35.43 2.2 0.0146 35.61 2013
10 1213 13 13 23.34 35.44 2.7 0.0153 29.87 2013 10 1214
14 14 23.33 35.35 2.8 0.0160 43.36 2013 10 1215 15 15
23.49 35.28 2.1 0.0188 46.86 2013 10 1216 16 16 23.33
35.39 2.7 0.0201 28.59 2013 10 1217 17 17 23.33 35.41
2.3 0.0229 32.36 2013 10 1218 18 18 23.31 35.31 1.8
0.0250 37.69 2013 10 1219 19 19 23.29 35.35 1.7 0.0278
28.53 2013 10 1220 20 20 23.27 35.30 2.0 0.0292 33.51 2013
10 1221 21 21 23.34 35.37 2.8 0.0306 34.84 2013 10 1222
22 22 23.38 35.47 1.4 0.0319 34.44 2013 10 1223 23 23
23.66 35.35 1.6 0.0333 51.00 2013 10 1224 24 24 23.36
35.47 2.4 0.0347 28.95 2013 10 1225 25 25 23.37 35.37
2.1 0.0354 30.17 2013 10 1226 26 26 23.71 35.34 2.1
0.0361 33.00 2013 10 1227 27 27 23.50 35.21 1.9 0.0375
50.66 2013 10 1228 28 28 23.36 35.44 3.6 0.0375 36.44 2013
10 1229 29 29 23.43 35.51 1.4 0.0417 39.86 2013 10 1230
30 30 23.53 35.14 1.9 0.0424 84.00 2013 10 1231 31 31
23.32 35.42 2.7 0.0438 34.16 2013 10 1232 32 32 23.37
35.41 2.6 0.0444 50.39 2013 10 1233 33 33 23.31 35.32
2.6 0.0451 27.26 2013 10 1234 34 34 23.26 35.34 1.7
0.0486 37.03 2013 10 1235 35 35 23.30 35.32 1.8 0.0500
27.80 2013 10 1236 36 36 23.26 35.38 1.8 0.0500 47.43 2013
10 1237 37 37 23.50 35.39 1.4 0.0521 40.79 2013 10 1238
38 38 23.39 35.38 1.5 0.0542 65.69 2013 10 1239 39 39
23.49 35.36 2.0 0.0563 40.30 2013 10 1240 40 40 23.23
35.31 1.4 0.0590 58.46 2013 10 1241 41 41 23.26 35.27
1.8 0.0597 56.33 2013 10 1242 42 42 23.30 35.37 2.1
0.0604 26.81 2013 10 1243 43 43 23.30 35.40 1.5 0.0618
38.99 2013 10 1244 44 44 23.32 35.39 1.7 0.0632 30.62 2013
10 1245 45 45 23.27 35.30 1.8 0.0646 48.67 2013 10 1246
46 46 22.98 35.52 2.4 0.0653 16.24 2013 10 1247 47 47
23.34 35.33 1.5 0.0660 45.26 2013 10 1248 48 48 23.32
35.29 1.9 0.0667 42.36 2013 10 1249 49 49 23.29 35.33
1.6 0.0681 38.23 2013 10 1250 50 50 23.36 35.41 1.6
0.0701 33.13 2013 10 1251 51 51 23.38 35.36 1.5 0.0715
47.19 2013 10 1252 52 52 23.27 35.31 2.0 0.0722 27.08 2013
10 1253 53 53 23.33 35.40 1.7 0.0729 42.67 2013 10 1254
54 54 23.25 35.42 2.9 0.0736 40.36 2013 10 1255 55 55
23.32 35.34 1.8 0.0757 20.38 2013 10 1256 56 56 23.31
35.41 1.6 0.0764 32.29 2013 10 1257 57 57 23.35 35.34
1.7 0.0778 53.87 2013 10 1258 58 58 23.28 35.33 1.6
0.0792 41.10 2013 10 1259 59 59 23.46 35.38 1.8 0.0792
39.68 2013 10 1260 60 60 23.34 35.41 1.4 0.0826 29.10 2013
10 1261 61 61 23.29 35.34 1.7 0.0833 23.91 2013 10 1262
62 62 23.56 35.37 2.4 0.0833 41.32 2013 10 1263 63 63
23.31 35.37 2.5 0.0847 41.17 2013 10 1264 64 64 23.53
35.36 1.6 0.0854 53.99 2013 10 1265 65 65 23.36 35.44
2.0 0.0854 32.22 2013 10 1266 66 66 23.39 35.40 1.3
0.0882 41.86 2013 10 1267 67 67 23.53 35.29 1.7 0.0896
63.00 2013 10 1268 68 68 23.36 35.32 1.7 0.0917 54.27 2013
10 1269 69 69 23.30 35.40 3.1 0.0944 34.64 2013 10 1270
70 70 23.41 35.34 1.2 0.0993 61.03 2013 10 1271 71 71
23.32 35.40 1.4 0.0993 45.82 2013 10 1272 72 72 23.31
35.29 1.1 0.1007 48.53 2013 10 1273 73 73 23.05 34.80
1.8 0.1056 43.31 2013 10 1274 74 74 23.29 35.21 2.0
0.1083 50.83 2013 10 1275 75 75 23.26 35.35 1.6 0.1097
46.08 2013 10 1276 76 76 23.30 35.42 1.5 0.1139 28.76 2013
10 1277 77 77 23.27 35.34 1.6 0.1208 43.02 2013 10 1278
78 78 23.25 35.31 1.7 0.1229 40.15 2013 10 1279 79 79
23.29 35.38 1.9 0.1271 32.71 2013 10 1280 80 80 23.30
35.37 1.9 0.1278 6.12 2013 10 1281 81 81 23.37 35.39
1.5 0.1306 36.27 2013 10 1282 82 82 23.29 35.30 1.8
0.1347 34.38 2013 10 1283 83 83 23.58 35.44 1.9 0.1389
44.60 2013 10 1284 84 84 23.34 35.44 1.1 0.1403 36.32 2013
10 1285 85 85 23.28 35.32 0.8 0.1431 46.46 2013 10 1286
86 86 23.76 35.49 1.2 0.1438 48.98 2013 10 1287 87 87
23.33 35.34 1.5 0.1486 31.62 2013 10 1288 88 88 23.31
35.41 2.3 0.1535 33.68 2013 10 1289 89 89 23.64 35.31
1.9 0.1646 24.00 2013 10 1290 90 90 23.59 35.30 0.8
0.1674 60.00 2013 10 1291 91 91 23.35 35.44 2.2 0.1694
33.20 2013 10 1292 92 92 23.34 35.33 1.9 0.1736 46.13 2013
10 1293 93 93 23.34 35.33 1.7 0.1778 53.84 2013 10 1294
94 94 23.35 35.28 1.4 0.1813 63.10 2013 10 1295 95 95
23.30 35.18 1.6 0.1875 51.83 2013 10 1296 96 96 23.37
35.45 1.7 0.1931 3.50 2013 10 1297 97 97 23.37 35.45
1.6 0.1958 28.36 2013 10 1298 98 98 23.30 35.41 1.4
0.2007 33.89 2013 10 1299 99 99 23.34 35.29 1.8 0.2035
50.60 2013 10 12100 100 100 23.32 35.44 1.9 0.2049 22.07 2013
10 12101 101 101 23.26 35.17 1.3 0.2063 51.05 2013 10 12102
102 102 23.58 35.37 1.7 0.2097 37.83 2013 10 12103 103 103
23.32 35.44 1.5 0.2111 30.49 2013 10 12104 104 104 23.37
35.46 1.9 0.2153 39.00 2013 10 12105 105 105 23.55 35.27
1.7 0.2208 57.00 2013 10 12106 106 106 23.57 35.37 1.7
0.2215 40.90 2013 10 12107 107 107 23.30 35.26 1.6 0.2326
48.05 2013 10 12108 108 108 23.27 35.39 2.3 0.2444 29.03 2013
10 12109 109 109 23.36 35.44 1.3 0.2493 39.60 2013 10 12110
110 110 23.40 35.33 1.6 0.2569 54.49 2013 10 12111 111 111
23.32 35.45 1.8 0.2583 64.86 2013 10 12112 112 112 23.31
35.44 1.7 0.2604 32.28 2013 10 12113 113 113 23.32 35.41
3.6 0.2674 36.26 2013 10 12114 114 114 22.94 35.32 3.2
0.2694 16.32 2013 10 12115 115 115 23.54 35.58 2.0 0.2701
65.79 2013 10 12116 116 116 23.31 35.50 2.0 0.2708 59.39 2013
10 12117 117 117 23.32 35.36 1.3 0.2715 40.23 2013 10 12118
118 118 23.34 35.45 1.8 0.2729 32.22 2013 10 12119 119 119
23.30 35.39 1.3 0.2743 7.88 2013 10 12120 120 120 23.50
35.82 1.4 0.2750 22.18 2013 10 12121 121 121 23.75 35.49
1.8 0.2799 39.00 2013 10 12122 122 122 23.43 35.45 0.6
0.2931 56.29 2013 10 12123 123 123 23.24 35.26 1.8 0.2951
36.13 2013 10 12124 124 124 23.30 35.35 1.8 0.3000 27.80 2013
10 12125 125 125 23.31 35.21 1.8 0.3056 53.34 2013 10 12126
126 126 23.59 35.27 1.3 0.3063 60.00 2013 10 12127 127 127
23.28 35.39 2.7 0.3215 39.96 2013 10 12128 128 128 23.33
35.35 1.5 0.3264 6.91 2013 10 12129 129 129 23.32 35.44
1.1 0.3264 26.95 2013 10 12130 130 130 23.23 35.25 2.1
0.3313 50.47 2013 10 12131 131 131 23.46 35.23 1.6 0.3333
61.48 2013 10 12132 132 132 23.67 35.20 1.6 0.3340 18.00 2013
10 12133 133 133 23.23 35.25 1.7 0.3403 53.93 2013 10 12134
134 134 23.29 35.40 2.1 0.3549 28.89 2013 10 12135 135 135
23.23 35.29 1.9 0.3583 25.75 2013 10 12136 136 136 23.31
35.32 2.1 0.3590 35.19 2013 10 12137 137 137 23.32 35.40
2.4 0.3660 34.13 2013 10 12138 138 138 22.89 35.31 1.8
0.3778 55.84 2013 10 12139 139 139 23.34 35.42 2.2 0.3785
25.90 2013 10 12140 140 140 23.28 35.42 1.8 0.3792 36.98 2013
10 12141 141 141 23.30 35.05 1.0 0.3833 39.70 2013 10 12142
142 142 23.32 35.42 1.5 0.3958 32.99 2013 10 12143 143 143
23.38 35.44 1.0 0.3986 41.17 2013 10 12144 144 144 23.31
35.41 2.0 0.4021 33.61 2013 10 12145 145 145 23.43 34.90
2.2 0.4028 82.38 2013 10 12146 146 146 23.27 35.37 1.4
0.4069 45.42 2013 10 12147 147 147 23.57 35.35 2.0 0.4090
42.00 2013 10 12148 148 148 23.29 35.15 1.7 0.4132 50.67 2013
10 12149 149 149 23.30 35.42 1.3 0.4139 28.51 2013 10 12150
150 150 23.33 35.41 2.8 0.4188 32.46 2013 10 12151 151 151
23.33 35.32 1.5 0.4403 51.67 2013 10 12152 152 152 23.34
35.41 1.7 0.4417 21.84 2013 10 12153 153 153 23.30 35.34
3.4 0.4514 29.21 2013 10 13154 154 154 23.25 35.32 1.1
0.4583 15.24 2013 10 13155 155 155 23.36 35.48 1.3 0.4674
43.65 2013 10 13156 156 156 23.27 35.36 1.4 0.4715 30.87 2013
10 13157 157 157 23.38 35.47 1.7 0.4715 42.01 2013 10 13158
158 158 23.54 35.35 0.9 0.4806 51.13 2013 10 13159 159 159
23.26 35.28 1.1 0.4826 39.46 2013 10 13160 160 160 23.27
35.38 0.9 0.4826 9.79 2013 10 13161 161 161 23.30 35.30
1.2 0.5257 49.46 2013 10 13162 162 162 23.37 35.44 1.1
0.5354 40.10 2013 10 13163 163 163 23.31 35.36 1.7 0.5389
33.93 2013 10 13164 164 164 23.33 35.38 1.4 0.5479 20.66 2013
10 13165 165 165 23.39 35.35 1.2 0.5604 50.43 2013 10 13166
166 166 23.26 35.29 1.1 0.5611 44.87 2013 10 13167 167 167
23.27 35.29 1.6 0.5618 41.08 2013 10 13168 168 168 23.35
35.46 1.2 0.5632 40.83 2013 10 13169 169 169 23.30 35.41
2.1 0.5771 39.61 2013 10 13170 170 170 23.35 35.40 1.5
0.5833 6.68 2013 10 13171 171 171 23.33 35.24 1.5 0.5861
59.79 2013 10 13172 172 172 23.26 35.25 2.2 0.5875 44.06 2013
10 13173 173 173 23.49 35.35 1.1 0.5882 63.00 2013 10 13174
174 174 23.31 35.41 3.0 0.6056 39.42 2013 10 13175 175 175
23.21 35.34 1.6 0.6139 21.59 2013 10 13176 176 176 23.73
35.51 0.6 0.6264 39.00 2013 10 13177 177 177 23.24 35.29
1.2 0.6333 48.37 2013 10 13178 178 178 23.60 35.21 1.5
0.6347 15.00 2013 10 13179 179 179 23.63 35.19 0.6 0.6438
9.00 2013 10 13180 180 180 23.34 35.06 1.2 0.6465 5.45 2013
10 13181 181 181 23.66 35.19 2.3 0.6514 15.00 2013 10 13182
182 182 23.31 35.38 1.6 0.6542 39.06 2013 10 13183 183 183
23.32 35.11 1.2 0.6556 13.85 2013 10 13184 184 184 23.33
35.36 1.2 0.6597 46.37 2013 10 13185 185 185 23.28 35.28
1.4 0.6604 37.61 2013 10 13186 186 186 23.65 35.20 0.8
0.6847 9.47 2013 10 13187 187 187 23.38 35.47 2.0 0.7139
26.39 2013 10 13188 188 188 23.27 35.35 1.8 0.7278 23.26 2013
10 13189 189 189 23.26 35.37 1.5 0.7403 35.55 2013 10 13190
190 190 23.31 35.39 1.3 0.7431 43.72 2013 10 13191 191 191
23.31 35.37 1.5 0.7451 41.36 2013 10 13192 192 192 23.32
35.43 1.4 0.7715 3.00 2013 10 13193 193 193 23.02 34.97
1.6 0.7924 17.81 2013 10 13194 194 194 23.96 35.26 0.6
0.8069 95.85 2013 10 13195 195 195 23.46 35.33 1.1 0.8229
57.14 2013 10 13196 196 196 23.62 35.32 1.4 0.8271 72.00 2013
10 13197 197 197 23.43 35.69 1.5 0.8549 3.04 2013 10 13198
198 198 23.40 35.47 1.6 0.8833 39.64 2013 10 13199 199 199
23.29 35.39 1.7 0.9035 21.86 2013 10 13200 200 200 23.48
35.23 1.2 0.9069 58.41 2013 10 13201 201 201 23.24 35.30
1.8 0.9174 16.41 2013 10 13202 202 202 23.38 35.39 1.4
0.9313 31.25 2013 10 13203 203 203 23.32 35.42 1.7 0.9493
19.93 2013 10 13204 204 204 23.30 35.30 2.0 0.9507 42.24 2013
10 13205 205 205 23.38 35.52 1.1 0.9611 37.07 2013 10 13206
206 206 23.35 35.41 1.6 0.9889 8.43 2013 10 13207 207 207
23.45 35.39 1.2 0.9944 50.62 2013 10 13208 208 208 23.33
35.09 3.2 1.0007 8.11 2013 10 13209 209 209 23.33 35.32
1.8 1.0410 22.20 2013 10 13210 210 210 23.24 35.30 1.5
1.0660 2.82 2013 10 13211 211 211 23.46 35.20 2.0 1.0681
53.39 2013 10 13212 212 212 23.28 35.30 2.0 1.0688 33.17 2013
10 13213 213 213 23.32 35.36 1.5 1.0708 5.87 2013 10 13214
214 214 23.33 35.13 1.6 1.0764 66.68 2013 10 13215 215 215
23.23 35.23 1.6 1.1299 35.20 2013 10 13216 216 216 23.48
35.31 1.8 1.1326 57.82 2013 10 13217 217 217 23.23 35.20
1.6 1.1507 65.23 2013 10 13218 218 218 23.45 35.45 1.9
1.1597 30.43 2013 10 13219 219 219 23.29 35.39 3.7 1.1889
38.55 2013 10 13220 220 220 23.36 35.37 2.0 1.1938 35.97 2013
10 13221 221 221 23.29 35.37 2.0 1.2167 42.33 2013 10 13222
222 222 23.20 35.42 1.6 1.2528 81.14 2013 10 13223 223 223
23.34 35.41 2.3 1.2549 42.24 2013 10 13224 224 224 23.36
35.30 1.6 1.2576 47.52 2013 10 13225 225 225 23.55 35.27
1.9 1.2701 54.00 2013 10 13226 226 226 23.48 35.19 1.3
1.2868 50.17 2013 10 13227 227 227 23.35 35.45 1.1 1.3618
22.56 2013 10 13228 228 228 23.29 35.33 1.7 1.3667 43.86 2013
10 13229 229 229 23.60 35.16 1.9 1.3708 16.20 2013 10 13230
230 230 23.33 35.42 1.6 1.3847 35.67 2013 10 13231 231 231
22.98 34.91 2.1 1.3979 17.89 2013 10 13232 232 232 23.67
35.17 2.5 1.4493 18.00 2013 10 13233 233 233 23.29 35.36
1.3 1.4688 23.43 2013 10 14234 234 234 23.30 35.36 2.1
1.4743 31.47 2013 10 14235 235 235 23.32 35.35 1.5 1.4986
26.11 2013 10 14236 236 236 23.25 35.20 0.9 1.6118 46.03 2013
10 14237 237 237 23.35 35.43 1.0 1.6944 31.98 2013 10 14238
238 238 23.32 35.42 2.7 1.7500 40.27 2013 10 14239 239 239
23.37 35.40 1.6 1.7799 6.70 2013 10 14240 240 240 23.45
35.41 1.7 1.7840 45.04 2013 10 14241 241 241 23.31 35.35
2.1 1.7951 33.00 2013 10 14242 242 242 23.33 35.37 1.4
1.8007 4.41 2013 10 14243 243 243 23.30 35.37 1.8 1.8424
8.78 2013 10 14244 244 244 23.29 35.36 1.8 1.8917 55.49 2013
10 14245 245 245 23.44 35.38 1.5 1.9125 16.14 2013 10 14246
246 246 23.35 35.40 1.5 1.9542 43.38 2013 10 14247 247 247
23.39 35.51 1.8 1.9903 13.43 2013 10 14248 248 248 23.43
35.20 1.6 2.0000 63.00 2013 10 14249 249 249 23.34 35.48
1.1 2.0021 30.71 2013 10 14250 250 250 23.28 34.86 1.9
2.0049 60.47 2013 10 14251 251 251 23.32 35.44 1.9 2.0069
27.19 2013 10 14252 252 252 23.28 35.29 1.6 2.0903 41.85 2013
10 14253 253 253 23.35 35.37 1.6 2.1604 43.74 2013 10 14254
254 254 23.23 35.34 1.5 2.2438 30.80 2013 10 14255 255 255
23.50 35.75 1.3 2.2674 16.13 2013 10 14256 256 256 23.41
35.53 1.8 2.2729 24.27 2013 10 14257 257 257 23.32 35.38
2.0 2.2854 33.40 2013 10 14258 258 258 23.41 35.13 1.7
2.2875 56.59 2013 10 14259 259 259 23.20 35.16 1.7 2.2917
60.16 2013 10 14260 260 260 23.70 35.25 1.2 2.3167 4.07 2013
10 14261 261 261 23.32 35.34 1.5 2.3632 26.49 2013 10 14262
262 262 23.34 35.33 1.3 2.4160 29.00 2013 10 14263 263 263
23.57 34.94 1.2 2.4417 25.86 2013 10 14264 264 264 23.23
35.33 1.8 2.4507 58.43 2013 10 15265 265 265 23.59 35.18
1.3 2.4646 15.38 2013 10 15266 266 266 23.45 35.88 1.1
2.5215 17.86 2013 10 15267 267 267 23.63 35.15 1.5 2.5306
16.39 2013 10 15268 268 268 23.54 35.31 1.3 2.5375 50.40 2013
10 15269 269 269 23.28 35.18 1.4 2.5597 56.07 2013 10 15270
270 270 23.66 35.26 0.8 2.5799 9.00 2013 10 15271 271 271
23.33 35.41 1.3 2.5875 3.21 2013 10 15272 272 272 23.46
35.90 1.1 2.6153 14.40 2013 10 15273 273 273 23.31 35.31
2.2 2.6493 26.93 2013 10 15274 274 274 23.31 35.38 2.1
2.6625 39.78 2013 10 15275 275 275 22.93 35.22 2.1 2.6889
55.32 2013 10 15276 276 276 23.67 35.17 0.8 2.7160 11.80 2013
10 15277 277 277 23.27 35.40 1.4 2.7472 2.88 2013 10 15278
278 278 23.59 35.33 1.1 2.7910 18.00 2013 10 15279 279 279
23.29 35.30 1.9 2.8132 23.54 2013 10 15280 280 280 23.48
35.25 1.2 2.8167 63.00 2013 10 15281 281 281 22.90 34.98
2.3 2.8549 29.77 2013 10 15282 282 282 23.37 35.31 1.1
2.8701 57.51 2013 10 15283 283 283 23.24 35.33 1.5 2.8826
40.17 2013 10 15284 284 284 23.36 35.39 0.8 2.9132 22.74 2013
10 15285 285 285 23.50 34.93 0.9 2.9167 32.71 2013 10 15286
286 286 23.30 35.40 0.6 2.9646 41.17 2013 10 15287 287 287
23.35 35.36 1.6 2.9681 4.62 2013 10 15288 288 288 22.55
35.39 1.5 2.9833 78.60 2013 10 15289 289 289 23.63 35.17
1.6 3.0194 18.00 2013 10 15290 290 290 23.29 35.39 1.7
3.0764 42.79 2013 10 15291 291 291 23.32 35.38 1.8 3.0875
48.76 2013 10 15292 292 292 23.25 35.28 1.8 3.1083 32.89 2013
10 15293 293 293 23.63 35.15 2.1 3.1278 39.00 2013 10 15294
294 294 23.58 35.29 1.3 3.1618 24.00 2013 10 15295 295 295
23.70 34.74 1.4 3.2271 10.37 2013 10 15296 296 296 23.24
35.30 1.7 3.2486 31.33 2013 10 15297 297 297 23.24 35.33
1.6 3.2535 26.83 2013 10 15298 298 298 23.31 35.39 1.4
3.3785 38.98 2013 10 15299 299 299 23.27 35.13 2.0 3.4437
2.98 2013 10 15300 300 300 23.42 35.44 2.6 3.5264 39.99 2013
10 16301 301 301 23.44 35.41 1.8 3.5944 32.93 2013 10 16302
302 302 23.39 35.45 2.6 3.6208 40.92 2013 10 16303 303 303
23.37 35.45 2.3 3.6243 48.51 2013 10 16304 304 304 23.11
35.31 2.1 3.6674 50.01 2013 10 16305 305 305 22.92 35.26
1.8 3.7569 52.81 2013 10 16306 306 306 23.29 35.34 2.0
3.7833 29.02 2013 10 16307 307 307 23.26 35.19 2.0 3.8472
34.99 2013 10 16308 308 308 23.48 35.23 2.1 3.9285 46.84 2013
10 16309 309 309 23.26 35.34 2.1 4.0194 81.16 2013 10 16310
310 310 23.40 35.42 3.4 4.0250 33.58 2013 10 16311 311 311
22.93 35.23 3.0 4.0271 33.72 2013 10 16312 312 312 23.19
35.28 2.0 4.0597 78.68 2013 10 16313 313 313 22.89 35.17
1.8 4.1194 16.32 2013 10 16314 314 314 23.22 35.31 1.4
4.1250 47.65 2013 10 16315 315 315 23.68 34.96 1.6 4.2194
56.19 2013 10 16316 316 316 23.31 35.34 2.0 4.2944 45.43 2013
10 16317 317 317 23.36 35.40 2.2 4.3729 39.19 2013 10 16318
318 318 23.53 35.27 1.6 4.3819 48.88 2013 10 16319 319 319
23.31 35.34 1.8 4.5861 29.15 2013 10 17320 320 320 23.23
35.24 2.1 4.6597 16.07 2013 10 17321 321 321 23.61 35.35
1.1 4.8458 33.14 2013 10 17322 322 322 23.65 35.08 2.3
5.2076 59.37 2013 10 17323 323 323 23.55 35.31 0.3 5.4132
24.00 2013 10 17324 324 324 22.72 35.66 2.0 5.4299 71.60 2013
10 17325 325 325 23.31 35.18 2.1 5.5319 53.16 2013 10 18326
326 326 23.45 35.42 1.1 5.7410 46.78 2013 10 18327 327 327
23.59 35.20 2.7 5.8063 16.48 2013 10 18328 328 328 23.27
35.21 2.3 5.9382 43.89 2013 10 18329 329 329 23.31 35.39
1.7 6.0375 12.37 2013 10 18330 330 330 23.33 35.36 1.6
6.0694 15.37 2013 10 18331 331 331 23.44 35.41 2.2 6.0937
42.90 2013 10 18332 332 332 23.30 35.40 1.9 6.2201 10.79 2013
10 18333 333 333 23.37 35.47 2.8 6.2660 36.45 2013 10 18334
334 334 23.31 35.41 1.7 6.2799 44.79 2013 10 18335 335 335
23.69 34.96 1.5 6.4208 22.95 2013 10 18336 336 336 23.70
34.94 1.2 6.4563 16.34 2013 10 19337 337 337 23.29 35.29
1.7 6.4736 40.11 2013 10 19338 338 338 23.28 35.12 1.8
6.5472 1.57 2013 10 19339 339 339 23.53 35.30 3.3 6.5653
47.30 2013 10 19340 340 340 23.27 35.37 2.8 6.6028 40.20 2013
10 19341 341 341 23.34 35.37 2.6 6.6028 26.26 2013 10 19342
342 342 23.38 35.43 3.0 6.6458 30.56 2013 10 19343 343 343
23.28 35.32 1.9 6.6556 21.72 2013 10 19344 344 344 23.32
35.35 1.9 6.6618 44.80 2013 10 19345 345 345 23.42 35.40
1.7 6.7514 3.05 2013 10 19346 346 346 23.26 35.11 1.3
6.7778 38.23 2013 10 19347 347 347 23.80 35.58 1.0 6.8528
33.00 2013 10 19348 348 348 23.42 35.53 1.7 6.9069 69.66 2013
10 19349 349 349 23.45 35.49 1.7 6.9910 9.17 2013 10 19350
350 350 23.70 35.44 1.1 7.0229 39.00 2013 10 19351 351 351
23.39 35.47 2.2 7.0590 38.56 2013 10 19352 352 352 23.41
35.58 2.2 7.1861 16.22 2013 10 19353 353 353 23.32 35.41
1.1 7.2715 36.76 2013 10 19354 354 354 23.31 35.41 1.5
7.3000 8.72 2013 10 19355 355 355 23.30 35.36 2.0 7.4507
35.99 2013 10 19356 356 356 23.99 35.09 1.6 7.6514 7.18 2013
10 20357 357 357 23.32 35.39 3.3 7.6806 32.25 2013 10 20358
358 358 23.94 35.06 1.7 7.7146 4.88 2013 10 20359 359 359
23.31 35.42 2.7 7.7569 26.03 2013 10 20360 360 360 23.26
35.29 2.3 7.8104 31.92 2013 10 20361 361 361 23.24 35.14
1.6 7.8243 4.45 2013 10 20362 362 362 23.52 35.34 1.9
7.9083 42.00 2013 10 20363 363 363 23.69 34.92 1.5 8.0493
40.44 2013 10 20364 364 364 22.54 35.77 2.4 8.1563 56.30 2013
10 20365 365 365 23.59 35.22 2.0 8.1757 15.00 2013 10 20366
366 366 23.61 35.33 0.8 8.1944 31.29 2013 10 20367 367 367
23.29 35.33 2.2 8.4049 45.36 2013 10 20368 368 368 23.31
35.39 2.4 8.5583 35.00 2013 10 20369 369 369 23.32 35.51
2.3 8.6458 39.12 2013 10 21370 370 370 23.18 35.98 1.2
9.1264 13.78 2013 10 21371 371 371 23.51 35.38 1.9 9.1917
39.44 2013 10 21372 372 372 23.33 35.12 1.7 9.2111 6.38 2013
10 21373 373 373 23.62 35.19 0.8 9.2146 18.00 2013 10 21374
374 374 23.34 35.37 1.2 9.2222 11.43 2013 10 21375 375 375
23.33 35.39 1.4 9.2313 48.93 2013 10 21376 376 376 23.33
35.41 1.9 9.2569 36.65 2013 10 21377 377 377 23.71 35.32
0.6 9.3410 7.44 2013 10 21378 378 378 23.76 35.42 0.5
9.4194 27.00 2013 10 21379 379 379 23.33 35.38 1.8 9.4528
38.67 2013 10 21380 380 380 23.26 35.32 1.6 9.4764 58.29 2013
10 22381 381 381 23.55 35.20 1.3 9.4785 18.00 2013 10 22382
382 382 23.39 35.43 1.0 9.5104 3.47 2013 10 22383 383 383
23.32 35.41 1.5 9.5833 37.93 2013 10 22384 384 384 23.35
35.38 1.4 9.6458 40.47 2013 10 22385 385 385 22.89 35.34
1.6 9.8174 34.86 2013 10 22386 386 386 22.92 35.97 1.0
9.9507 44.69 2013 10 22387 387 387 23.62 35.20 0.9 10.1132
15.00 2013 10 22388 388 388 23.42 34.97 1.6 10.1299 2.98 2013
10 22389 389 389 23.31 35.44 1.7 10.1424 45.64 2013 10 22390
390 390 23.66 35.18 2.1 10.1688 9.00 2013 10 22391 391 391
23.67 35.19 0.8 10.3493 9.00 2013 10 22392 392 392 23.02
35.64 1.1 10.4139 58.33 2013 10 22393 393 393 23.35 35.37
1.4 10.4160 34.40 2013 10 22394 394 394 23.25 35.31 1.4
10.5569 36.36 2013 10 22395 395 395 23.32 35.30 1.6 11.2521
26.49 2013 10 23396 396 396 32.24 35.90 1.7 11.5333 12.50 2013
10 23397 397 397 23.37 35.42 1.9 11.5799 39.90 2013 10 23398
398 398 23.36 35.39 1.2 11.5847 8.49 2013 10 24399 399 399
23.31 35.38 1.9 11.6813 42.41 2013 10 24400 400 400 23.25
35.32 1.4 11.6979 39.37 2013 10 24401 401 401 23.35 35.40
2.3 11.7660 5.40 2013 10 24402 402 402 23.32 35.37 1.4
11.9347 19.20 2013 10 24403 403 403 23.59 35.17 1.5 12.0813
16.42 2013 10 24404 404 404 23.44 35.58 1.1 12.0965 11.09 2013
10 24405 405 405 23.31 35.38 1.3 12.1104 8.21 2013 10 24406
406 406 23.34 35.39 2.5 12.1701 2.62 2013 10 24407 407 407
23.19 35.21 2.0 12.2215 46.69 2013 10 24408 408 408 23.75
35.16 1.5 12.2333 18.00 2013 10 24409 409 409 23.28 35.20
1.7 12.4465 77.29 2013 10 24410 410 410 23.33 35.35 1.3
12.8389 10.83 2013 10 24411 411 411 23.31 34.70 2.0 12.9160
5.93 2013 10 24412 412 412 23.27 35.16 1.9 13.0076 32.41 2013
10 25413 413 413 23.36 35.39 1.9 13.1660 34.79 2013 10 25414
414 414 23.55 35.86 1.4 13.3111 26.56 2013 10 25415 415 415
23.27 35.29 2.3 13.3444 35.93 2013 10 25416 416 416 23.61
35.20 1.1 13.5153 15.00 2013 10 25417 417 417 23.44 34.98
1.9 13.8931 16.28 2013 10 25418 418 418 23.39 34.99 1.8
14.1903 28.00 2013 10 26419 419 419 23.35 35.41 1.6 14.1910
59.26 2013 10 26420 420 420 23.33 35.06 1.5 14.3750 7.63 2013
10 26421 421 421 23.61 35.20 1.4 14.4833 18.00 2013 10 26422
422 422 22.62 35.76 1.8 14.7875 36.50 2013 10 26423 423 423
23.35 35.43 1.8 15.1333 25.33 2013 10 27424 424 424 22.92
35.34 1.9 15.2944 0.00 2013 10 27425 425 425 23.26 35.21
1.9 15.3910 6.93 2013 10 27426 426 426 23.34 35.38 1.0
15.4729 5.26 2013 10 27427 427 427 22.53 35.16 2.3 15.5250
8.60 2013 10 27428 428 428 23.37 34.97 1.9 15.5326 20.62 2013
10 27429 429 429 23.34 35.40 2.4 15.6368 36.69 2013 10 28430
430 430 23.33 35.39 1.6 15.7382 30.83 2013 10 28431 431 431
23.34 35.43 1.5 15.7431 26.27 2013 10 28432 432 432 23.88
35.04 1.0 15.7583 20.20 2013 10 28433 433 433 23.39 35.51
1.6 15.8576 26.04 2013 10 28434 434 434 22.90 35.42 1.6
15.8743 29.02 2013 10 28435 435 435 23.49 35.84 1.8 15.9833
16.40 2013 10 28436 436 436 23.36 35.06 1.5 16.0042 28.00 2013
10 28437 437 437 23.30 35.39 1.3 16.1750 21.63 2013 10 28438
438 438 23.97 35.18 1.0 16.1889 6.84 2013 10 28439 439 439
23.35 35.41 1.2 16.2681 16.10 2013 10 28440 440 440 23.47
35.33 0.5 16.3785 42.78 2013 10 28441 441 441 23.42 35.44
1.4 16.5069 5.73 2013 10 28442 442 442 24.00 34.88 1.3
16.5132 9.00 2013 10 28443 443 443 23.27 35.33 1.7 16.5153
30.58 2013 10 28444 444 444 23.29 35.30 1.4 16.5382 35.60 2013
10 29445 445 445 23.70 35.20 0.8 16.5625 15.11 2013 10 29446
446 446 23.63 35.20 0.8 16.5701 18.00 2013 10 29447 447 447
22.94 35.44 1.8 16.5806 50.16 2013 10 29448 448 448 23.71
35.21 0.7 16.6007 15.00 2013 10 29449 449 449 23.63 35.19
0.7 16.6819 15.00 2013 10 29450 450 450 23.72 35.52 1.0
16.8340 12.00 2013 10 29451 451 451 23.74 35.85 0.8 16.8583
16.44 2013 10 29452 452 452 22.87 35.62 2.3 16.9125 73.00 2013
10 29453 453 453 23.69 35.23 1.1 16.9194 27.00 2013 10 29454
454 454 23.67 35.21 0.6 16.9208 18.00 2013 10 29455 455 455
23.62 35.79 1.1 16.9264 32.00 2013 10 29456 456 456 22.91
34.97 1.9 17.0014 82.72 2013 10 29457 457 457 23.38 35.45
1.7 17.0437 33.60 2013 10 29458 458 458 23.36 35.36 1.0
17.1944 36.11 2013 10 29459 459 459 23.06 35.91 1.1 17.3389
47.81 2013 10 29460 460 460 23.39 35.34 2.8 17.3778 35.87 2013
10 29461 461 461 23.20 35.20 1.2 17.4278 7.53 2013 10 29462
462 462 23.30 35.40 1.2 17.5403 40.55 2013 10 29463 463 463
23.62 35.35 0.9 17.5618 35.59 2013 10 29464 464 464 22.89
35.32 3.1 17.6160 13.39 2013 10 29465 465 465 23.32 35.02
1.2 17.6222 33.39 2013 10 30466 466 466 23.34 35.00 2.2
17.6347 36.80 2013 10 30467 467 467 23.94 35.00 1.3 17.6813
27.00 2013 10 30468 468 468 23.78 35.00 1.0 17.8132 39.98 2013
10 30469 469 469 23.01 35.97 1.5 17.9444 9.00 2013 10 30470
470 470 23.28 35.27 0.9 18.1618 30.11 2013 10 30471 471 471
23.28 35.29 1.4 18.2090 28.90 2013 10 30472 472 472 23.30
35.36 2.2 18.2250 50.64 2013 10 30473 473 473 23.29 35.31
1.0 18.2257 34.01 2013 10 30474 474 474 23.97 34.92 1.3
18.4569 9.10 2013 10 30475 475 475 23.34 35.42 3.0 18.8882
38.92 2013 10 30476 476 476 23.34 35.37 2.5 18.9639 39.66 2013
10 30477 477 477 22.83 35.63 1.5 19.0549 43.85 2013 10 31478
478 478 23.72 35.51 1.2 19.3847 12.09 2013 10 31479 479 479
23.42 35.65 1.7 19.4222 21.66 2013 10 31> df <-
data.frame(data)>
coord.type=c("m","m","s","s","m","t","m","m","m","m")>
X <- ppx(data=df,
coord.type=c("m","m","s","s","m","t","m","m","m","m"))>
XMultidimensional point pattern479 points 2-dimensional space coordinates
(longitude,latitude)1-dimensional time coordinates (time)7 columns of marks:
‘X’, ‘No’, ‘magnitude’, ‘depth’, ‘year’, ‘month’ and ‘day’ > plot(X)Loading
required package: scatterplot3d> win <- owin(c(32, 36), c(22, 24))>
plot(X[, 2:3])> plot(win, add=TRUE)> proj <- long2flat(X, win)>
tperiod <- c(0, 20)> param01 <- c(0.0,+ 63.348, + 0.038209,+ 2.6423,+
1.0169)> res <- etas(proj$X, proj$win, tperiod, m0=1.6, param0=param01,
no.itr=1)Error in if (tstart2 < tstart || tstart2 > tend || tlength +
tstart > : missing value where TRUE/FALSE needed>
Thank you in advance,Katerina
[[alternative HTML version deleted]]