Missing rand procedure for MvNormalMeanScalePrecision

ReactiveBayes / ReactiveMP.jl

High-performance reactive message-passing based Bayesian inference engine

MIT License

101 stars 14 forks source link

Missing rand procedure for MvNormalMeanScalePrecision #387

Open albertpod opened 6 months ago

albertpod commented 6 months ago

Following the discussion, we have found that sampling from MvNormalMeanScalePrecision is not implemented, i.e.

julia> dist = MvNormalMeanScalePrecision(zeros(2), 2.0)
MvNormalMeanScalePrecision{Float64, Vector{Float64}}(μ=[0.0, 0.0], γ=2.0)

julia> rand(dist)
ERROR: MethodError: no method matching length(::MvNormalMeanScalePrecision{Float64, Vector{Float64}})

bvdmitri commented 6 months ago

Maybe we should move this distribution in ExponentialFamily.jl?