Как использовать `grm` в пакете` ltm`?

Question

Я пытаюсь запустить grm в ltm пакете.Мой сценарий выглядит следующим образом:

library (ltm)
library (msm)
library (polycor)

dim(data)
head(data)
str(data)
descript(data)

options(max.print=1000000) 
rcor.test(data, method = "pearson")

data_2 <- data
data_2[] <- lapply(data_2, factor)
out <- grm(data_2)
out2 <- grm(data_2, constrained = TRUE)
anova(out2,out)

 margins(out)

Однако, когда я запускаю margins(out), я получаю эту ошибку: Error in exp[ind] <- n * colSums(GHw * pp) : subscript out of bounds

Кто-нибудь, пожалуйста, объясните это?И как я могу решить эту проблему?

У меня есть 35 пунктов в моей анкете и 576 респондентов.Вот пример данных (первые 6 респондентов и первые 6 пунктов).

  pespd_qa1 pespd_qa2 pespd_qa3 pespd_qa4 pespd_qa5 pespd_qa6 
1         9         5         7         4         1         3                 
2         5         0         9         6         0         8                  
3         5         3         5         6         3         5                  
4         7         5         4         3         1         1                  
5         2         3         0         0         0         0                  
6        10         1         8         2         2         5