gistfile1.jl

julia> quote
       for i in 1:n
         x[i] ~ Bernoulli(p)
       end
       
       p ~ Beta(alpha, beta)
       alpha = 1.0
       beta = 1.0
       end
quote 	#  line 2:
    for i = 1:n	#  line 3:
        ~(x[i], Bernoulli(p))
    end	#  line 6:
    ~(p, Beta(alpha, beta))	#  line 7:
    alpha = 1.0	#  line 8:
    beta = 1.0
end

julia> quote
       for i in 1:n
         x[i] ~ Normal(mu, sigma)
       end
       
       mu ~ Normal(0.0, 10000.0)
       sigma = 1.0
       end
quote 	#  line 2:
    for i = 1:n	#  line 3:
        ~(x[i], Normal(mu, sigma))
    end	#  line 6:
    ~(mu, Normal(0.0, 10000.0))	#  line 7:
    sigma = 1.0
end

julia> quote
       for i in 1:n
         mu[i] = beta[1] * x[i] + beta[2]
         y[i] ~ Normal(mu[i], sigma)
       end
       
       for i in 1:n
         beta[i] ~ Normal(0.0, 10000.0)
       end
       
       sigma ~ Gamma(1.0, 10000.0)
       end
quote 	#  line 2:
    for i = 1:n	#  line 3:
        mu[i] = +(*(beta[1], x[i]), beta[2])	#  line 4:
        ~(y[i], Normal(mu[i], sigma))
    end	#  line 7:
    for i = 1:n	#  line 8:
        ~(beta[i], Normal(0.0, 10000.0))
    end	#  line 11:
    ~(sigma, Gamma(1.0, 10000.0))
end