You need to enable JavaScript to run this app.
优惠活动
大模型
产品
解决方案
定价
更多

如何为glmmLasso输入模型规格获预期迭代?遇聚类变量警告

使用glmmLasso时随机效应分组变量的警告与NA结果问题

我在使用glmmLasso分析数据时,将subj(受试者ID)设为随机效应的分组变量,该用法参考了glmmLasso的公开示例及CRAN包自带的足球数据集,但运行代码时始终弹出如下警告:

Warning message:
In est.glmmLasso.RE(fix = fix, rnd = rnd, data = data, lambda = lambda, :
Cluster variable should be specified as a factor variable!

尽管除subj外我并未使用其他分类变量。

我的示例代码如下:

library(lme4)
library(tidyverse)
library(tools)
library(glmmLasso)

ChoiceData <- read_csv("FatigueData.csv")

dfs <- ChoiceData
for(i in c(2,3,5)){
  dfs[,i] <- scale(ChoiceData[,i])
}

LassodModel = glmmLasso(C ~ E_sq + R + trialNum, 
                        rnd = list(subj =~ 1),
                        data = na.omit(dfs), lambda=5, family = binomial())

summary(LassodModel)

我尝试在随机效应中对subj使用as.factor(.),却出现语法错误:

LassodModel = glmmLasso(C ~ E_sq + R + trialNum, rnd = list(as.factor(subj) =~ 1),
Error: unexpected '=' in:
"LassodModel = glmmLasso(C ~ E_sq + R + trialNum, 
                        rnd = list(as.factor(subj) =" 

调整lambda值也没有任何变化。设置rnd=NULL时代码可运行,但未执行任何迭代,系数的StdErr、z.value及p.value均为NA:

Fixed Effects:

Coefficients:
Estimate StdErr z.value p.value
(Intercept) -0.338451 NA NA NA
E_sq -1.173294 NA NA NA
R 1.291759 NA NA NA
trialNum -0.043156 NA NA NA

No random effects included!

运行glmmLasso包内的demo也得到类似结果,已重装更新包,仍无法解决问题。

内容的提问来源于stack exchange,提问作者trytryagain404

相关产品推荐
方舟 Agent Plan

超全模态模型 × Harness 升级,最新支持 Deepseek-V4.1-Flash、GLM-5.3 系列、Doubao-Seedream-5.0-pro、Kimi-K3 (部分), 限时 9.9 元起

最近更新时间:2026.07.29 13:07:41