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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
