# Unable to fit example Turing model

**URL:** https://discourse.julialang.org/t/unable-to-fit-example-turing-model/37764
**Category:** Probabilistic Programming
**Created:** [April 17, 2020, 2:44pm UTC](https://discourse.julialang.org/t/unable-to-fit-example-turing-model/37764 "2020-04-17T14:44:43Z")
**Posts on this page:** 2
**Page:** 1

<div class="post-metadata">

### Author: ![jacobcvt12](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/jacobcvt12/32/2695_2.png) [@jacobcvt12](https://discourse.julialang.org/u/jacobcvt12)
#### Post date: [April 17, 2020, 2:44pm UTC](https://discourse.julialang.org/t/unable-to-fit-example-turing-model/37764/1 "2020-04-17T14:44:43Z")

</div>

The following code, copied from the Turing website throws an error

```julia
using Turing
using StatsPlots

# Define a simple Normal model with unknown mean and variance.
@model gdemo(x, y) = begin
  s ~ InverseGamma(2, 3)
  m ~ Normal(0, sqrt(s))
  x ~ Normal(m, sqrt(s))
  y ~ Normal(m, sqrt(s))
end

# Run sampler, collect results
chn = sample(gdemo(1.5, 2), HMC(0.1, 5), 1000)

```

```julia
ERROR: MethodError: no method matching AxisArrays.AxisArray(::Array{Float64,3}; iter
=1:1:1000, var=["m", "s", "lp", "n_steps", "is_accept", "acceptance_rate", "log_dens
ity", "hamiltonian_energy", "hamiltonian_energy_error", "step_size", "nom_step_size"
], chain=1:1)

```

I am on Julia 1.4.1 and Turing 0.10.1.

---

<div class="post-metadata">

### Author: ![jacobcvt12](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/jacobcvt12/32/2695_2.png) [@jacobcvt12](https://discourse.julialang.org/u/jacobcvt12)
#### Post date: [April 17, 2020, 2:49pm UTC](https://discourse.julialang.org/t/unable-to-fit-example-turing-model/37764/2 "2020-04-17T14:49:43Z")

</div>

I duplicated an issue on the Turing GitHub [here](https://github.com/TuringLang/Turing.jl/issues/1220). The solution was to update AxisArrays to 0.4.2.
