Visualization
LikelihoodProfiler.get_interval function outputs estimated confidence interval along with other data as LikelihoodProfiler.ParamInterval structure.
LikelihoodProfiler provides a @recipe for Plots.jl to visualize confidence interval estimation and plot parameter profile based on LikelihoodProfiler.ParamInterval.
using LikelihoodProfiler
# Likelihood function
f(x) = 5.0 + (x[1]-3.0)^2 + (x[1]-x[2]-1.0)^2 + 0*x[3]^2
# Calculate parameters intervals for x[1], x[2], x[3]
res = [
get_interval(
[3., 2., 2.1],
i,
f,
:CICO_ONE_PASS;
loss_crit = 9.
) for i in 1:3]
# Plot parameter profile x[2]
using Plots
plotly()
plot(res[2])
To make a smooth plot compute more profile points using LikelihoodProfiler.update_profile_points! which internally uses PlotUtils.adapted_grid
update_profile_points!(res[2])
plot(res[2])