Hi Phil,
i use fscores(..., method = "EAPsum") to get theta estimates for each sum score. however, at theta = 4 there is little increase in theta when the sum score increases and i suspect it is impossible for the estimate to go beyond 4. i tried to adjust with the theta_lim argument but the estimation seems uninfluenced by that (even if i narrow the band to theta = c(-1,1) the expected thetas are the same). I guess, EAPsum just doesn't take theta_lim into account, could you look into that?
Thanks a lot, best wishes, Felix