我经常通过将我需要的部分粘贴在一起来为lm或lme等模型拟合函数指定公式参数，就像@DWin对这个问题的回答：。(as.formula(f), random=~1|id, data=d, method="ML")m1 <- getm("x1")但是，使用nlme包中的lme，比较以<e
浏览 1提问于2011-10-06得票数 9
回答已采纳
2回答
我正在拟合一个关于缺席-在场数据的模型，我想检查随机因素是否显着。要做到这一点，应该将glmm与glm进行比较，并检查LR-test，如果我理解正确的话，哪一个是最重要的。但是如果我执行ANOVA(glm，glmm)，我得到的是Deviance Table的分析，没有比较模型的输出。先谢谢你，科恩
浏览 1提问于2013-04-02得票数 0
1回答
我有一个连续的反应，三个范畴变量，没有随机效应。我用零调整的伽马分布拟合了一个伽马物体。链接到数据，因为它很大。我无法将数据放到一个更标准的发行版中，如果可能的话，我希望避免使用障碍类型的分析方法，以便更容易地解释。本文描述了如何将gamlss对象转换为nlme::lme对象(，第17页)，这似乎是一种很好的方法，因为统计数据包更广泛地接受lme模型。一旦混合GAMLSS被拟合并存储在一个R对象中，
浏览 4提问于2021-12-13得票数 -1
我在R中运行以下代码行：summary4.157933 0.0000
Number of Groups: 9
但是，我找不到随机变量：deployID的p
浏览 3提问于2021-06-24得票数 0
回答已采纳
2回答
我试图通过RPy分析用lme4.lmer()生成的一组线性模型的偏差。展示了一个完整的示例，我导入了我的dep，加载了我的文件，运行了lme4.lmer()，并且没有让anova在它们上运行。为了您的方便，这里再次是一个粘贴的线，正在失败，我想看到工作。of length 3 not meaningful In an
浏览 1提问于2013-11-12得票数 1
回答已采纳
我一直在使用来自anova.merMod包的lme4函数，通过科学出版物的似然比测试获得固定效果的p值(大多数审查员仍然在我的领域中要求p值)。我注意到，anova.merMod函数使用REML = FALSE重新计算lmer函数(参见下面的示例)，这是一个非常好的特性，迫使不熟悉的用户正确地进行测试。但是，我一直试图阅读lme
浏览 4提问于2014-02-03得票数 2
回答已采纳
1回答
我很难找到一个特别的测试来破译我在什么“会话”或在会话影响内进行治疗的时间。这是我的数据：1 1 1 3.285714286 123 2 1.714285714 101 3 4 0.714285714 12model.b =
浏览 10提问于2017-07-09得票数 0
具体来说，我将它与LME4结合使用。每当我试图使用bfFixefLMER_F.fnc (在LMER模型中固定效果的重新拟合)或pamer.fnc (计算方差或偏差的上、下界p-值)来分析通过LME4拟合的LMER模型时，无论数据集是什么，我都会遇到“在pf(anova.tableterm，"F值”，anova.tableterm，"Df"，nrow(模型@frame)-
浏览 1提问于2020-09-21得票数 0
回答已采纳
2回答
我知道如何将广义线性模型(GLMs)和广义线性混合模型(GLMs)与lme4软件包中的glm和glmer进行拟合。作为统计学专业的学生，我有兴趣学习如何拟合GLM和GLMM<code>E 29</code>以下逐级公式基R码。如果你能指出这方面的任何资源和/或参考，我将不胜感激。提前谢谢。EDiT
我想用公式一步一步地做、GLM、和GLMM，就像我们用矩阵方法做LM一样。是否有使用这种方法<
浏览 2提问于2011-08-03得票数 1
1回答
通常，在aov()上使用summary()函数后，您可以得到剩余值。
但是，当我使用重复测量的方差和公式是不同的时候，我怎么能得到残差？
浏览 4提问于2014-10-02得票数 3
回答已采纳
1回答
在pymer4 (或类似的用于线性混合效果回归模型的python模块)中，有没有一种方法来比较两种不同的模型？在R中，lme4包的命令anova执行以下操作： model1 = lmer(DV ~ IV1 + (1|RV), data = data)
anova(model1, model2) pymer4的API页面似乎没有提到这一点。(我可以用R来做这件事，但我更愿
浏览 23提问于2021-02-04得票数 0
1回答
背景：mod = lme(fixed = y ~ w1 + w2 + x1,我的问题是，我希望能够使用an
浏览 2提问于2016-08-22得票数 1
回答已采纳
2回答
这里的新手。我在lmerTest和stargazer的工作中遇到了问题。我在这里学习本教程，让stargazer在R中使用lme4。
anovaM2 <- lmer(angle ~ factor(temperature) + (1 | replicate) + (
浏览 4提问于2015-07-09得票数 9
回答已采纳
1回答
实验的设计涉及10名参与者。所有患者均经过A、B、C、D条件治疗，而参与者1-5的治疗条件为E、F，参与者6-10的治疗条件为G、H。我正在使用带有lme功能的nlme包来处理丢失的数据，并防止按列表删除参与者。测量变量= DV，固定效应=条件，随机效应=参与者)。当一切都刚刚过去的时候，这就是我所拥有的：
lme(DV~cond, random =~1|ppt, data = OutcomeData, method = "
浏览 1提问于2014-07-27得票数 0
回答已采纳
我正在用lm()修复一个线性回归，比如其中a是具有多个因子水平的阻塞变量。0.4228, Adjusted R-squared: 0.3937 我试图显示与使用ggplot绘制原始数据和回归线时使用摘要( R2 )得到的方程和模型相同的方程和模型，但因为我实际上没有提供一个，所以没有考虑stat_poly_eq()的<
浏览 1提问于2020-10-31得票数 0
1回答
在一个使用lme进行混合效果建模的项目中，我试图比较具有不同相关结构和相同固定部分的模型。由于我将构建许多这些模型(针对不同的因变量)，我尝试编写一个函数来生成具有不同相关结构的模型列表，如下面的示例所示(我确实试图将其保持在最小的工作示例中)。如果我在这个列表中的元素上运行anova()，这是可行的，但是如果在我的全局环境中，则只运行。为什么是这种情况？有什么办法可以避免这个问题
浏览 2提问于2015-09-30得票数 1
回答已采纳
2回答
我有后续的重复测量设计:前测-干预-后测。我怀疑的是，是否必须在随机部分中指定pre_post +期望值+ pre_post:expectations +(1区参与者)+ (1|pre_post)，data
浏览 3提问于2020-05-25得票数 1
当对以下数据(链接：)在R中执行双向rm方差分析时，我遇到了问题：
ANOVA_RTS <- ezANOVA( , type = 2 , return_aov = FALSEANOVA_RTS
结果是正确的(我重复检查使用统计)。然而，当我使用lme函数执行r
浏览 0提问于2018-06-13得票数 1
回答已采纳
1回答
比较速率-曲线和退化
、、、、
我有9个退化曲线，我想比较，并希望就如何最好地这样做的建议。我最初的想法是比较非线性回归。我会先解释这个问题，然后再详细说明实验设计：
我在3个地块(A，B，C)中放置了3种类型的有机质(X，Y，Z)。我知道每个样品的原始有机物(包括总量和干物质百分比)的含量。在每周3个时间点(T1、T2、T3)
浏览 6提问于2015-12-18得票数 0
1回答
在创建这样的mppm模型时：fit1 <- mppm(Points ~ group, simba, interaction=Str, 使用anova.mppm进行似然比测试表明，交互作为一个整体是非常重要<
浏览 6提问于2022-01-02得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
