如何在 lme 中分别适配两个随机效应？
我正在使用nlme包进行REML的线性混合效应模型拟合。以下是适用于我的代码： # Linear mixed-effects model fit by REML (intercept and not slope) x <- lme (DV ~ IV1 + IV2 + IV1*IV2, ...
非线性回归与随机效应以及lsoda
我遇到了一个无法解决的问题。我想使用nlme或nlmODE来执行具有随机效应的非线性回归，以固定系数的二阶微分方程的解作为模型（阻尼振荡器）。 我成功地使用nlme进行简单模型的建立，但似乎使用deSolve生成微分方程的解会导致问题。以下是一个示例和我面临的问题。 数据和函数 这是使用...
如何从nlme调用中获取Hessian
library(nlme) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, ...
允许GLS中的相关参数依赖于分组因素
下面是我遇到的问题的最小工作示例。 我正在使用nlme包中的Orthodont数据集，该数据集由27个儿童（16个女孩，11个男孩）的4个测量组成。 为了建立相关性模型，我使用未结构化协方差结构，通过指定correlation = corSymm（form = ~1 | Subject）来实现...
在nlme和lme4中拟合相同的模型
数据来自这里。 library(nlme) dat0 <- read.table("aids.dat2",head=T) dat1 <- dat0[dat0$day<=90, ] # use only first 90-day data dat2 <- dat1[...
在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和lme4中指定不同的随机效应？
我希望使用nlme::lme在模型中指定不同的随机效应（底部有数据）。这些随机效应是：1）intercept和position在subject上变化；2）intercept在comparison上变化。使用lme4::lmer很容易实现： lmer(rating ~ 1 + position...
在R中，lmer（或lme）中的随机效应变量是否自动被视为因子？
我了解在混合效应模型中，将连续或数值变量作为随机效应并不太合理（例如，请参见此处）。 但我想知道的是，R中的lme4::lmer或nlme::lme是否有意防止你这样做... 具体来说，我的问题是：如果我将任何非因子（非分类）变量作为随机效应提供给lmer（或lme），该函数会自动将其视为...
解决错误信息“在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, ...
在lme模型拟合中提取每个单位的系数及其标准误差
在类似于以下的线性混合模型中，我如何提取每个实验单位（plot）的系数（b0和b1）及其相应的标准误差： Better fits for a linear model 使用与此数据集（df）相同的数据，以及拟合的模型（fitL1）：我如何获得像这样的数据帧... plot b...