LoadError с использованием приблизительных байесовских критериев

Я получаю сообщение об ошибке, которое меня смущает.

using DifferentialEquations
using RecursiveArrayTools # for VectorOfArray
using DiffEqBayes
f2 = @ode_def_nohes LotkaVolterraTest begin
  dx = x*(1 - x - A*y)
  dy = rho*y*(1 - B*x - y)
end A B rho

u0 = [1.0;1.0]
tspan = (0.0,10.0)
p = [0.2,0.5,0.3]
prob = ODEProblem(f2,u0,tspan,p)
sol = solve(prob,Tsit5())
t = collect(linspace(0,10,200))
randomized = VectorOfArray([(sol(t[i]) + .01randn(2)) for i in 1:length(t)])
data = convert(Array,randomized)

priors = [Uniform(0.0, 2.0), Uniform(0.0, 2.0), Uniform(0.0, 2.0)]
bayesian_result_abc = abc_inference(prob, Tsit5(), t, data,
priors;num_samples=500)

Возвращает ошибку

ERROR: LoadError: DimensionMismatch("first array has length 400 which does not match the length of the second, 398.")
while loading..., in expression starting on line 20.

Мне не удалось найти ни одного массива размером 400 или 398. Спасибо за помощь.


person Gregory    schedule 20.08.2018    source источник
comment
Для справки: решение строит массив 2x200, но из-за ошибки не сохраняет время начала, поэтому это массив 2x199. Это 400 и 398.   -  person Chris Rackauckas    schedule 22.08.2018
comment
Спасибо Крис! Я сходил с ума.   -  person Gregory    schedule 22.08.2018


Ответы (1)


Взгляните на https://github.com/JuliaDiffEq/DiffEqBayes.jl/issues/52 из-за ошибки при передаче t. Это было исправлено на мастере, поэтому вы можете использовать это или подождать некоторое время, скоро у нас будет новый выпуск с обновлениями 1.0, в которых это также будет исправлено.

Спасибо!

person Vaibhav Dixit    schedule 21.08.2018
comment
Я просто опубликую это прямо сейчас и исправлю проблему. Pkg.update() примерно через полчаса. - person Chris Rackauckas; 22.08.2018