glmmLasso 错误和警告

Question

我正在尝试使用 glmmLasso 在广义线性混合模型中执行变量选择，但出现了我无法解决的错误和警告。数据集不平衡，一些参与者（PTNO）的样本比其他参与者多；没有丢失的数据。我的因变量是二进制的，所有其他变量（除了 ID 变量 PTNO）都是连续的。 我怀疑正在发生一些非常普遍的事情，但显然看不到它，并且在文档或网络上没有找到任何解决方案。 该代码基本上只是改编自 glmmLasso 足球示例：

glm8 <- glmmLasso(Group~NDUFV2_dCTABL+GPER1_dCTABL+ ESR1_dCTABL+ESR2_dCTABL+KLF12_dCTABL+SP4_dCTABL+SP1_dCTABL+  PGAM1_dCTABL+ANK3_dCTABL+RASGRP1_dCTABL+AKT1_dCTABL+NUDT1_dCTABL+                   POLG_dCTABL+   ADARB1_dCTABL+OGG_dCTABL+ PDE4B_dCTABL+  GSK3B_dCTABL+ APOE_dCTABL+  MAPK6_dCTABL, rnd = list(PTNO=~1),  
    family = poisson(link = log), data = stackdata, lambda=100, 
    control = list(print.iter=TRUE,start=c(1,rep(0,29)),q.start=0.7)) 

错误信息如下所示。具体来说，我不相信数据集中有任何 NA，我不确定有关因子变量的警告的含义。

迭代 1 grad.lasso[b.is.0]

包含必要变量的缩写数据集以 R 格式提供，可以下载 here。 我希望可以指导我如何继续进行分析。如果数据集有任何问题或您无法下载，请告诉我。非常感谢任何帮助。

Answer 1

只是为了跟进上面的@Kristofersen 评论。确实是 start 向量搞砸了您的分析。

如果我跑

glm8 <- glmmLasso(Group~NDUFV2_dCTABL+GPER1_dCTABL+ ESR1_dCTABL+ESR2_dCTABL+KLF12_dCTABL+SP4_dCTABL+SP1_dCTABL+  PGAM1_dCTABL+ANK3_dCTABL+RASGRP1_dCTABL+AKT1_dCTABL+NUDT1_dCTABL+                   POLG_dCTABL+   ADARB1_dCTABL+OGG_dCTABL+ PDE4B_dCTABL+  GSK3B_dCTABL+ APOE_dCTABL+  MAPK6_dCTABL, 
                  rnd = list(PTNO=~1), 
                  family = binomial(), 
                  data = stackdata, 
                  lambda=100,     
                  control = list(print.iter=TRUE))

然后一切都很好而且很花哨（即，它会收敛并产生一个解决方案）。您已经复制了带有泊松回归的示例，您需要根据您的情况调整代码。我不知道输出是否有意义。

快速说明：我在上面的代码中使用二项式分布，因为您的结果是二元的。如果估计相对风险是有意义的，那么泊松可能是合理的（并且它也会收敛），但是您需要重新编码您的结果，因为这两个组被定义为1 和2，这肯定会弄乱泊松回归.

换句话说，做一个

stackdata$Group <- stackdata$Group-1

在您运行分析之前。