I used to work in r and now i’m switching to Julia. In r, I was using DifferentialEvolution package where its very straight forward if you want to give multiplicative or additive noise with respect to data as shown:
likelihood1 = sum(dnorm(x= prediction, mean= data, sd=(data+0.01)*SD,log=TRUE))
I tried same logic in DiffenrentialEvolutionMCMC in julia as it has a direct loglikelihood function but its not working.
for j in 1:3
for i in 1:6
LL += logpdf(MvNormal(predicted[:,i], data[j][:,i] .*σ ), data[j][:,i])
end
end
Also I’m trying in turing.jl
for j in 1:3
for i in 1:6
data[j][:,i] ~ MvNormal(predicted[:,i], data[j][:,i] .*σ)
end
end
I’m pretty sure there is something basic i’m missing here… It will be really nice if someone can help… Thank You