例子?如何在lmer中指定θ和β的起始值?

例子?如何在lmer中指定θ和β的起始值?,r,lme4,R,Lme4,我对R相对来说是新手,仍然在努力理解一些文档。我正在使用glmer分析,我的目标是指定θ(截距)和β(固定效果)的起始值 任何人都可以提供一个示例,说明一个beta和多个beta的代码应该是什么样子的吗 我已经看过了glmer文档,上面写着“开始 模型中参数的起始值的命名列表或数字向量。数字起始参数将用作θ的起始值。如果start是列表,则θ元素(数字向量)将用作第一个优化步骤的起始值(对于较低Cholesky系数的对角线元素,默认值为1;对于非对角线元素,默认值为0);从第一步开始的θ的拟合值

我对R相对来说是新手,仍然在努力理解一些文档。我正在使用glmer分析,我的目标是指定θ(截距)和β(固定效果)的起始值

任何人都可以提供一个示例,说明一个beta和多个beta的代码应该是什么样子的吗

我已经看过了glmer文档,上面写着“开始
模型中参数的起始值的命名列表或数字向量。数字起始参数将用作θ的起始值。如果start是列表,则θ元素(数字向量)将用作第一个优化步骤的起始值(对于较低Cholesky系数的对角线元素,默认值为1;对于非对角线元素,默认值为0);从第一步开始的θ的拟合值加上开始[[“fixef”]],用作第二个优化步骤的起始值。如果start同时包含fixef和θ元素,则跳过第一个优化步骤。有关优化的更多详细信息或更精细的控制,请参阅。”

我已经成功地指定了θ开始值,但我不知道如何同时指定θ和β(或多个β)


Glucose1我从未使用过glmer,但在非线性回归中,代码是:

Richards = nls(Hd ~ a/(1+b*exp(-c*Id)), data = Capacidade_Prod, start = list(a = 20, b = 4, c = 0.01))

下面是一个
glmer
示例,其中添加了start参数,听起来您希望使用它:

library(lme4)

glmer(
  cbind(incidence, size - incidence) ~ period + (1 | herd),
  data = cbpp,
  family = binomial,
  start = list(theta = 0.6, fixef = c(-1.4, -1, -1.1, -1.6))
)
因此,在你的情况下:

Glucose1b <-
  glmer(
    GlucoseYN ~ 1 + (1 | Subject),
    data = LongFormat,
    family = "binomial",
    control = glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 10000000)), 
    start = list(theta = 3.276187, fixef = -0.3515557)
  )
Glucose1b是密切相关的,放在这里作为参考。
Glucose1b <-
  glmer(
    GlucoseYN ~ 1 + (1 | Subject),
    data = LongFormat,
    family = "binomial",
    control = glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 10000000)), 
    start = list(theta = 3.276187, fixef = -0.3515557)
  )