Tilde

julia> m = @model (x, s) begin
        σ ~ Exponential()
        @inbounds x[1] ~ Normal(σ = σ)
        n = length(x)
        @inbounds for j = 2:n
            x[j] ~ StudentT(1.5, x[j - 1], σ)
        end
    end;
    
julia> x = zeros(3);

julia> r = rand(m(x,10))
(σ = 9.096155145583953, x = [14.876934886768867, 6.612967845518229, 2.045770246490428])

julia> x
3-element Vector{Float64}:
 14.876934886768867
  6.612967845518229
  2.045770246490428

julia> ℓ = logdensityof(m(x, 1.0) | (;x), (σ = 1.0,))
-122.91114458882001