Calculate score function for each dichotomous item or each nonzero category for polytomous items Only applicable to saturated model ofr joint attribute distribution
score(object, parm = "delta")an object of class GDINA
Either delta or prob indicating score function for delta parameters and
success probabily parameters
a list where elements give the score functions for each item or category
if (FALSE) {
dat <- sim10GDINA$simdat
Q <- sim10GDINA$simQ
fit <- GDINA(dat = dat, Q = Q, model = "GDINA")
score(fit)
}