Skip to content

Commit

Permalink
Re-run LV yet again
Browse files Browse the repository at this point in the history
  • Loading branch information
nathanaelbosch committed Nov 13, 2023
1 parent afb7845 commit 851984a
Show file tree
Hide file tree
Showing 19 changed files with 6,842 additions and 6,818 deletions.
16 changes: 11 additions & 5 deletions benchmarks/lotkavolterra.jmd
Original file line number Diff line number Diff line change
Expand Up @@ -251,7 +251,7 @@ plot(wp, x=:L2, color=[1 1 1 1 2 2 2 2])
## EK0 vs. EK1: Calibration
Final time-point:
```julia
plot(wp, x=:final, y=:chi2_final, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:final, y=:chi2_final, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (final)")

# Should be distributed according to a Chi-squared distribution:
low, high, mid = quantile(Chisq(2), [0.01, 0.99])..., mean(Chisq(2))
Expand All @@ -262,15 +262,15 @@ hline!([mid], linestyle=:solid, color=:black, label="")

Discrete time-series:
```julia
plot(wp, x=:l2, y=:chi2_steps, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:l2, y=:chi2_steps, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (discrete steps)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
```

Interpolation:
```julia
plot(wp, x=:L2, y=:chi2_interp, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:L2, y=:chi2_interp, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (dense)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
Expand Down Expand Up @@ -316,7 +316,10 @@ plot(wp, color=[2 2 3 3 4 4 5 5])

Calibration:
```julia
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3 4 4 5 5])
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3 4 4 5 5], yguide="Chi-squared (final)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
```

### EK1
Expand Down Expand Up @@ -353,7 +356,10 @@ plot(wp, color=[2 2 3 3])

Calibration:
```julia
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3])
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3], yguide="Chi-squared (final)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
```

## Comparison of the different initialization schemes
Expand Down
734 changes: 369 additions & 365 deletions docs/src/benchmarks/figures/lotkavolterra_10_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
914 changes: 457 additions & 457 deletions docs/src/benchmarks/figures/lotkavolterra_11_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
846 changes: 423 additions & 423 deletions docs/src/benchmarks/figures/lotkavolterra_12_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
762 changes: 381 additions & 381 deletions docs/src/benchmarks/figures/lotkavolterra_13_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
696 changes: 348 additions & 348 deletions docs/src/benchmarks/figures/lotkavolterra_14_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,081 changes: 543 additions & 538 deletions docs/src/benchmarks/figures/lotkavolterra_15_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
610 changes: 304 additions & 306 deletions docs/src/benchmarks/figures/lotkavolterra_16_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,061 changes: 533 additions & 528 deletions docs/src/benchmarks/figures/lotkavolterra_17_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
2,232 changes: 1,116 additions & 1,116 deletions docs/src/benchmarks/figures/lotkavolterra_18_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
266 changes: 133 additions & 133 deletions docs/src/benchmarks/figures/lotkavolterra_2_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
608 changes: 305 additions & 303 deletions docs/src/benchmarks/figures/lotkavolterra_3_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
598 changes: 300 additions & 298 deletions docs/src/benchmarks/figures/lotkavolterra_4_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
602 changes: 302 additions & 300 deletions docs/src/benchmarks/figures/lotkavolterra_5_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
628 changes: 314 additions & 314 deletions docs/src/benchmarks/figures/lotkavolterra_6_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
616 changes: 306 additions & 310 deletions docs/src/benchmarks/figures/lotkavolterra_7_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
616 changes: 306 additions & 310 deletions docs/src/benchmarks/figures/lotkavolterra_8_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
758 changes: 380 additions & 378 deletions docs/src/benchmarks/figures/lotkavolterra_9_1.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
16 changes: 11 additions & 5 deletions docs/src/benchmarks/lotkavolterra.md
Original file line number Diff line number Diff line change
Expand Up @@ -290,7 +290,7 @@ plot(wp, x=:L2, color=[1 1 1 1 2 2 2 2])
## EK0 vs. EK1: Calibration
Final time-point:
```julia
plot(wp, x=:final, y=:chi2_final, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:final, y=:chi2_final, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (final)")

# Should be distributed according to a Chi-squared distribution:
low, high, mid = quantile(Chisq(2), [0.01, 0.99])..., mean(Chisq(2))
Expand All @@ -305,7 +305,7 @@ hline!([mid], linestyle=:solid, color=:black, label="")

Discrete time-series:
```julia
plot(wp, x=:l2, y=:chi2_steps, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:l2, y=:chi2_steps, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (discrete steps)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
Expand All @@ -317,7 +317,7 @@ hline!([mid], linestyle=:solid, color=:black, label="")

Interpolation:
```julia
plot(wp, x=:L2, y=:chi2_interp, color=[1 1 1 1 2 2 2 2])
plot(wp, x=:L2, y=:chi2_interp, color=[1 1 1 1 2 2 2 2], yguide="Chi-squared (dense)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
Expand Down Expand Up @@ -371,7 +371,10 @@ plot(wp, color=[2 2 3 3 4 4 5 5])

Calibration:
```julia
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3 4 4 5 5])
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3 4 4 5 5], yguide="Chi-squared (final)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
```

![](figures/lotkavolterra_15_1.svg)
Expand Down Expand Up @@ -416,7 +419,10 @@ plot(wp, color=[2 2 3 3])

Calibration:
```julia
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3])
plot(wp, x=:final, y=:chi2_final, color=[2 2 3 3], yguide="Chi-squared (final)")
hline!([low, high], linestyle=:dash, color=:black, label="",
fill_between=true, fillcolor=:green, fillalpha=0.15)
hline!([mid], linestyle=:solid, color=:black, label="")
```

![](figures/lotkavolterra_17_1.svg)
Expand Down

0 comments on commit 851984a

Please sign in to comment.