if (FALSE) {
# PD matrix
mod <- mirt(Science, 1, SE=TRUE)
info <- solve(vcov(mod)) ## observed information
secondOrderTest(info)
secondOrderTest(info, method = 'chol')
secondOrderTest(info, method = 'det')
# non-PD matrix
mat <- matrix(c(1,0,0,0,1,1,0,1,1), ncol=3)
mat
secondOrderTest(mat)
secondOrderTest(mat, method = 'chol')
secondOrderTest(mat, method = 'det')
}
Run the code above in your browser using DataLab