validation check

philchalmers · Jan 28, 2025 · 1b2ac32 · 1b2ac32
1 parent bf1e50c
commit 1b2ac32
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 3 deletions.
diff --git a/R/SingleGroup-methods.R b/R/SingleGroup-methods.R
@@ -1884,8 +1884,7 @@ mirt2traditional <- function(x, vcov, nfact){
     x@par <- par
     names(x@est) <- names(par)
     x@parnames <- names(x@par)
-    if(length(vcov) == 0L || (is.na(vcov[1L,1L]) || !(cls %in%
-                                                      c('dich', 'graded', 'gpcm', 'nominal', 'nestlogit', 'Luo2001')))){
+    if(!valid_vcov(vcov) || !(cls %in% c('dich', 'graded', 'gpcm', 'nominal', 'nestlogit', 'Luo2001'))){
         x@SEpar <- numeric()
     } else {
         vpick <- subset_vcov(xorg, vcov)

diff --git a/R/utils.R b/R/utils.R
@@ -420,7 +420,7 @@ SE.Lambdas <- function(pars, acov, nfact){
         out[1, index]
     }
 
-    if(length(acov) == 1 && is.na(acov)) return(NULL)
+    if(!valid_vcov(acov)) return(NULL)
     nitems <- length(pars) - 1
     SE.F <- matrix(NA, nitems, nfact)
     for(i in 1:nitems){
@@ -2919,6 +2919,12 @@ is.latent_regression <- function(mod){
     !is.null(mod@Data$covdata)
 }
 
+valid_vcov <- function(obj){
+    ret <- if(length(obj) == 0 || (length(obj) == 1 && is.na(obj[1]))) FALSE
+      else TRUE
+    ret
+}
+
 # subset of vcov for delta method
 subset_vcov <- function(obj, vcov){
     nms <- colnames(vcov)