Tag: lmer

用先前的估计值重新开始混合效应模型估计

我使用包lme4 lmer()来估计混合效果模型。 这个效果很好，但是现在我想运行估计过程进行固定次数的迭代，然后通过指定最后估计过程计算的开始值来恢复过程。 根据?lmer的帮助，这是可能的，通过设定参数： start – 这些是新的起始值，根据帮助，可以从一个拟合模型中提取ST的值，并使用这些值，即使用x@ST maxiter – 作为control的命名参数提供 所以，举个例子，假设我想用iris数据来适应lme ，可以试试这个： library(lme4) # Fit model with limited number of iterations frm <- "Sepal.Length ~ Sepal.Width | Species" x <- lmer(frm, data=iris, verbose=TRUE, control=list(maxIter=1), model=FALSE) # Capture starting values for next set of iterations start <- list(ST=x@ST) # Update model twoStep <- lmer(frm, data=iris, verbose=TRUE, […]