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