Replace bernoullli with bernoulli_logit

[JaeyeongYang]

Since bernoulli function in Stan is based on log1m, it seems to be better to use bernoulli_logit(...) instead of bernoulli(inv_logit(...)). I presume it can reduce initialization issues on hBayesDM.

• Ref: https://discourse.mc-stan.org/t/initialization-between-2-2-failed-after-100-attempts/3941