# estimate model
d4 <- din(sim.dino, q.matr = sim.qmatrix)
# compute tetrachoric correlations
skill.cor(d4)
## estimated tetrachoric correlations
## $cor.skills
## V1 V2 V3
## V1 1.0000000 0.2567718 0.2552958
## V2 0.2567718 1.0000000 0.9842188
## V3 0.2552958 0.9842188 1.0000000
Run the code above in your browser using DataLab