问题标签 [lme4]

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

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

我尝试使用bfFixefLMER_t.fnc或fitLMER.fnc从LMERConvenienceFunctions包中。在这两种情况下，我都会收到“输入模型不是 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

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

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

我有一个数据集 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/ 如何解决它（我尝试重新安装旧版本但我在包之间存在兼容性问题）的任何想法？

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

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

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

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

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

原始示例

这是错误消息

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

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

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

对于我的完整模型，我将命令设置如下：

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

谢谢！

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

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

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

以下是一些数据（称为 dat）：

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

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

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

非常感谢！