# an IPV that comprises the honesty/humility and the agreeableness factor of
# the HEXACO (reduced to first 4 items per facet and first 1000 observations
# to reduce runtime)
res <- ipv_est(
HEXACO[1:500, grep("^H_.*[1-4]$|^A_.*[1-4]$", names(HEXACO))],
"HA")
nested_chart(res$est)
Run the code above in your browser using DataLab