Pedro
I think this is a bug related to having defined parameters and the new parameter labels. As the HPD is calculated from the MCMC draws, but the defined parameters are calculated outside of it.
For now, you can get the hpd for the all the parameters, except the defined ones from blavaan like this
blavInspect(fit, what = "hpd", add.labels = FALSE)
Or you can calculate them from the posterior draws and get the HPD like this
parTable(fit)[,c("lhs","op","rhs","label","pxnames")]
mc_out <- data.frame(as.matrix(blavInspect(fit, "mcmc", add.labels = FALSE)))
dim(mc_out)
head(mc_out)
mc_out$ab <- mc_out$bet_sign.3. * mc_out$bet_sign.1.
mc_out$total <- (mc_out$bet_sign.3. * mc_out$bet_sign.1.) + mc_out$bet_sign.2.
head(mc_out)
coda:::HPDinterval(as.mcmc(mc_out) )
Hope this helps