问题标签 [nlme]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
251 浏览

r - texreg 与 lmer 和 lme 对象 - 差异不同

伙计们,我正在尝试打印包含差异的 lmer() {nlme} 和 lme() {lme4} 模型的 {texreg} 表。然而,两个模型之间的差异显着不同(几个数量级)。似乎 lme() 方差是 lmer() 方差的平方根。哪些是正确的?

0 投票
0 回答
256 浏览

r - R:使用 nLme 和 multcomp 分解交互

我试图了解混合设计研究的交互。

variable是重复测量并且status是主体间因子

有什么不同吗?:

发现任务有所不同,所以我想我会看看一些主要影响事后:

发现了一个重要的交互,那么为什么不看一些交互事后:

问题:

那么,有没有人需要查看使用这些工具的交互?

0 投票
1 回答
834 浏览

r - 绘制 lmList 拟合中系数的各个置信区间

我正在尝试使用 R 解决 Pinheiro 和 Bates 书中的练习(S 和 S-PLUS 中的混​​合效果模型)。它使用来自 nlme 包的数据像素。

练习说:使用 lmList 在一天中为每只狗拟合一个单独的二次模型。绘制 lmList 拟合中系数的各个置信区间。验证只有截距和线性系数似乎随 Dog 显着变化。

我安装了这样的模型:

但是当我尝试使用此函数绘制置信区间时:

R 给出了这个错误:

有没有人知道如何解决任何问题?

0 投票
0 回答
324 浏览

r - 用 visreg() 说明反向转换的 GLS

我有一个带有对数转换因变量的 GLS。使用 visreg() 进行可视化会给出错误代码:

有趣的是,同样适用于 LM:

0 投票
0 回答
648 浏览

r - 警告消息:1:在 logLik.reStruct(object, conLin)

我正在 lme对模拟数据的函数运行多元线性混合效应模型:

只是为了清除一些变量: cs_d0 和 va_d0 是虚拟变量。我使用了一个 melt 函数来组合数据集。
值具有 cs_d0 和 va_d0 的值。

我面临的问题是:

我不知道这个警告信息是什么意思!!!另外,我认为结果是不明智的

因此,在固定效果中:cs_d0 和 va_do 的截距相同,斜率也相同!!!只是我在这里问一下,如果有人可以帮助我,好吗?

0 投票
0 回答
58 浏览

r - 减少模拟数据的变化

我正在对双变量纵向数据进行模拟。我将描述我的代码,我希望减少第 0 组csse的变化。因此,有两个类Y = (0,1),并且对于每个类,我生成线性混合效果模型。每组有 4 个随机参数(截距和斜率),也就是说协方差矩阵是 4x4,并且随机效应之间存在相关性。此外,残差是矩阵 2x2(每个变量都有误差)并且存在相关性。我假设每组中的每个变量都有 4 个固定效应参数。
因此,我的代码是:

总结如下:

只关注最后的摘要,我的问题是如何减少Y==0不包括 group 1 范围的 group 0 上 csse 的变化。请问有人有答案吗?谢谢

0 投票
1 回答
862 浏览

r - 函数 y = at^b * exp(-ct) 有名称吗?可以线性化吗?我如何估计a,b,c?

我正在尝试拟合非线性模型,但在网上找不到任何好的示例。

这个函数有名字吗?

可以线性化吗?

我试图估计参数a,bc随机效应g(如在组中)作为 time 的函数t,如下所示。nls我可以在没有随机效应的情况下拟合模型,但无法让模型收敛。欢迎提出建议(最好在 R 中,但任何合适的包都可以)?

0 投票
1 回答
5375 浏览

r - 混合建模 - lme 和 lmer 函数之间的不同结果

我目前正在阅读 Andy Field 的书《使用 R 发现统计数据》。第 14 章是关于混合建模的,他使用lme了包中的函数nlme

他使用快速约会数据创建的模型是这样的:

我尝试使用包中的lmer函数重新创建一个类似的模型lme4;但是,我的结果是不同的。我以为我有正确的语法,但也许不是?

此外,当我运行这些模型的系数时,我注意到它只会为每个参与者产生随机截距。然后我试图创建一个同时产生随机截距和斜率的模型。我似乎无法让这两个函数的语法正确。任何帮助将不胜感激。

0 投票
1 回答
1247 浏览

r - 用 gls 预测

我正在学习如何使用 gls 进行预测,但是当我尝试这样做时,我不断收到错误消息。这是我正在使用的代码。

到目前为止,这工作正常。现在,我想预测特定日期特定风速的臭氧。所以我用。

最后一行代码给了我以下错误:

eval 中的错误(expr、envir、enclos):找不到对象“臭氧”

然后我尝试了以下代码来预测:

这是我得到的错误,尽管该nlme库已包含在内。

错误:找不到函数“predict.gls”

谁能告诉为什么会出现这些错误?另外,在这里使用什么代码进行预测的正确方法是什么?

0 投票
1 回答
211 浏览

r - nlme 和 lme4 中的不同随机效应

我适合这个模型nlme

我可以得到系数

现在我适合相同的模型lme4

我得到完全相同的固定效果:

我现在想要每个工人的随机效果,它们不会显示在摘要中,但它必须是这样的:

为什么它们与nlme结果不同,而固定效果相同?