# transform probability of success to delta
# based on saturated GDINA model
# assuming an item with 2 attributes has the following
# probability of success for 00, 10, 01 and 11
pj <- c(0.2,0.4,0.6,0.8)
Mj <- designmatrix(2)
# delta parameters can be calculated in this way:
deltaj <- solve(Mj)%*%pj
# for reduced CDMs, OLS or WLS may be used
Run the code above in your browser using DataLab