diff --git a/README.md b/README.md index 3f9fffd..b013dd9 100644 --- a/README.md +++ b/README.md @@ -276,14 +276,14 @@ containers_bwd = CSV.read(joinpath("data", "containers-bwd.csv"), DataFrame); first(containers_bwd, 6) #> 6×6 DataFrame #> Row │ timestamp obs sensor_id receiver_x receiver_y radius -#> │ String31 Int64 Int64 Float64 Float64 Int64 -#> ─────┼─────────────────────────────────────────────────────────────────────── -#> 1 │ 2016-03-17 01:52:00 1 26 7.09242e5 6.25311e6 1500 -#> 2 │ 2016-03-17 01:54:00 1 26 7.09242e5 6.25311e6 1500 -#> 3 │ 2016-03-17 01:56:00 1 26 7.09242e5 6.25311e6 1500 -#> 4 │ 2016-03-17 01:58:00 1 26 7.09242e5 6.25311e6 2250 -#> 5 │ 2016-03-17 02:00:00 1 26 7.09242e5 6.25311e6 1500 -#> 6 │ 2016-03-17 02:02:00 1 26 7.09242e5 6.25311e6 1500 +#> │ String31 Int64 Int64 Float64 Float64 Float64 +#> ─────┼──────────────────────────────────────────────────────────────────────── +#> 1 │ 2016-03-17 01:52:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 2 │ 2016-03-17 01:54:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 3 │ 2016-03-17 01:56:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 4 │ 2016-03-17 01:58:00 1 26 7.09242e5 6.25311e6 2250.0 +#> 5 │ 2016-03-17 02:00:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 6 │ 2016-03-17 02:02:00 1 26 7.09242e5 6.25311e6 1500.0 # Read archival (depth) observations archival = CSV.read(joinpath("data", "archival.csv"), DataFrame); @@ -339,27 +339,27 @@ datasets_bwd = [acoustics, containers_bwd, archival] #> 12959 │ 2016-03-18T01:48:00 31 0 7.09742e5 6.25331e6 #> 12960 │ 2016-03-18T01:48:00 33 0 7.07342e5 6.25341e6 #> 3 columns and 12945 rows omitted -#> 360×6 DataFrame -#> Row │ timestamp obs sensor_id receiver_x receiver_y radius -#> │ DateTime Int64 Int64 Float64 Float64 Int64 -#> ─────┼─────────────────────────────────────────────────────────────────────── -#> 1 │ 2016-03-17T01:52:00 1 26 7.09242e5 6.25311e6 1500 -#> 2 │ 2016-03-17T01:54:00 1 26 7.09242e5 6.25311e6 1500 -#> 3 │ 2016-03-17T01:56:00 1 26 7.09242e5 6.25311e6 1500 -#> 4 │ 2016-03-17T01:58:00 1 26 7.09242e5 6.25311e6 2250 -#> 5 │ 2016-03-17T02:00:00 1 26 7.09242e5 6.25311e6 1500 -#> 6 │ 2016-03-17T02:02:00 1 26 7.09242e5 6.25311e6 1500 -#> 7 │ 2016-03-17T02:04:00 1 26 7.09242e5 6.25311e6 2250 -#> 8 │ 2016-03-17T02:06:00 1 26 7.09242e5 6.25311e6 1500 -#> ⋮ │ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ -#> 354 │ 2016-03-17T13:14:00 1 26 7.09242e5 6.25311e6 95250 -#> 355 │ 2016-03-17T13:16:00 1 26 7.09242e5 6.25311e6 96000 -#> 356 │ 2016-03-17T13:18:00 1 26 7.09242e5 6.25311e6 96750 -#> 357 │ 2016-03-17T13:20:00 1 26 7.09242e5 6.25311e6 97500 -#> 358 │ 2016-03-17T13:22:00 1 26 7.09242e5 6.25311e6 98250 -#> 359 │ 2016-03-17T13:24:00 1 26 7.09242e5 6.25311e6 99000 -#> 360 │ 2016-03-17T13:26:00 1 26 7.09242e5 6.25311e6 99750 -#> 345 rows omitted +#> 731×6 DataFrame +#> Row │ timestamp obs sensor_id receiver_x receiver_y radius ⋯ +#> │ DateTime Int64 Int64 Float64 Float64 Float64 ⋯ +#> ─────┼────────────────────────────────────────────────────────────────────────── +#> 1 │ 2016-03-17T01:52:00 1 26 7.09242e5 6.25311e6 1500.0 ⋯ +#> 2 │ 2016-03-17T01:54:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 3 │ 2016-03-17T01:56:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 4 │ 2016-03-17T01:58:00 1 26 7.09242e5 6.25311e6 2250.0 +#> 5 │ 2016-03-17T02:00:00 1 26 7.09242e5 6.25311e6 1500.0 ⋯ +#> 6 │ 2016-03-17T02:02:00 1 26 7.09242e5 6.25311e6 1500.0 +#> 7 │ 2016-03-17T02:04:00 1 26 7.09242e5 6.25311e6 2250.0 +#> 8 │ 2016-03-17T02:06:00 1 26 7.09242e5 6.25311e6 1500.0 +#> ⋮ │ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱ +#> 725 │ 2016-03-18T01:36:00 1 26 7.09242e5 6.25311e6 373500.0 ⋯ +#> 726 │ 2016-03-18T01:38:00 1 26 7.09242e5 6.25311e6 374250.0 +#> 727 │ 2016-03-18T01:40:00 1 26 7.09242e5 6.25311e6 375000.0 +#> 728 │ 2016-03-18T01:42:00 1 26 7.09242e5 6.25311e6 375750.0 +#> 729 │ 2016-03-18T01:44:00 1 26 7.09242e5 6.25311e6 376500.0 ⋯ +#> 730 │ 2016-03-18T01:46:00 1 26 7.09242e5 6.25311e6 377250.0 +#> 731 │ 2016-03-18T01:48:00 1 26 7.09242e5 6.25311e6 378000.0 +#> 716 rows omitted #> 720×5 DataFrame #> Row │ timestamp sensor_id obs depth_sigma depth_deep_eps #> │ DateTime Int64 Float64 Int64 Int64 @@ -427,7 +427,7 @@ fwd = particle_filter(timeline = timeline, xinit = xinit, yobs = yobs_fwd, model_move = model_move, - n_record = 500, + n_record = 1000, direction = "forward"); # Simulate initial states for the backward filter @@ -447,7 +447,7 @@ bwd = particle_filter(timeline = timeline, xinit = xinit, yobs = yobs_bwd, model_move = model_move, - n_record = 500, + n_record = 1000, direction = "backward"); ``` @@ -468,7 +468,7 @@ and after each step): ``` julia # (optional) Set vmap to improve speed here -n_particle = 100; +n_particle = 750; smo = two_filter_smoother(timeline = timeline, xfwd = fwd.state[1:n_particle, :], xbwd = bwd.state[1:n_particle, :], diff --git a/docs/figures/README-unnamed-chunk-13-2.png b/docs/figures/README-unnamed-chunk-13-2.png index d789899..4a0db75 100644 Binary files a/docs/figures/README-unnamed-chunk-13-2.png and b/docs/figures/README-unnamed-chunk-13-2.png differ