lmer()的后期测试出错: multcomp()和emmeans()
lmer()是一个用于拟合线性混合效应模型的函数，常用于统计学中的多层次建模和混合效应分析。它是R语言中lme4包的一部分。
在进行lmer()函数的后期测试时，可能会出现multcomp()和emmeans()的错误。multcomp()是一个用于多重比较校正的函数，用于对线性混合效应模型进行后续的统计推断和比较。emmeans()是一个用于估计边际均值的函数，用于计算不同组别或条件下的均值。
出现这个错误可能是由于以下原因之一：
- 包版本不兼容：multcomp()和emmeans()函数可能需要特定版本的R语言或相关包才能正常工作。建议检查并更新相关包的版本。
- 数据不符合要求：multcomp()和emmeans()函数对输入数据有一定的要求，例如数据的格式、缺失值处理等。请确保数据符合函数的要求。
- 函数调用错误：可能是由于函数的调用方式不正确导致的错误。请仔细检查函数的参数和语法是否正确。
针对这个错误，可以尝试以下解决方法：
- 更新包版本：使用最新版本的R语言和相关包，确保包的兼容性。
- 检查数据格式：确保输入数据符合multcomp()和emmeans()函数的要求，例如数据类型、缺失值处理等。
- 检查函数调用：仔细检查multcomp()和emmeans()函数的参数和语法是否正确，确保函数的正确调用。
- 查阅文档和示例：查阅相关函数的文档和示例，了解函数的使用方法和常见问题的解决方案。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云计算服务：https://cloud.tencent.com/product
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云服务器：https://cloud.tencent.com/product/cvm
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云物联网：https://cloud.tencent.com/product/iot
- 腾讯云移动开发：https://cloud.tencent.com/product/mobdev
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云元宇宙：https://cloud.tencent.com/product/vr
相关·内容
我正在运行这个混合模型：
dflm1 <- lmer(formula=df[,16] ~ trt + (1|id), data=df)
emm <- emmeans(dflm1,'trt')
library(multcomp)
cld(emm)
Error in UseMethod("cld") :
no applicable method for 'cld' applied to an object of class "data.frame"
不幸的是，我不能上传数据。我已经看过关于这个的其他问题，但它不
浏览 3提问于2020-02-20得票数 0
2回答
不工作在lmer中的动态协变量
、、、、
我试图在lmer模型中使用动态协变量(取自用户输入)。这很好，但是应用emmeans不再起作用了。
示例：
fit1 = lmer(get(uservar) ~ Intervention + period + sequence + age + (1|ID),
data = mydata, REML=T)
和
emmeans(fit1, list(pairwise ~ Intervention), adjust = "tukey")
毫无问题地奔跑。
lmer(get(uservar) ~ Intervention + period + sequence
浏览 1提问于2019-05-10得票数 0
lhiry1 <- lmer(ACC ~ LR1st +(1|Subject),data = learndata_long3) # hierarchy 1*
lhiry2 <- lmer(ACC ~ LR2nd +(1|Subject),data = learndata_long3) # hierarchy 2*
onelmer=c(lhiry1,lhiry2)
onelmer1=c('LR1st','LR2nd')
for (i in 1:length(onelmer))
{onelmercontr[i]=lapply(onelmer,
fu
浏览 2提问于2020-03-05得票数 1
回答已采纳
由于缺乏更好的措辞:我需要从现有的数据帧创建一个新的数据帧，以运行lm/lmer模型来获得模型预测。
假设我有：
x <- as.numeric(rep(1:6,5))
y <- as.numeric(rep(1:5,6))
int1 <- factor(rep(c("a","b"),15))
int2 <- factor(rep(c("11","12","13"),10))
g <- rep(c("f","m")
浏览 6提问于2019-07-31得票数 0
回答已采纳
我对此真的很陌生。如果能帮上忙我会很感激！我导入的dataframe是一个包含所有结果的大表，我真的希望使用相同的lmer函数遍历每一列。 for (i in d[c(6:45)]){
para = colnames(d[c(6:45)])[i] #para is parameter
d.lmer <- lmer(get(para) ~ Group*Sex*CS.NCS + (1|Dam), data = d)
d.anova <- anova(d.lmer)
library(emmeans)
d.compare <- emmeans(d.
浏览 12提问于2020-08-25得票数 0
我有三个因素:单词、类型和语域。在SPSS中，很容易在SPSS中进行成对比较(或简单比较)，语法是：
/EMMEANS=TABLES(word*register*type) COMPARE(type) ADJ (BONFERRONI)
它会给我一个这样的结果：
但是我如何在R中使用Multcomp包实现这一点呢？假设我有一个这样的anova结果：
library (multcomp)
anova <- aov (data, min~word*register*type)
summary (anova)
如何进行类似于SPSS的成对比较？
浏览 3提问于2016-06-22得票数 0
我使用混合模型来计算特定位置(4类)在R中选定频率(连续)的速度估计。
到目前为止，这就是我所拥有的：
lmer.df = lmer(log(velocity) ~ location + frequency + (1 | heart), data = df)
emm = emmeans(lmer.df, specs = ~ location, type = "response")
summary(emm, infer = TRUE) #back-transformed estimated marginal means (for mean frequency)
summar
浏览 10提问于2022-07-06得票数 0
回答已采纳
1回答
我正在尝试使用下面的数据执行混合模型拟合。
df.urbana <- structure(list(Location = structure(c(2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Monmouth",
"Urbana"), class = "factor"), treatment = structure(c(1L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L,
浏览 10提问于2019-11-09得票数 0
我有一个在不同位置测量"Y“的数据集，我试图通过运行lmer()模型并分析结果来确定变量Y如何受到变量A、B和D的影响。然而，当我到达post步骤时，我在尝试分析时收到一个错误。 以下是我的数据示例： table <- " ID location A B C D Y
1 1 AA 0 0.6181587 -29.67 14.14 168.041
2 2 AA 1 0.5816176 -29.42 14.21 200.991
3 3 AA 2 0.4289670 -28.57
浏览 50提问于2021-09-28得票数 0
根据，我得出了以下结论：
library(papaja)
library(datasets)
library(stats)
library(multcomp)
bla <- glm('MPG.highway ~ DriveTrain*Origin',data = Cars93)
contrast.matrix <- rbind('main' = c(0,1,0,0,0.5,0), 'int' = c(0,0,0,0,1,0))
blac <- glht(bla, linfct = contrast.matrix)
apa_
浏览 3提问于2020-07-14得票数 2
1回答
从R中的摘要对象中提取df
、、、、
在下面的R代码中，我想知道如何从df对象中提取ems$emmeans值？
library(lme4)
library(emmeans)
dat <- read.csv('https://raw.githubusercontent.com/hkil/m/master/z.csv')
dat$year <- as.factor(dat$year)
m1 <- lmer(y~year*group + (1|stid), data = dat)
ems <- emmeans(m1, pairwise ~ group*year, infer = c(T, T)
浏览 2提问于2020-05-19得票数 1
回答已采纳
1回答
我试图解决这个错误：
在if ((misc$estType ==“== ")和(粘贴(c(”"，by)，折叠=“，”))中出现错误!=：缺失值，其中真/假需要
我的代码：
model1 <- lmer (Y ~ X1 * X2 + (1|"FACTOR of 13 levels"), data=data, REML = FALSE)
library (emmeans)
post_hoc <- emmeans (model1, ~ X1*X2)
pairs (post_hoc, adjust="tukey")
但它返回上述错误。我
浏览 4提问于2021-11-24得票数 0
回答已采纳
1回答
为了满足lmer中线性模型的假设，我必须转换变量响应(例如变量1)，使用这里建议的重尾数据的方法，并演示如下：
TransformVariable1 <- sqrt(abs(Variable1 - median(Variable1))
然后，我将数据与以下示例模型相匹配：
fit <- lmer(TransformVariable1 ~ x + y + (1|z), data = dataframe)
接下来，我更新引用网格，以便按照这里的建议对转换进行解释( )
rg <- update(ref_grid(fit), tran = "TransformVaria
浏览 12提问于2022-06-12得票数 1
回答已采纳
1回答
我有一个连续响应变量yld和一个分类预测器check (具有3个级别)。我做了一次单向的方差分析和一次临时测试，看哪种水平是不同的。
mdl<-aov(sqrt(var$yld) ~ var$check); summary(mdl);TukeyHSD(mdl)
Df Sum Sq Mean Sq F value Pr(>F)
var$check 2 5162 2581.2 13.51 1.46e-06 ***
Residuals 2775 530395 191.1
浏览 5提问于2016-11-02得票数 2
回答已采纳
我正在分析一个双向重复测量数据集，建模如下：
model= lmer(result ~ treatment * time + (1|subject), data=df)
..。在那里每个受试者都会接受每一次治疗并在每次都进行测试。然而，当分析对比时，似乎没有对球度进行校正。在这里，我使用emmeans来测试每个处理和对照处理之间的差异，在每个“时间”的水平上。
emm <- emmeans(model,c("treatment","time")
contrast(emm, "trt.vs.ctrl", ref="Control
浏览 6提问于2021-03-07得票数 0
1回答
我在使用emmeans评估所有预测的平均值(或加权平均值)时遇到了问题。例如，混合模型： library(emmeans)
library(lme4)
m1 <- lmer(mpg ~ 1 + wt + (1|cyl),data=mtcars) 固定效果wt成功： emmeans(m1,specs="wt")
wt emmean SE df lower.CL upper.CL
3.22 20.2 1.71 1.83 12.1 28.3 然而，为了计算预测的平均值，以下方法以前有效(大约12个月前)，但现在失败了： emmeans(m1
浏览 31提问于2020-08-07得票数 0
回答已采纳
我使用emmeans包，以便从我使用lme4包生成的线性混合模型中生成估计均值。下面是我用来生成线性混合模型的代码(我选择排除dput()，因为它太大，太麻烦了)：
countsleepdataoff25WAKEANGLMM <- lmer(unlist(propnssANG)~
as.factor(intoffidhrtime)+
Daysincedisturb+
as.factor(intoffidhrtime)*Daysincedisturb+
(1|Batch)+(1|Batch_Bird_ID)+(1|fdfsrtoff)+
(1|fdfsrtof
浏览 21提问于2022-11-16得票数 0
回答已采纳
这只是一个一般的问题，获得信任区间的交互手段，我已经阅读了所有的普通教程，但我不知道如何为双向和3路互动。这是一个3路交互的例子.
X= continous variable
A= 3 levels, nested in L
L= 2 levels
G= 2 levels, crossed with L (and thus also A)
ID= participants nested in G
然后我用随机截距拟合了一个线性混合模型：
fit1 <- lmer(X~G*A*L+(1|ID),data=df)
emms=emmeans(fit1,specs=pairwi
浏览 0提问于2020-01-14得票数 0
回答已采纳
1回答
我有这样的档案：
我用这个数据集来预测一个线性混合模型，我想用函数emmeans来计算我的条件的估计平均值。我使用的代码如下：
newtab5 <- read.csv(file="sliding_LMM_test.csv", header=T, sep=",")
head(newtab5)
library(lmerTest)
model <- lmer(formula = data1 ~ flabel + (1 | whichFragments), data = newtab5, control=lmerControl(check.nl
浏览 1提问于2019-04-10得票数 0
回答已采纳
1回答
假设我有这些数据 library(MASS)
m<-lmer(Y~N*V + (1|B),data=oats) 如何在emmeans中创建手动对比度？例如 Victoria_0.2cwt 1
Victoria_0.4cwt -1
Marvellous_0.2cwt -1
Marvellous_0.4cwt 1
浏览 21提问于2019-02-28得票数 0
回答已采纳
2回答
我一直在尝试计算我在R中的lmer & glmer的边际均值。我找到了emmeans函数，我一直在努力理解它，并将其应用到我的模型中。我发现很难获得交互的方法，所以我从加性预测开始，但函数的工作方式与示例中显示的方式不同(例如，这里的)
emmeans(Oats.lmer, "nitro")
nitro emmean SE df lower.CL upper.CL
0.0 78.89207 7.294379 7.78 61.98930 95.79484
0.2 97.03425 7.136271 7.19 80.25029 113
浏览 14提问于2018-11-23得票数 1
1回答
双变量线性混合模型的后自组织检验
、、、、
我建立了一个线性混合模型，并对其进行了后特别测试。固定的因素是相数(时间)和组。
statistic_of_comp <- function (x, df) {
x.full.1 <- lmer(x ~ phase_num + group + (1|mouse), data=df, REML = FALSE)
x_phase.null.1 <- lmer(x ~ group + (1|mouse), data=df, REML = FALSE)
print(anova (x.full.1, x_phase.null.1))
summary(glht(x.full.1, li
浏览 0提问于2019-02-15得票数 1
在这个例子中，我希望在不同的治疗方法之间进行多个比较。以下是数据：
data.1 <-read.csv(text = "
location,treat,response
loc1,T1,120
loc1,T2,60
loc1,T3,59
loc1,T4,10
loc2,T1,129
loc2,T2,55
loc2,T3,59
loc2,T4,8
loc3,T1,134
loc3,T2,60
loc3,T3,58
")
我就是这样做的：
library(lme4)
library(lmerTest)
library(emmeans)
library(multcomp)
m
浏览 17提问于2022-01-26得票数 1
回答已采纳
我有以下情况：
我的固定效应模型发现了Relation_PenultimateLast在被称为“作曲家”的参与者群体中的主要效应。因此，我想找出Relation_PenultimateLast的水平在统计上与其他级别不同。
f.e.model.composers = lmer(Score ~ Relation_PenultimateLast + (1|TrajectoryType) + (1|StimulusType) + (1|Relation_FirstLast) + (1|LastPosition), data=datasheet.complete.composers)
摘要(
浏览 2提问于2018-01-13得票数 1
回答已采纳
我有一个假的数据集，描述空气温度对两种植物(a和b)生长的影响。
data1 <- read.csv(text = "
year,block,specie,temperature,growth
2019,1,a,0,7.217496163
2019,1,a,1,2.809792001
2019,1,a,2,16.09505635
2019,1,a,3,24.52673264
2019,1,a,4,49.98455022
2019,1,a,5,35.78568291
2019,2,a,0,8.332533323
2019,2,a,1,16.5997836
2019,2,a,2,11
浏览 8提问于2022-08-23得票数 1
回答已采纳
1回答
我的全部代码运行得非常好。当我试图使用Knitr生成PDF/ HTML时，我得到了一个错误：
Error in Anova(mod1, type = "III") : could not find function "Anova"
对我能做什么有什么不同的建议吗？
Relevant code:
library(lsmeans)
library(lme4)
df <- read.csv(file = "file.csv", header = TRUE)
library(lmerTest)
mod1 <- lm(Yield ~ Trea
浏览 4提问于2022-02-14得票数 1
回答已采纳
2回答
我试图在一个较大的数据集上运行emmeans函数，但它不起作用。以下是我的数据： structure(list(Date = structure(c(16578, 16578, 16578, 16578,
16578, 16578), class = "Date"), Time = c(7, 7, 7, 9, 11, 11),
Turtle = c("R3L12", "R3L12", "R3L12", "R3L12", "R3L12", "R3L12"
), T
浏览 74提问于2021-10-13得票数 0
1回答
我想测试一个2x3的析因设计，并像这样对比变量 library(lme4)
library(emmeans)
my.helmert = matrix(c(2, -1, -1, 0, -1, 1), ncol = 2)
contrasts(Target3$mask) = my.helmert
contrasts(Target3$length) 因此，对于掩码，我想将第一组与其他两组的平均值进行比较，并在第二步中将第二组与第三组进行比较。 这在我的LMM中运行良好 Target3.2_TT.lmer = lmer(logTotalTime ~ mask*
浏览 36提问于2019-01-20得票数 0
回答已采纳
1回答
你好，我正在尝试用R来可视化一些数据。然而，数据在横坐标中显示的顺序是错误的。我试图修复这个函数scale_x_discrete(limits=c("Fallow", "Mustard", "Mix4", "Mix12"))，这部分地解决了这个问题，因为一些小写字母(a、ac、bc、c)的排序不正确。正确的顺序是(a，ab，bc，c)。
这是情节的代码：
library(tidyverse)
library(lme4)
library(emmeans)
library(multcomp)
# set theme for ggpl
浏览 2提问于2020-11-21得票数 0
回答已采纳
嗨，我有一个lmer模型，在这个模型中，我试图看看mouse_hit_target (命中/误)和click_number (第一个/秒)对QE_duration的预测有多大。
duration <- lmer(QE_duration ~ mouse_hit_target + click_number + (1|pid), data=lmm.data_2)
emmeans(duration, list(pairwise ~ click_number : mouse_hit_target), adjust = "tukey" )
当尝试比较不同条件下的方法时，对于“第一次命
浏览 3提问于2022-05-11得票数 1
1回答
我在线性混合效果模型的posthoc比较中遇到了一个问题。我将试着用一个快速构建的不完美的例子来解释它：
下面是我的示例数据：
Variable<-as.factor(rep(c(1,2,3),5))
Random<-rep(c(1,2,2),5)
Result<-rnorm(15,mean=10,sd=2)
Data<-as.data.frame(cbind(Variable,Random,Result))
我的模型中实际上包含了几个固定的和随机的效果，但这足以说明我的问题：
library(lme4)
LME=lmer(Result~Variable+(1|Ran
浏览 24提问于2021-02-19得票数 0
回答已采纳
1回答
我使用简单的平方变换转换了一个变量(例如leaf_area)，然后对包含交互作用的以下模型进行了拟合：
fit <- lmer(leaf_area^2 ~genotype*soil_type + date_measurement + light + (1|repetition) + (1|y_position) + (1|x_position), data = dataset)
为了获得每个测量日期的平均累积基因型和土壤类型，我进一步使用以下命令：
fit.emm <- emmeans(fit, ~ genotype*soil_type + date_measurement, ty
浏览 6提问于2022-07-11得票数 2
1回答
使用Tukey's测试校正效果
、、、、
library(lme4)
library(emmeans)
m <- lmer(angle ~ recipe*temp + (1|replicate), data=cake)
我感兴趣的是每个食谱的简单效果(不是成对比较)，所以基于这个
out.emtrends <- emtrends(m, pairwise~recipe, var="temp")
emtr <- as.data.frame(out.emtrends$emtrends)
tvalues <- emtr$temp.trend/emtr$SE
dfs <- emtr$df
pvalu
浏览 14提问于2019-08-28得票数 1
我正在运行的对比的p值没有被正确地转换为data.frame。这是为什么，我该怎么解决呢？
电子设备的控制台输出：
> pairs(emmeans(lmer.mod, ~ Status*Stim*Treatment), simple = "each")
$`simple contrasts for Status`
Stim = 1, Treatment = None:
contrast estimate SE df t.ratio p.value
Control - Subclinical -0.24213 0.0571 57
浏览 2提问于2021-03-14得票数 3
回答已采纳
2回答
我使用glmmTMB分析负二项式广义线性混合模型(GLMM)，其中因变量是计数数据(CT)，它是过分散的。
在相关数据帧中有115个样本(行)。有两个固定效应(F1，F2)和一个随机截距(R)，其中嵌套了另一个随机效应(NR)。还有一个偏移量，由每个样本中总计数的自然对数(LOG_TOT)组成。
数据框架df的一个例子是：
CT F1 F2 R NR LOG_TOT
77 0 0 1 1 12.9
167 0 0 2 6 13.7
289 0 0 3 11 13.9
253 0 0 4 16 13.9
125 0 0
浏览 1提问于2018-07-13得票数 2
回答已采纳
1回答
我已经看到了几个例子，如何可能选择理想的成对比较，但不幸的是，不知道如何应用于我的数据。
以下是我的简化数据集：
# FIXED effects: age and brain_area
df$age <- factor(df$age)
df$brain_area <- factor(df$brain_area)
# RANDOM effects: subject_ID and section
df$subject_ID <- factor(df$subject_ID)
df$section <- factor(df$section)
# dependent va
浏览 4提问于2022-01-21得票数 1
回答已采纳
1回答
我想知道为什么map不能在我下面的可重现的例子上运行？ library(purrr)
library(emmeans)
library(lme4)
dat <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/dat2.csv')
fit <- lmer(math ~ time*ses + I(time^2)*ses+ (time|id), data = dat)
map(0:7,
~emtrends(fit, ~ses, var= "time"
浏览 5提问于2020-11-17得票数 0
回答已采纳
1回答
我有一个纵向研究的数据，并计算回归使用lmer 4：：lmer函数。之后，我计算了这些数据的对比，但我很难解释我的结果，因为它们是出乎意料的。我想我可能是在密码上弄错了。不幸的是，我不能用一个例子来复制我的结果，但是我会在下面发布失败的例子和我的实际结果。
我的研究结果：
library(lme4)
library(lmerTest)
library(emmeans)
#regression
regmemory <- lmer(memory ~ as.factor(QuartileConsumption)*Age+
(1 + Age | ID) +
浏览 4提问于2022-04-27得票数 0
回答已采纳
1回答
我建立了一个线性混合模型来分析我的数据，并使用软件包multcomp进行配对比较。我已经用模型输出创建了一个ggplot，我希望将输出(字母)放在栏之上。
但是我不知道如何从cld()输出中提取字母。这是一个包含许多对象和向量的清单，坦白地说，我不完全理解其中的一半(而且，虽然我读过几次，做了很多基本的练习，但我仍然很难理解清单)。
在R中运行cld()时，它会准确地显示控制台中所需的内容，这就是我想要提取的内容。到目前为止，我的互联网搜索还没有找到任何解决方案。
代码
# packages
library(lme4)
library(multcomp)
library(ggplot2)
#
浏览 0提问于2018-12-19得票数 0
回答已采纳
2回答
R函数的最后一行阻止绘图
、、、、
我真的想知道为什么当我将最后一行(sigma)保留在函数foo1中，我的plot调用就停止工作了，但是当我删除foo2中的最后一行时，plot调用工作得很好吗？！
注意:我的理解是，情节应该显示在函数中的任何位置，例如：
foo3 <- function(x = 1:3){; plot(x); return(x); }; foo3()
我需要保持最后一行，但我也需要情节，这是固定的吗？
library(lme4)
library(emmeans)
h <- read.csv('https://raw.githubusercontent.com/hkil/m/master/
浏览 0提问于2020-05-21得票数 3
回答已采纳
我正在用图书馆重复测量R中的变异数：
library(ordinal)
library(car)
library(RVAideMemoire)
我有两组:月份和距离，因变量是CO2：
距离月CO2
0米7月234
我为CO2建立了一个clmm模型，解释了距离、月份以及月份和距离之间的交互作用：
model_CO2 = clmm(CO2.f ~ month + distance + month:distance + (1|nest),
data = field_data,
threshold = "equidistant")
浏览 2提问于2020-04-24得票数 1
回答已采纳
1回答
我试图用emmeans包中的emmeans函数在R中进行后期测试。然而，我无法找出我应该在规范的论点中写些什么。据我所知，这是我想要对比的变量的位置(我的自变量)。然而，当我把我的静脉注射/试剂盒，它会产生错误。我把我的代码和错误放在下面：
X= as.factor(rep(c("A", "B", "C"), each= 50))
Y= as.factor(c("K", "L", "M"), times= 50)
Z= rnorm(150)
Model= lm(Z ~ X+Y+X:Y)
emmean
浏览 2提问于2018-05-12得票数 0
回答已采纳
1回答
这是一个被嵌套的数据集，当我想要看到饮食和周之间的两两比较时，这就成了问题。错误“尝试从‘by’中删除嵌套因子”是什么？卑劣？
form <- as.formula(paste(colnames(df)[8],'~ Diet + period +week*Diet +(1|id)')) #get data for interactions
dflmer <- lmer(form, data=df)
a <- Anova(dflmer, type=3)
library(emmeans)
emm <- emmeans(df
浏览 6提问于2020-04-23得票数 0
回答已采纳
我正在尝试将我的excel数据集连接到R以进行统计分析。我在OSX Sierra (10.12.6)、R studio (1.0.153)和Java8(更新144)上运行。
一周前，函数"read_excel“能够打开我的excel文档。当我将excel和R文档一起移动到另一个文件夹时，它不再起作用。重新加载库没有任何效果。经过多次尝试(并重新启动R studio和计算机)，终于可以工作了，但再也找不到函数"lmer“。重新加载库"lme4“后，"read_excel”不再起作用！
我也尝试过使用"read.xlsx“和"readWorkshee
浏览 0提问于2017-10-12得票数 5
回答已采纳
3回答
我正在分析5个处理的受精试验，即：
1_Control，2_Fertilizer_A，3_Fertilizer_B，4_Fertilizer_C和5_Fertilizer_D
在不同剂量(40、80、120和160 Kg公顷-1)下，对肥料A~D进行了评价。
我的目标是评估是否有任何组合的“肥料x剂量”是有效的克服“控制”(控制不包括肥料)。
我对以下模型进行了拟合，得到了二次响应曲线，并找到了最优响应率：
ModelFert <- lmer (yield ~ treatment + dose + dose2 +
(1 | environment/block
浏览 9提问于2022-07-07得票数 1
1回答
在pymer4 (或类似的用于线性混合效果回归模型的python模块)中，有没有一种方法来比较两种不同的模型？ 在R中，lme4包的命令anova执行以下操作： model1 = lmer(DV ~ IV1 + (1|RV), data = data)
model2 = lmer(DV ~ IV1*IV2 + (1|RV), data = data)
anova(model1, model2) pymer4的API页面似乎没有提到这一点。(我可以用R来做这件事，但我更愿意用python来做，因为我必须通过一个for循环在不同的组中运行这个测试，而python为这个循环提供了灵活性。)
浏览 23提问于2021-02-04得票数 0
我运行了一个混合模型与lmerTest，我需要一个后特别的测试。
下面是模型，每一次试用都是连续的。然后，我使用number，但得到以下错误(可能是由于大量的观察)。
我会感激帮助或其他方式运行后挂钩，因为其他方式，我也尝试了不起作用。当我试图添加参数'lmerTest.limit = 13812‘时，它根本不起作用。
lmerTest :: lmer (RT ~ condition * pronoun * objectification_center +
(1+ pronoun| subject_ID),
data = data)```
emm1 = emmeans(mixed_
浏览 3提问于2021-05-31得票数 0
回答已采纳
这是我的数据：
subject arm treat bline change
'subject1' 'L' N 6.3597 4.9281
'subject1' 'R' T 10.3499 1.8915
'subject3' 'L' N 12.4108 -0.9008
'subject3' 'R' T 13.2422 -0.7357
'subject4' 'L' T 8.7383 2.7
浏览 5提问于2021-04-23得票数 0
回答已采纳
我用我的数据集(不相等的样本大小，不相等的方差)在r中运行了双向anova :1个变量在3个物种中测量(每个物种中有男性和女性)。这在物种之间产生了显着的结果，所以我想知道哪些成对比较产生了显着的意义。我知道在R中有用于执行后期测试的包中的函数:例如
来自的邓内特的post测试。需要的包："multcomp"，"mvtnorm"，"survival"，"splines“
library(multcomp)
test.dunnett=glht(anova_results,linfct=mcp(method="D
浏览 9提问于2013-09-18得票数 2
我正在尝试将我在其他地方计算的p值添加到我的ggplot barplot中。
d <- read.table("DS4 R Analysis.csv",header=TRUE,sep = ",")
d.lmer <- lmer(Corticosterone ~ Group*Sex*CS.NCS + (1|Dam), data = d)
d.anova <- anova(d.lmer)
d.compare <- emmeans(d.lmer, ~ Sex*Group*CS.NCS)
d.pairs <- pairs(d.compar
浏览 0提问于2020-09-02得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
