LRtest-分割向量不允许NA

Question

我目前正在运行LRtest（Rasch-Model，eRm-R-package），并且需要使用外部分割标准（每平方公里人口密度的中位数^ 2）分割我的小组。 一些人没有回答，R错误是：

lrt8 <- LRtest (res, splitcr = splitresid)

  

LRtest.Rm（res，splitcr = splitresid）中的错误：分割向量应该   不包含NA

我该如何排除这些NA或将其替换为中位数等？

Ps：我的测试运行时间超长（对于400个人中有100个项目的3小时测试），这是正常现象吗？

Answer 1

一个选项是

library(eRm)
library(zoo)
res <- raschdat1_RM_fitted   # using the example from `?LRtest`
splitvec <- sample(1:2, 100, replace = TRUE)
splitvec[2:3] <- NA # insert some NA
lrres <- LRtest(res, splitcr = splitvec)

  

LRtest.Rm（res，splitcr = splitvec）中的错误：         分割向量中不得包含NA

使用NA中的median将na.aggregate替换为zoo

lrres <- LRtest(res, splitcr = as.integer(na.aggregate(splitvec, median)))