在lme模型拟合中提取每个单位的系数及其标准误差
在类似于以下的线性混合模型中，我如何提取每个实验单位（plot）的系数（b0和b1）及其相应的标准误差： Better fits for a linear model 使用与此数据集（df）相同的数据，以及拟合的模型（fitL1）：我如何获得像这样的数据帧... plot b...
使用ggplot绘制非线性回归列表
作为非线性回归分析的输出图表，来自此链接 https://stats.stackexchange.com/questions/209087/non-linear-regression-mixed-model 使用此数据集： zz <-(" iso temp diam Itiqui...
如何在nlme和lme4中指定不同的随机效应？
我希望使用nlme::lme在模型中指定不同的随机效应（底部有数据）。这些随机效应是：1）intercept和position在subject上变化；2）intercept在comparison上变化。使用lme4::lmer很容易实现： lmer(rating ~ 1 + position...
解决错误信息“在NLS步骤中，步长减半因子降至最低值”：调整nlsTols无效。
我将使用gnls函数来拟合一些数据的逻辑增长曲线。 数据： structure(list(Nest = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, ...
在nlme和lme4中拟合相同的模型
数据来自这里。 library(nlme) dat0 <- read.table("aids.dat2",head=T) dat1 <- dat0[dat0$day<=90, ] # use only first 90-day data dat2 <- dat1[...
汽车包无法加载，nlme版本错误。
当我尝试加载“car”软件包时，出现以下错误： library(car) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : namespace 'nlm...
在R中，lmer（或lme）中的随机效应变量是否自动被视为因子？
我了解在混合效应模型中，将连续或数值变量作为随机效应并不太合理（例如，请参见此处）。 但我想知道的是，R中的lme4::lmer或nlme::lme是否有意防止你这样做... 具体来说，我的问题是：如果我将任何非因子（非分类）变量作为随机效应提供给lmer（或lme），该函数会自动将其视为...
使用交叉随机效应的nlme
我正在尝试使用nlme包拟合一个交叉非线性随机效应模型，就像这个问题和这篇邮件列表帖子中提到的那样，将其作为线性随机效应模型。但无论我尝试什么，都会出现错误。以下是一个示例： library(nlme) ##### # simulate data set.seed(18112003) na...
在na.fail.default中出现错误：对象中缺少值 - 但没有缺少的值。
我正在尝试使用这些数据运行一个线性混合效应模型：tot_nochc=runif(10,1,15) cor_partner=factor(c(1,1,0,1,0,0,0,0,1,0)) age=runif(10,18,75) agecu=age^3 day=factor(c(1,2,2,3,3,...
如何从nlme调用中获取Hessian
library(nlme) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, ...