R中的交叉重复效应和AR1协方差结构的线性混合模型
我拥有来自参与者（part）的受试者内生理数据，他们都在三轮（round）中查看刺激物（阅读报纸），每个轮次有五份报纸（paper），每个报纸中都有不同数量的访问（visit）。我有两个固定因素（CONDhier和CONDSbund）加上交互作用来预测生理状态（例如EDA），通常是自回归的。我...
在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,...
如何使用lme4比较没有随机效应的模型和有随机效应的模型？
我可以使用nlme包中的gls()函数来建立没有随机效应的mod1模型。 然后，我可以使用lme()函数来建立包含随机效应的mod2模型，并通过AIC对比mod1和mod2。 mod1 = gls(response ~ fixed1 + fixed2, method="REML", data...
如何从nlme调用中获取Hessian
library(nlme) fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), data = Loblolly, fixed = Asym + R0 + lrc ~ 1, ...
在nlme和lme4中拟合相同的模型
数据来自这里。 library(nlme) dat0 <- read.table("aids.dat2",head=T) dat1 <- dat0[dat0$day<=90, ] # use only first 90-day data dat2 <- dat1[...
在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, ...
允许GLS中的相关参数依赖于分组因素
下面是我遇到的问题的最小工作示例。 我正在使用nlme包中的Orthodont数据集，该数据集由27个儿童（16个女孩，11个男孩）的4个测量组成。 为了建立相关性模型，我使用未结构化协方差结构，通过指定correlation = corSymm（form = ~1 | Subject）来实现...
汽车包无法加载，nlme版本错误。
当我尝试加载“car”软件包时，出现以下错误： library(car) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : namespace 'nlm...
非线性回归与随机效应以及lsoda
我遇到了一个无法解决的问题。我想使用nlme或nlmODE来执行具有随机效应的非线性回归，以固定系数的二阶微分方程的解作为模型（阻尼振荡器）。 我成功地使用nlme进行简单模型的建立，但似乎使用deSolve生成微分方程的解会导致问题。以下是一个示例和我面临的问题。 数据和函数 这是使用...