Open willwerscheid opened 5 years ago
I think fixing this is as simple as changing alpha = outer(-x, m$a,FUN="+")/s to alpha = outer(x, -m$b,FUN="+")/s in comp_postmean2 (since we send -b and -a to e2truncgamma) but I need to check that this change wouldn't break other likelihoods.
alpha = outer(-x, m$a,FUN="+")/s
alpha = outer(x, -m$b,FUN="+")/s
comp_postmean2
-b
-a
e2truncgamma
I think fixing this is as simple as changing
alpha = outer(-x, m$a,FUN="+")/s
toalpha = outer(x, -m$b,FUN="+")/s
incomp_postmean2
(since we send-b
and-a
toe2truncgamma
) but I need to check that this change wouldn't break other likelihoods.