with the same data set, I get different omega reliability estimates (when looking for omega hierarchical from both) with the psych and the MBESS packages.
ci.reliability(data = (full_data[, my_vector]),
type = "hierarchical", interval.type = "bca", B = 10000, conf.level = 0.95)
both return the output exactly as I would expect in terms of content, no errors etc, but the values are very different
Has anyone else encountered this? Any ideas why?