问题标签 [lme4]

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 投票
3 回答
5020 浏览

r - 使用 Lattice(或其他东西)在 R 中绘制 lme4 的回归结果

由于先前的答案,我已经使用 lme4 进行了回归。现在我对每个州都有一个回归拟合,我想使用 lattice 为每个州绘制 QQ 图。我还想以格子格式绘制每个状态的误差图。如何使用 lme4 回归的结果制作格子图?

下面是一个使用两种状态的简单示例(是的,我喜欢一个很好的头韵)。我想制作一个由对象制成的两个面板格子适合。

0 投票
1 回答
1240 浏览

r - LanguageR/lme4:对具有随机相关参数的模型的 p 值有任何帮助吗?

我将包languageR用于混合效果模型,其语法位于本文末尾。我可以pvals.fnc用来获取模型 1 和 3(和)的p值。将此与模型二一起使用会给出以下错误消息:hd_lmer1hd_lmer2

p2 = pvals.fnc(hd_lmer2) pvals.fnc(hd_lmer2) 中的错误:对于具有随机相关参数的模型,MCMC 采样尚未在 lme4_0.999375 中实现

如果有人能帮助我了解如何获得此类模型的 p 值,我将不胜感激。

楷模:

0 投票
1 回答
1697 浏览

r - 将 sim() 与 lmer() 一起使用

我已经使用相同的预测变量运行了两个多级逻辑回归,但有两个不同的响应:

Income是离散的c(-2,-1,0,1,2)HDI在 0 和 1 之间是连续的,Country是分类的,并且响应都是 1/0。

为了绘制置信带,我使用 arm 包中的 sim() 函数运行模拟:

第一个计算得很好。第二个返回以下错误:

我实际上是用 8 种不同的反应来做这件事的。其中六个工作正常,其中两个返回此错误。

有谁知道如何纠正这个?

0 投票
1 回答
1685 浏览

r - 从 lmList 对象中提取 RMSE

我有以下代码从 lmList 对象中提取系数。

此代码工作正常。但是,这是提取系数的有效方法吗?

我的主要问题是,如何从同一个 lm.model 对象中提取 RMSE?

0 投票
3 回答
7313 浏览

r - 有什么方法可以从 lme4 mer 模型拟合对象生成 LaTeX 表?

有谁知道一种从 lme4mer对象生成高质量出版质量的 LaTeX 表的方法?xtable方法(包xtable)和latex方法(包)都不Hmisc知道如何处理mer对象。

例如,给定这种拟合:

是否有任何选项可以生成固定和随机效应的系数估计值的漂亮 LaTeX 表?

编辑:

因为这有点隐藏在下面的评论线程中,请注意社区 wiki 正在为 R LaTeX 表开发:Tools for making latex tables in R

0 投票
3 回答
13437 浏览

r - lme4中的多元线性混合模型

我想知道如何将多元线性混合模型与lme4. 我用以下代码拟合了单变量线性混合模型:

我想知道如何将多元线性混合模型与lme4. 数据如下:

提前感谢您的时间和合作。

0 投票
1 回答
826 浏览

r - 随机模型公式对象

我想将公式放入随机模型中,但我认为以下错误是由于错误的公式对象(?),但无法修复。

错误:length(formula <- as.formula(formula)) == 3 is not TRUE

0 投票
2 回答
3701 浏览

r - 在 nlme 或 lme4 中获得固定效应的广义最小二乘均值

最小二乘均值及其aov对象的标准误差可以通过model.tables函数获得:

nlme我想知道如何从或lme4对象的标准误差中获得广义最小二乘均值:

任何评论和提示将不胜感激。谢谢

0 投票
2 回答
1081 浏览

r - lme4 和 nlme 包之间的干扰

我有这个脚本:

两周前有效。今天我得到

我觉得是lme4和nlme之间的干扰,但我不知道如何删除nlme。如果有任何帮助,我试过这个:

但这对我来说毫无意义,除了加强我对两个包裹之间的干扰感。有任何想法吗?我已经退出 R,重新开始,加载包,更新包,我能想到的一切。删除它们基本上不起作用。

0 投票
1 回答
145 浏览

r - 为什么 log() 语句会导致 lmList 失败?

在解决这个问题时:lme4 和 nlme 之间的干扰将 OP 打包,我确定将响应向量包装在对导致失败log()的调用中。lme4::lmListlmList

例如,

给出回应

但是,如果您首先记录响应和预测器,它就可以正常工作,即

有人可以解释为什么在第一个示例中将响应包装在 log() 中不起作用吗?在我看来(至少直觉上)它应该这样做。