R:找不到lme4 4的nlmer错误变量。

Question

以下示例中的数据来自这里

library(tidyverse)
library(lme4)
dat <- read.table("aids.dat2",head=T) %>%
  filter(day <= 90) %>%
  mutate(log10copy = log10(lgcopy)) %>%
  na.omit()

> head(dat)
  patid day    cd4   lgcopy    cd8 log10copy
2 11542   2 159.84 4.361728 619.38 0.6396586
3 11542   7 210.60 3.531479 666.90 0.5479566
4 11542  16 204.12 2.977724 635.04 0.4738844
5 11542  29 172.48 2.643453 407.68 0.4221716
6 11542  57 270.94 2.113943 755.78 0.3250933
8 11960   2 324.72 3.380211 856.08 0.5289438

运行以下代码会给出错误：Error in eval(expr, envir, enclos) : object 'log10copy' not found，但log10copy显然是数据集中的列之一？

lme4.fit <- lme4::nlmer(log10copy ~ exp(p1-b1*day) + exp(p2-b2*day + 1) + 
                          (1|p1) + (1|b1) + (1|p2) + (1|b2), data = dat)

我想要对p1、b1、p2、b2和4种随机效应在相同的参数集上拟合一个模型。

Answer 1

你这里有几个问题..。

1)起始值必须是命名向量。

2) nlmer中的data参数应该接收dat作为值，而不是像您的示例那样接收aids.dat

start <- c(p1 = 10, b1 = 0.5, p2 = 6, b2 = 0.005)
lme4.fit <- lme4::nlmer(log10copy ~ exp(p1-b1*day) + exp(p2-b2*day + 1) ~ 
                            (p1|patid) + (b1|patid) + (p2|patid) + (b2|patid), data = dat,
                        start = start)

这将触发以下错误：

Erreur : is.matrix(gr <- attr(val, "gradient")) is not TRUE

如文件所述：

目前，农林(.)公式部分不仅必须返回数字向量，而且还必须具有“梯度”属性，即矩阵。函数SSbiexp、SSlogis等，请参阅selfStart，提供此功能(以及更多)。或者，您也可以使用deriv()自动生成这样的函数或表达式。

然后，您可以修改文档提供的示例：

## a. Define formula
nform <- ~ exp(p1-b1*input) + exp(p2-b2*input + 1)
## b. Use deriv() to construct function:
nfun <- deriv(nform, namevec=c("p1", "b1", "p2", "b2"),
              function.arg=c("input","p1", "b1", "p2", "b2"))
lme4.fit <- lme4::nlmer(log10copy ~ nfun(day, p1, b1, p2, b2) ~ 
                            (p1|patid) + (b1|patid) + (p2|patid) + (b2|patid), data = dat,
                        start = start)

然后，您将出现以下错误

Error in fn(nM$xeval()) : prss failed to converge in 300 iterations

这可能意味着你的模型对你的数据来说太复杂了…或者我在规范中犯了一个错误，因为我对nlmer不太了解(我只是尝试应用文档.)我也不知道你的模型/问题。

当你改变优化器时，收敛问题似乎就消失了.

有关使用这里进行“故障排除”(包括收敛问题)的建议，请参见lme4

lme4.fit <- lme4::nlmer(log10copy ~ nfun(day, p1, b1, p2, b2) ~ 
                            (p1|patid) + (b1|patid) + 
                            (p2|patid) + (b2|patid), 
                        data = dat,
                        start = start, 
                        nlmerControl(optimizer = "bobyqa"))