
; From Harley Montgomery, a MCMC transition probability for a Dirichlet Mixture Model

(herbie-test ((c_p (positive default)) (c_n (positive default)) t s)
  "Harley's example"
  (/ (* (expt (/ (+ 1 (exp (- s)))) e+) (expt (- 1 (/ (+ 1 (exp (- s))))) e-))
     (* (expt (/ (+ 1 (exp (- t)))) e+) (expt (- 1 (/ (+ 1 (exp (- t))))) e-)))
  (* (expt (/ (+ 1 (exp (- t))) (+ 1 (exp (- s)))) e+)
     (expt (/ (+ 1 (exp t)) (+ 1 (exp s))) e-)))


; From Harley Montgomery, a MCMC transition probability for a Dirichlet Mixture Model

(lambda ([c_p (positive default)] [c_n (positive default)] t s)
  #:name "Harley's example"
  (/ (* (expt (/ (+ 1 (exp (- s)))) c_p) (expt (- 1 (/ (+ 1 (exp (- s))))) c_n))
     (* (expt (/ (+ 1 (exp (- t)))) c_p) (expt (- 1 (/ (+ 1 (exp (- t))))) c_n)))
  #:target
  (* (expt (/ (+ 1 (exp (- t))) (+ 1 (exp (- s)))) c_p)
     (expt (/ (+ 1 (exp t)) (+ 1 (exp s))) c_n)))

Fix variable names in Harley's example
Fix variable names in Harley's example
