emmeans中所有预测的平均值

文章/答案/技术大牛

发布

1回答

mean、lm、lme4、emmeans

我在使用emmeans评估所有预测的平均值(或加权平均值)时遇到了问题。例如，混合模型： library(emmeans)m1 <- lmer(mpg ~ 1 + wt + (1|cyl),data=mtcars) 固定效果wt成功： emmeansspecs="wt") 3.22 20.2 1.71 1.83

浏览 31提问于2020-08-07得票数 0

回答已采纳

1回答

for .变量水平的估计

r、regression、emmeans

在这种模式下：mod <- lm(conc~source + percent, data=pigs) percentemmean SE df lower.CL upper.CL我正在使用emmeans包，它给出了percent平均值的估计值。但是，如果我想预测percent值8、10和1

浏览 13提问于2022-10-31得票数 1

回答已采纳

1回答

平均累加变量的意外输出

lme4、emmeans

我进一步使用以下命令：尽管如此，对于变量date_measurement，The均值是平均的。如下面的例子所示，30均值是土壤MT和测量日期27.4中x、y和z基因型的平均值，但测量日期实际上发生在21、23、28、30和35 das上。，运行for均值时，测量数据仍为平均值。date_meas

浏览 6提问于2022-07-11得票数 2

1回答

R中in函数的误差。怎么修呢？

r、emmeans

我有这样的档案：我用这个数据集来预测一个线性混合模型，我想用函数emmeans来计算我的条件的估计平均值。我使用的代码如下：head(newtab5) means <- emmeans(model, "flabel&qu

浏览 1提问于2019-04-10得票数 0

回答已采纳

1回答

我想用Type_f和Type_space解释实验中的Exhaustion_product率和定量变量Age。gl=glm(Type_f ~ Type_space + Exhaustion_product + Age , family=binomial,data=res) emmp <- emmeans( gl, pairwise ~ Exhaustion_product + Age) summary( emmp, infer=TR

浏览 31提问于2018-09-18得票数 2

回答已采纳

1回答

在虚拟变量中使用emmeans

r、emmeans

我有这样的数据： score = c(0, 0, -3, -5, 0, -4,formula = score ~ visit1 + visit2 + visit4 + visit2:trt + visit4:trt + us(ID | visit), data = dat) 我感兴趣的是，在最近一次访问(访问4/V4)的治疗组的LS方法和这方面的对比。使用下面的模型fit2 <-

浏览 13提问于2022-11-10得票数 0

1回答

为什么绘制的预测值和边际效应不同？

r、ggplot2、predict

我正在使用一个具有以下规范的LMER模型： previous_vote_share + coal + coalshare*9033.925

浏览 0提问于2020-05-11得票数 0

2回答

Emmeans函数-参考网格中没有变量

r、lme4、emmeans

我试图在一个较大的数据集上运行emmeans函数，但它不起作用。以下是我的数据： structure(list(Date = structure(c(16578, 16578, 16578, 16578, 16578, 16578), class = "Date"<- emmeans(model1, "Turtle") 以下是我得到的错误： To enable adjustments, add the argument 'pbkrtest.limi

浏览 74提问于2021-10-13得票数 0

2回答

负二项式模型偏移似乎产生了一个2级因子。

r、offset、glm、emmeans、mass

我试图将一些数据拟合成负二项式模型，并使用emmeans进行两两比较。数据有两个不同的样本大小，15和20 (在下面的示例中是num_sample)。我建立了两个数据框架：good.data使用15到20之间的随机样本大小产生offset()的预期结果，而bad.data使用的样本大小为15或20，这似乎产生了一个因子为15或20。与good.data相比，bad.data成对比较产生了太多的比较，尽管它们应该产生相同的数字？->mod.bad

浏览 17提问于2022-08-01得票数 0

回答已采纳

1回答

交流(统计建议)

lme4、interaction

我需要帮助理解并跟踪使用glmer()从lme4获得的交互。这些是我得到的结果在第一次互动之后，我使用电子邮件的方式进行了跟踪： emmeans(model1, list(pairw

浏览 0提问于2019-03-25得票数 0

回答已采纳

2回答

如何在Julia中从连续变量中创建范围？

julia

我正在使用效果包运行Julia的平均边际效应。我的目标是看看不同年龄的男人和女人的体重如何变化。正如你在下面的输出中所看到的那样，它对男性和女性来说，每个年龄的平均边际效应。不过，我想取一个年龄变量的范围，而不是每年单独考虑。例如，我希望有0:5，5:10，10:15等年龄范围。这必须在运行回归模型之后进行，而不是事先进行。我试着自己动手，但我对朱莉娅不够流利。因此，唯一需要纠正的问题是： d1 = Dict(:sex => ["mal

浏览 9提问于2022-05-11得票数 3

回答已采纳

2回答

glmmTMB，后特别测试和glht

r、lme4、mixed-models、posthoc、emmeans

还有一个偏移量，由每个样本中总计数的自然对数(LOG_TOT)组成。固定效应F1的值为0、1、2、3、4和6，固定效应F2的值为0、8、25和75。我把固定效应看作是连续的，而不是顺序的，因为我想识别因变量CT中单调的单向变化，而不是上下变化。emmeans。如果某个善良的人能告诉我如何正确地使用glht和linfct中的公式到glmmTMB的emmeans场景<

浏览 1提问于2018-07-13得票数 2

回答已采纳

1回答

具有R的一般混合线性模型截取的假设检验

r、testing、regression

文件位于：中library(nlme)library(car)library(emmeans)libraryvalue adjustment: bonferroni method for 6 tests 问题如何检验每种基因型的截取年龄是否等于零如何检验每种基因型不

浏览 0提问于2018-10-18得票数 0

1回答

用公共变量进行多对交互时的差异和手工比较

r、regression、interaction、emmeans

但是，我无法获得与emmeans相同的结果，也无法在前面的文章中找到解决方案。情况如下:我有三个变量。var1和var2分为两个层次(A和B，高和低)。另一个变量已经进行了缩放(平均值为0，标准差为1)，因此估计值表示该变量平均时的效果。var1与var2和var3_z相互作用，然后我将估计值与as的输出进行比较，特别是第一个，因为解释是直接的。正如您可以看到的(下面的代码)，当var2为“低”时，A和B之间的区别是1.36，但to意味着它是1.3

浏览 13提问于2022-11-20得票数 0

1回答

映射是从R中的线性模型列表中得到的

r、dictionary、linear-regression、tidyverse、emmeans

我有一个超过100个线性模型的列表，我想为每个模型取估计的平均值和标准误差。 ) 使用emmeans，我可以很容易地得到每个模型的估计均值和标准误差。mod <- lm(hp ~ cyl, data = df)但如果我有一份模特名单呢？is needed 理想情况下，我想要的是能为

浏览 0提问于2018-09-18得票数 3

回答已采纳

1回答

如何从对数(x+1)变量中求绝对差估计和置信区间

r、transform、confidence-interval、emmeans

我有一个混合效应模型，带有一个对数(x+1)转换的响应变量。emmeans的输出类型为"response“，提供了我正在比较的两组的平均值和置信区间。然而，我想要的是组之间差异的平均值和CI (即估计值)。emmeans只提供比率(使用type=“响应”)或对数比率(使用type=“链接”)，我不确定如何将其更改为绝对值。如果您在不使用对数(x+1)转换的情况下运行模型，则emme

浏览 10提问于2019-09-08得票数 0

1回答

当我试图用get获得边际(调整)的意思时，错误消息

r、linear-regression、emmeans

嗨，我在R中做线性回归，并希望为我的子组创建边际(调整)方法。Error: The rows of your requested reference grid would be 19800, which

浏览 15提问于2022-03-17得票数 1

1回答

带not的Box-Cox反向变换不能正常工作

r、linear-regression、lsmeans、emmeans

我正在对需要转换的数据进行线性回归，为此，我使用Box-Cox幂转换，然后使用原始规模的反向转换来编写报告。我一直在尝试使用emmeans包这样做，我遵循了emmeans package vignette中描述的步骤，但是，我发现估计平均值的汇总结果与未转换的数据一点也不相似。事实上，输出根本不会被转换。下面是一个使用emmeans包中的示例的可重现的</

浏览 15提问于2019-03-09得票数 0

3回答

界定不同治疗水平之间的对比

r、contrast、emmeans

我的目标是评估是否有任何组合的“肥料x剂量”是有效的克服“控制”(控制不包括肥料)。为此，我运行以下代码：CLD = cld(marginal, alpha=0.05,reversed=TRUE, Letters=LETTERS)这里的问题是，与100公斤公顷-1的“控制”处理相比，这实际

浏览 9提问于2022-07-07得票数 1

1回答

Weights=“单元”和weights=“比例”在r包中的区别是什么？

r、emmeans

我想用离散预测器和不平衡数据来计算glm模型中的边际均值。利用函数emmeans的emmeans包得到边缘均值，给出了设置weights="cell"和weights="proportional"的不同结果。包装文档称，"proportional"使用的权重与平均的因子组合的频率(在原始数据中)成比例，而"cells"则

浏览 5提问于2021-03-22得票数 1

回答已采纳

点击加载更多