问题标签 [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 投票
1 回答
1443 浏览

r - R中的频率权重(使用lme4的多级)

我正在使用 lme4 构建一个多级模型(两个级别:具有教育数据的班级和学生级别)。我正在使用 ICCS 2009 数据:http ://www.iea.nl/iccs_2009.html )。ICCS 2009 建议使用权重。我做了一个“校内学生体重”作为班级和学生水平权重因子的乘积。我还使权重的总和等于样本量。我的问题:如何在 R 中指定我想将此权重用作频率权重。我必须使用“调查”包吗?

0 投票
0 回答
1069 浏览

r - R- represent glmm (Package MuMin and lme4) with categorical and numeric variable

I have a table with Ancylostoma's infection, vs sex (2 factor), location (2 factor), year, management (2 factor), ancestry (4 factor) and viremia like categorical variable, and the I have HL an age like numeric variable.**

I made a glmm:

And I got this result

I would like represent the probability of infection (y) vs age (x), but using the estimate of my model.avg**

I tried with this script:

But I haven't got graph...or I have error

Someone know any form to represent my model.avg with categorical and variable numeric?, But taking into account that I only want represent probability of infection vs age, with two line: localizacion1 and localizacion2...(localization had 2 factors).**

my original date would be this table:

#

#

And with that date I make my model's estimates. Then I would like do prediction

Thank you, I don't sure if I am explaining well

I so sorry for my english

0 投票
1 回答
959 浏览

r - 后向和前向拟合函数的 LMERConvenienceFunctions 错误:模型不是 mer 对象

我尝试使用bfFixefLMER_t.fncfitLMER.fncLMERConvenienceFunctions包中。在这两种情况下,我都会收到“输入模型不是 mer 对象”的错误。

我尝试了http://artax.karlin.mff.cuni.cz/r-help/library/LMERConvenienceFunctions/html/00Index.html中的示例。我得到同样的错误。

例如,当我从示例中运行时

这是我得到的结果。

有人对这些功能有过这种经验吗?

有一些函数可以反向拟合固定效应和正向拟合随机效应。有没有办法对模型的固定效应进行前向拟合glmer?还是这在统计上毫无意义?我正在做生态建模,所以我对高级统计的理解并不多,所以,如果有人可以用外行的话更好地解释一下

R 版本 3.0.1 (2013-05-16) 平台:x86_64-pc-linux-gnu (64-bit)

语言环境:[1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_GB.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_GB.UTF-8
[ 7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C

附加的基础包:[1] stats graphics grDevices utils datasets methods base

其他附加包:[1] LMERConvenienceFunctions_2.0 lme4_0.99999911-8
[3] RcppEigen_0.3.1.2.1 Rcpp_0.10.4
[5] Matrix_1.0-12 lattice_0.20-23
[7] LCFdata_1.0

通过命名空间加载(未附加):[1] grid_3.0.1 MASS_7.3-28 minqa_1.2.1 nlme_3.1-111 rpart_4.1-2
[6] splines_3.0.1 tools_3.0.1

0 投票
1 回答
989 浏览

r - 如何修改插槽 lme4 >1.0

我一直在使用下面的代码成功地修改了模型的“Zt”、“L”和“A”插槽,使用 lme4 版本 <1.0。今天刚更新到lme4 1.0-4,发现model对象不一样。任何人都可以提供有关如何修改新 lmer 模型对象中的这些插槽的见解/指导吗?

任何有关如何修改这些插槽的帮助或见解将不胜感激。与此同时,我想我会坚持为这些模型使用旧版本的 lme4。

0 投票
3 回答
3293 浏览

r - lme4 和 languageR 兼容性错误:“输入模型不是 mer 对象”

我有一个数据集 dat2,我想在其上拟合线性混合效应模型。我过去使用 lmer()(包 lme4)作为 pvals.fnc 的补充来计算相关的 p 值。

但是,由于我使用新的 lme4 (1.0.4) 和 languageR (1.4) 软件包重新安装了 R 3.0.2 版本,因此我收到有关 lmer 函数输出的错误。它表示输出不是 mer 对象。事实上,它的类是 lmeRmod。

这是我使用的代码:

我得到这个错误:

事实上,当我查看 mod1 时,我发现它是一个 lmeRmod 对象,而不是一个 mer 对象。

这个问题与这里描述的问题非常相似。关于 1/ 可能是什么问题(为什么我不输出 mer 对象)和 2/ 如何解决它(我尝试重新安装旧版本但我在包之间存在兼容性问题)的任何想法?

任何帮助都会很棒!谢谢 !

0 投票
0 回答
160 浏览

r - 此消息是什么意思:“isS4”不是 BUILTIN 函数?

我安装了lme4软件包并尝试运行example(glmer)但收到错误消息:

但是,当我运行时builtins()[590],R 返回"isS4",这让我相信,无论 BUILTIN 函数是什么,都isS4应该是一个。我在网上搜索了错误消息,并在 LinkedIn 上发现了一个关于它的问题,这次是来自试图安装 R Commander 的人,但没有答案。我注意到我们都在使用 R 2.15.2 并且想知道这个问题是否与不同版本的 R 包的兼容性有关。

我真的很想知道是否有办法解决这个错误,以及究竟是什么原因造成的。

0 投票
1 回答
895 浏览

r - 由于使用了新版本的 lme4,从以前的答案复制结果不起作用

我试图从这个问题的答案中重现结果“Estimating random effects and apply user defined correlation/covariance structure with R lme4 or nlme package” https://stats.stackexchange.com/questions/18563/estimating-random-效果和应用用户定义的相关协方差结构

亚伦伦达尔的代码

原始示例

这是错误消息

我将不胜感激任何帮助?谢谢

0 投票
1 回答
12054 浏览

r - 使用 glmer 处理嵌套数据

我有跨纬度感染特定宿主物种的病原体多样性的数据。该设计涉及在不同纬度的 4 个地点的 3 个地点收集 20 个人,因此我有 20 个人,嵌套在 3 个地点,嵌套在 4 个地点。

鉴于我的病原体多样性数据是具有许多零的计数数据,这就是为什么我一直在探索使用 GLMM 和lme4::glmerR 中的命令来分析数据的原因。对于分析,我想将纬度视为数字固定因子,将站点视为与位置嵌套的随机因子。

对于我的完整模型,我将命令设置如下:

这是我所描述的正确语法吗?

谢谢!

0 投票
1 回答
3696 浏览

r - lme4 升级产生错误消息 Error in `[[<-.data.frame`(`*tmp*`, i, value = integer(0))

我上周用 lmer (lme4) 构建了一个广义混合线性模型,效果很好:

然而,更新 lme4 包后,模型不再工作,并给出错误消息:

使用 glmer 而不是 lmer 并没有解决第一条错误消息。任何建议将不胜感激!

以下是一些数据(称为 dat):

提前谢谢了!图斯·范诺德韦克

0 投票
1 回答
21035 浏览

r - 增加新版本 lmer 的迭代次数?

我刚刚将 lme4 更新到 1.0-4 版本,当我运行 lmer() 之前正在收敛的混合效果模型时,现在会打印以下警告:

所以,我想尝试增加迭代次数,看看我是否可以解决这个问题。(我必须说我不知道​​是什么导致了警告,因为消息的第一部分听起来有点不透明)。无论如何,我在文档中读到现在我应该使用lmerControl(),但我无法实现它。有人可以给我一个具体的例子来说明你如何做到具体吗?(帮助文件没有帮助)。这是我的模型:

非常感谢!