Fix 173

src/samplers/sampler.jl

@@ -83,7 +83,7 @@ assume{T<:Distribution}(spl::Void, dists::Vector{T}, vn::VarName, var::Any, vi::
     end
   end
 
-  acclogp!(vi, sum(logpdf(dist, rs, istrans(vi, vns[1]))))
+  acclogp!(vi, sum(logpdf_with_trans(dist, rs, istrans(vi, vns[1]))))
 
   var
 end
@@ -95,5 +95,9 @@ observe{T<:Distribution}(spl::Void, dists::Vector{T}, value::Any, vi::VarInfo) =
   @assert length(dists) == 1 "[observe] Turing only support vectorizing i.i.d distribution"
   dist = dists[1]
   @assert isa(dist, UnivariateDistribution) || isa(dist, MultivariateDistribution) "[observe] vectorizing matrix distribution is not supported"
-  acclogp!(vi, sum(logpdf(dist, value)))
+  if isa(dist, UnivariateDistribution)  # only univariate distributions support broadcast operation (logpdf.) by Distributions.jl
+    acclogp!(vi, sum(logpdf.(dist, value)))
+  else
+    acclogp!(vi, sum(logpdf(dist, value)))
+  end
 end

src/samplers/support/transform.jl

@@ -87,7 +87,7 @@ link{T<:Real}(d::RealDistribution, x::Union{T,Vector{T}}) = x
 
 invlink{T<:Real}(d::RealDistribution, x::Union{T,Vector{T}}) = x
 
-logpdf_with_trans{T<:Real}(d::RealDistribution, x::Union{T,Vector{T}}, transform::Bool) = logpdf(d, x)
+logpdf_with_trans{T<:Real}(d::RealDistribution, x::Union{T,Vector{T}}, transform::Bool) = logpdf.(d, x)
 
 
 #########
@@ -103,7 +103,7 @@ link{T<:Real}(d::PositiveDistribution, x::Union{T,Vector{T}}) = log(x)
 invlink{T<:Real}(d::PositiveDistribution, x::Union{T,Vector{T}}) = exp.(x)
 
 logpdf_with_trans{T<:Real}(d::PositiveDistribution, x::Union{T,Vector{T}}, transform::Bool) = begin
-  lp = logpdf(d, x)
+  lp = logpdf.(d, x)
   transform ? lp + log.(x) : lp
 end