如何在 lme 中分别适配两个随机效应？
我正在使用nlme包进行REML的线性混合效应模型拟合。以下是适用于我的代码： # Linear mixed-effects model fit by REML (intercept and not slope) x <- lme (DV ~ IV1 + IV2 + IV1*IV2, ...
在lme模型拟合中提取每个单位的系数及其标准误差
在类似于以下的线性混合模型中，我如何提取每个实验单位（plot）的系数（b0和b1）及其相应的标准误差： Better fits for a linear model 使用与此数据集（df）相同的数据，以及拟合的模型（fitL1）：我如何获得像这样的数据帧... plot b...
如何在nlme和lme4中指定不同的随机效应？
我希望使用nlme::lme在模型中指定不同的随机效应（底部有数据）。这些随机效应是：1）intercept和position在subject上变化；2）intercept在comparison上变化。使用lme4::lmer很容易实现： lmer(rating ~ 1 + position...
使用ggplot绘制非线性回归列表
作为非线性回归分析的输出图表，来自此链接 https://stats.stackexchange.com/questions/209087/non-linear-regression-mixed-model 使用此数据集： zz <-(" iso temp diam Itiqui...
使用交叉随机效应的nlme
我正在尝试使用nlme包拟合一个交叉非线性随机效应模型，就像这个问题和这篇邮件列表帖子中提到的那样，将其作为线性随机效应模型。但无论我尝试什么，都会出现错误。以下是一个示例： library(nlme) ##### # simulate data set.seed(18112003) na...