BMD Single Dichotomous values are not matching

[Sooyeong]

Summary of single model fit results of BMD are not matching with modelfit$bmd result.

BMDL,BMDU seems correct, but BMD estimates are not matching.

See the log below

mData <- matrix(c(0, 2,50,
1, 2,50,
3, 10, 50,
16, 18,50,
32, 18,50,
33, 17,50),nrow=6,ncol=3,byrow=TRUE)
D <- mData[,1]
Y <- mData[,2]
N <- mData[,3]

res<-single_dichotomous_fit(mData[,1],mData[,2],mData[,3],model_type ="hill",fit_type = "mcmc", BMR = 0.1)

res$bmd
BMD BMDL BMDU
3.372234 1.319615 7.559162

summary(res)
Summary of single model fit (MCMC) using ToxicR
Model: Hill

BMD: 2.72 (1.32, 7.56) 90.0% CI

Convergence Diagnostics on BMD

Effective Sample Size: 353.61

Geweke Z-score that mean of first 30% of
MCMC chain is different from last 40%
Z-Score: 3.859 P-value 0.000