我的结构如下:
'data.frame': 1041 obs. of 4 variables:
$ RST_AST: num 0.45 0.576 0.45 0.45 0.675 ...
$ subj : chr "4" "13" "24" "25" ...
$ domain : chr "perceptual" "perceptual" "perceptual" "perceptual" ...
$ ce : chr "compromise" "compromise" "compromise" "compromise" ...现在,我想使用这样的lme4包运行一个多级模型:
library(lme4) #load package
ce_model <- lmer(RST_AST ~ domain + ce+ (1|subj), data = df_RST_AST)然后我得到错误消息:
boundary (singular) fit: see ?isSingular我不知道怎么回事?
有人能帮我吗?
发布于 2022-02-02 03:29:34
很可能,您的一个变量没有被评估为造成差异的原因。最有可能的是,“罪魁祸首”是随机变量(在您的情况下是子j)。如果您查看模型的输出(只需打印ce_model并点击enter,或汇总(Ce_model)),就应该有一个变量具有零或非常接近零的方差和标准差。有几种方法来处理这个问题(请参阅上面的任何链接)。我建议分析错误的变量,它有什么样的分布,它可以被转换,...etc。也许您需要对数据执行一个或多个广义模型。查看您的模型也可能有帮助;但是,如果没有更多的信息,则很难判断。
https://stackoverflow.com/questions/70880851
复制相似问题