cronbach_alpha computes Cronbach's alpha internal consistency reliability
cronbach_alpha
cronbach_alpha(responses)
responses
the oberved responses, 2d matrix
cronbach_alpha(model_3pl_gendata(1000, 20)$u)