# Turing.jl setting different updates for arrays of variables

**URL:** <https://discourse.julialang.org/t/turing-jl-setting-different-updates-for-arrays-of-variables/92330>\
**Category:** Probabilistic Programming\
**Created:** [December 30, 2022, 7:14pm UTC](https://discourse.julialang.org/t/turing-jl-setting-different-updates-for-arrays-of-variables/92330 "2022-12-30T19:14:30Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![dariogasbarra](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/dariogasbarra/32/223397_2.png) [@dariogasbarra](https://discourse.julialang.org/u/dariogasbarra)\
**Post date:** [December 30, 2022, 7:14pm UTC](https://discourse.julialang.org/t/turing-jl-setting-different-updates-for-arrays-of-variables/92330/1 "2022-12-30T19:14:30Z")

</div>

Hello I have a model which has an array of discrete binary parameters J[t]  
and an array of continuous parameters X[t]. I want to use NUTS() or  
HMC() updates for the continuous variables and MH() updates  
for the discrete parameters in the McMC. It should be possible  
by using GIbbs(), however it is not clear how  
to use it when the random variables are indexed arrays, I could not find  
not such examples in the documentation.

Here is the changepoint model

Turing.@model function poisson\_model(N,T,p,rho)  
J = tzeros(Int,T)  
X = tzeros(Real,T)  
J[1] ~ Dirac(1)  
X[1] ~ Exponential(rho)  
theta0= X[1]  
N[1] ~ Poisson(theta0)  
for t in 2:T  
J[t] ~ Bernoulli(p)  
X[t] ~ Laplace(0,rho)  
theta = sum( (J[s]\*X[s]) for s in 1:t)  
if theta \< 0.0  
Turing.@addlogprob! -Inf  
return  
end  
N[t] ~ Poisson(theta)  
end  
return(J,X)  
end

poisson\_posterior= poisson\_model(N,T,p,rho)

chain = Turing.sample(poisson\_posterior,MH( ),nMcMC)

works but it is mixing very slowly. Then if I try to combine  
different updates using GIbbs()

chain = Turing.sample(poisson\_posterior,Gibbs(HMC(0.1,5,:X),MH(:J)),nMcMC)

the last command gives errors, while it should work when X and J are scalars.

Thanks !

---

<div class="post-metadata">

**Author:** ![Dan](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/dan/32/42581_2.png) [@Dan](https://discourse.julialang.org/u/Dan)\
**Post date:** [December 30, 2022, 8:06pm UTC](https://discourse.julialang.org/t/turing-jl-setting-different-updates-for-arrays-of-variables/92330/2 "2022-12-30T20:06:13Z")

</div>

It is better to surround code with triple quotes (```) to make it format better.  
Also, it is better to have code which runs when copy pasted into REPL (known as a MWE - minimal working example). In this case:

```julia
using Turing
Turing.@model function poisson_model(N, T, p, rho)
    J = tzeros(Int, T)
    X = tzeros(Real, T)
    J[1] ~ Dirac(1)
    X[1] ~ Exponential(rho)
    theta0 = X[1]
    N[1] ~ Poisson(theta0)
    for t = 2:T
        J[t] ~ Bernoulli(p)
        X[t] ~ Laplace(0, rho)
        theta = sum((J[s] * X[s]) for s = 1:t)
        if theta < 0.0
            Turing.@addlogprob! -Inf
            return
        end
        N[t] ~ Poisson(theta)
    end
    return (J, X)
end
rho = 0.5
T = 3
N = [1, 2, 2]
p = 0.5
poisson_posterior = poisson_model(N, T, p, rho)
nMCMC = 1000
using StatsPlots
chain = Turing.sample(poisson_posterior, Gibbs(HMC(0.1, 5, :X), MH(:J)), nMCMC)
plot(chain)

```

As for the question, the above code ran without errors on Julia 1.8.3 with Turing 0.22.0. What is the error on your setup?

---

<div class="post-metadata">

**Author:** ![dariogasbarra](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/dariogasbarra/32/223397_2.png) [@dariogasbarra](https://discourse.julialang.org/u/dariogasbarra)\
**Post date:** [December 30, 2022, 8:44pm UTC](https://discourse.julialang.org/t/turing-jl-setting-different-updates-for-arrays-of-variables/92330/3 "2022-12-30T20:44:46Z")

</div>

Thanks a lot Dan, next time I will use ‘’’ and put a minimal working example. I was not sure whether the syntax  
Gibbs( NUTS(100,0.65,:X),MH(:J)) would be correct when X[t] and J[t] are arrays of random variables, but yes it is correct , the code is running now. The problem was perhaps in too large HMC stepsize for my data  
N = [182,165,182,136,135,127,136,133,  
132,112,114,116,99,80,129,111,118,107,94,90,101,103,  
97,120,93,126,148,125,120,171,155,135,142,126,158,  
130,137,149,134,134,159,132,138,158,160,154,173,166,  
165,150,170,142,125,145,139,154,133,99,130,105,107,  
116,119,110,102,99,73,78,78,70,63,61]