Skip to content

Commit 0c9535c

Browse files
committed
reference DiffEqFlux and mentione final parameter set.
1 parent 0154d6b commit 0c9535c

File tree

1 file changed

+4
-1
lines changed

1 file changed

+4
-1
lines changed

docs/src/tutorials/parameter_estimation.md

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,5 @@
11
# Parameter Estimation
2-
The parameters of a model, generated by Catalyst, can be fitted using various packages available in the Julia ecosystem. Refer [here](https://diffeq.sciml.ai/stable/analysis/parameter_estimation/) for a more detailed description. Below follows a quick tutorial of how a parameter set for a model can be fitted to data.
2+
The parameters of a model, generated by Catalyst, can be estimated using various packages available in the Julia ecosystem. Refer [here](https://diffeq.sciml.ai/stable/analysis/parameter_estimation/) for a more detailed description. Below follows a quick tutorial of how [DiffEqFlux](https://diffeqflux.sciml.ai/dev/) can be used to fit a parameter set to data.
33

44
First, we fetch the required packages.
55
```julia
@@ -88,6 +88,9 @@ plot!(sol_estimate,color=[:darkblue :darkred], linestyle=:dash,lw=3,label=["X es
8888
```
8989
![parameter_estimation_plot4](../assets/parameter_estimation_plot4.svg)
9090

91+
The final parameter set becomes `[0.9996559014056948, 2.005632696191224]` (the real one was [1.0, 2.0]`).
92+
93+
9194
The reason we chose to fit the model on a smaller interval to begin with, and then extend the interval, is to avoid getting stuck in a local minimum. Here specifically, we chose our initial interval to be smaller than a full cycle of the oscillation. If we had chosen to fit a parameter set on the full interval immediately we would have received an inferior solution.
9295
```julia
9396
p_estimate = optimise_p([5.,5.],30.).minimizer

0 commit comments

Comments
 (0)