1

我multipleGroup()对 R 中的函数有一些问题，我真的希望有人可以帮助我。
我想首先创建一个独立的、不受约束的模型，其中两个组的所有参数都可能不同，但我收到以下错误消息：

错误：如果没有适当的约束，将无法识别多组模型（组包含缺失的数据模式，其中项目响应已完全省略，或者每个组中的类别数不等于类别总数）

下面你会看到一个运行良好的示例代码（也适用于多头数据），接下来我会展示我的代码和一些关于我的数据的信息。我看不到我的代码中缺少什么......有没有人可以给我一点小费？

示例代码如下：

lsat <- expand.table(LSAT7)
set.seed(1234)
group <- sample(c('group.1', 'group.2'), 1000, replace=TRUE)
xtabs(~group)

mod <- multipleGroup(lsat, model=1, group=group)

我的代码产生了上述错误消息，是这样的：

Tdata <- import(here::here("test_data.xlsx"))
group <- c(rep('group1',25),rep('group2',25))
xtabs(~group)
Tdata <- Tdata %>% dplyr::select(item1, item2, item3, item4, item5)

multipleGroup(Tdata, model = 1, group = group)

向您显示我的数据的详细信息：

> head(Tdata)
  item1 item2 item3 item4 item5
1     1     2     1     2     1
2     2     2     1     2     3
3     3     4     3     3     3
4     3     4     5     4     5
5     2     1     2     1     2
6     2     2     1     2     2

> summary(Tdata)
     item1          item2          item3          item4          item5    
 Min.   :1.00   Min.   :1.00   Min.   :1.00   Min.   :1.00   Min.   :1.0  
 1st Qu.:2.00   1st Qu.:2.00   1st Qu.:2.00   1st Qu.:2.00   1st Qu.:2.0  
 Median :2.00   Median :3.00   Median :3.00   Median :3.00   Median :3.0  
 Mean   :2.70   Mean   :2.88   Mean   :2.82   Mean   :2.92   Mean   :2.9  
 3rd Qu.:3.75   3rd Qu.:4.00   3rd Qu.:4.00   3rd Qu.:4.00   3rd Qu.:4.0  
 Max.   :5.00   Max.   :5.00   Max.   :5.00   Max.   :5.00   Max.   :5.0  

> class(Tdata)
[1] "data.frame"

> class(Tdata$item1)
[1] "numeric"

> any(is.na(Tdata))
[1] FALSE

如果您需要任何进一步的信息，请告诉我。非常感谢你！

4

0 回答 0