使用agricolae::HSD.test使用purrr::map分析lme4::lmer主题
agricolae::HSD.test是一个R语言中的函数，用于进行多重比较分析。它可以帮助我们在统计学中进行方差分析后，对不同组之间的差异进行比较。
该函数的主要参数包括：
model：方差分析模型对象，通常是通过aov()函数创建的。trt：处理组的因子变量。group：分组方式，可以是数字或字符向量。console：是否在控制台中打印结果。alpha：显著性水平，默认为0.05。
使用agricolae::HSD.test函数可以进行多个组之间的两两比较，以确定是否存在显著差异。该函数返回一个包含比较结果的数据框，其中包括组别、均值、标准误差、上下界限、显著性水平等信息。
在云计算领域中，agricolae::HSD.test函数可以应用于数据分析和统计建模的过程中，用于比较不同处理组之间的差异。例如，在云计算平台上进行用户行为分析时，可以使用该函数对不同用户群体的行为数据进行比较，以确定是否存在显著差异。
腾讯云提供了一系列与数据分析和统计建模相关的产品和服务，可以与agricolae::HSD.test函数结合使用。其中，腾讯云的数据仓库产品TencentDB for PostgreSQL可以用于存储和管理大规模数据集，而腾讯云的弹性MapReduce（EMR）服务可以用于在大规模数据集上进行并行计算和分析。此外，腾讯云还提供了一系列人工智能相关的产品和服务，如腾讯云机器学习平台（Tencent ML-Platform）和腾讯云自然语言处理（Tencent NLP），可以用于数据分析和模型构建过程中的特征工程和算法优化。
更多关于腾讯云数据分析和人工智能产品的详细介绍和使用方法，可以参考以下链接：
需要注意的是，以上答案仅供参考，具体的产品选择和使用方法应根据实际需求和情况进行决策。
相关·内容
我希望从包中do.call一个(导出的)函数，并需要将what参数中的包和函数指定为字符串(也就是说，当包未加载到我的环境中或存在潜在的函数名称冲突时)。
函数名是一个字符串，如"lme4::lmer"，它既指定，也指定包和函数的。
例如，有条件地调用一个函数(类似于)：
FUN = if(TRUE) {
"lme4::lmer"
} else {
"nlme::nmle"
}
args.list = list(Reaction ~ Days + (Days | Subject),
quote(lme
浏览 1提问于2016-08-16得票数 5
回答已采纳
我很难理解如何使用map来应用一个函数来分析我数据中各个主题的数据
library(lme4)
library(tidyverse)
library(agricolae)
这是我开始使用的模型
fm <- lmer(Reaction ~ Days + (1 | Subject), data = sleepstudy)
我可以查看个人的数据：
glimpse(ss <- sleepstudy %>%
nest(-Subject) %>%
mutate(fit = map(data, ~ lm(Reaction ~ Days, data = .)),
浏览 16提问于2017-07-13得票数 1
1回答
我正在尝试使用和multilevelmod从tidymodels中提取lme4随机拦截。我可以使用下面的lme4来完成这个任务：
使用R和lme4的lme4:
library("tidyverse")
library("lme4")
# set up model
mod <- lmer(Reaction ~ Days + (1|Subject),data=sleepstudy)
# create expanded df
expanded_df <- with(sleepstudy,
data.frame(
浏览 2提问于2021-08-31得票数 2
回答已采纳
3回答
我想知道如何用lme4拟合多元线性混合模型。我用以下代码拟合了单变量线性混合模型：
library(lme4)
lmer.m1 <- lmer(Y1~A*B+(1|Block)+(1|Block:A), data=Data)
summary(lmer.m1)
anova(lmer.m1)
lmer.m2 <- lmer(Y2~A*B+(1|Block)+(1|Block:A), data=Data)
summary(lmer.m2)
anova(lmer.m2)
我想知道如何用lme4拟合多元线性混合模型。数据如下：
Block A B Y1 Y2
1 1 1
浏览 0提问于2011-10-20得票数 15
回答已采纳
5回答
我想知道是否有一种方法可以编写逻辑测试(TRUE/FALSE)来显示来自lme4包的模型是否已经收敛？
下面显示了一个示例，如果有任何模型附带收敛警告(即Model failed to converge)消息，我想要捕获吗？
library(lme4)
dat <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/nc.csv')
m <- lmer(math ~ ses*sector + (ses | sch.id), data = dat)
Warning message:
I
浏览 5提问于2020-09-29得票数 1
回答已采纳
1回答
lmer的替代优化算法
、、、
默认情况下，lmer包中的lme4函数使用来自minqa包的bobyqa作为优化算法。
根据下面的文章，也可以在minqa包中使用其他优化算法
如何使用uobyqa或newuoa作为lmer的优化算法
library(lme4)
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy, control=lmerControl(optimizer="bobyqa"))
浏览 3提问于2014-08-05得票数 2
回答已采纳
1回答
考虑lme4包中的睡眠学习数据，如下所示。其中18名受试者在不同的时间重复测量了他们的反应。假设有一个附加变量Days2，如下所示：
library("lme4")
sleepstudy$Days2<- rep(10:18, rep(20,9))
sleepstudy
Reaction Days Subject Age Days2
1 249.5600 0 308 20 10
2 258.7047 1 308 20 10
3 250.8006 2 308 20 10
...
178 343.2199
浏览 3提问于2021-11-20得票数 1
回答已采纳
1回答
我刚开始在R中打包mice，但我正在尝试从popmis中计算5个数据集，然后分别拟合一个lmer()模型with()，最后在它们之间进行pool()。 我认为mice()中的pool()函数不能与lme4包中的lmer()调用一起工作，对吧？ 如果是这样的话，有没有一种方法可以编写一个定制的函数，在我下面的例子中就像pool()一样？ library(mice)
library(lme4)
imp <- mice(popmis, m = 5) # `popmis` is a dataset from `mice`
fit <- with(imp, lme4::lmer(pop
浏览 19提问于2020-11-08得票数 3
回答已采纳
我在编写循环方面是新手，在这方面我有一些困难。我已经看过了其他问题，但没有找到我的具体问题的答案。
因此，让我们创建一个随机数据集，给出列名，并将变量设置为字符：
d<-data.frame(replicate(4,sample(1:9,197,rep=TRUE)))
colnames(d)<-c("variable1","variable2","trait1","trait2")
d$variable1<-as.character(d$variable1)
d$variable2<-as.characte
浏览 8提问于2018-12-18得票数 1
回答已采纳
1回答
大家好，
当'lme4‘中的'lmer’函数运行时，它会生成一个带有大量插槽的S4对象。我感兴趣的是其中一个插槽，即模型@X，以及这个“X”槽输出是如何产生的。我想尝试为我使用的另一个模型函数(glmmPQL)再现这个输出，它不会自动产生这个'X‘输出(FYI 'lmer’生成一个类'mer‘的对象，而槽'X’是一个固定效果的模型矩阵)。下面的代码显示了我所说的.
我想知道的是这个“X”数据是如何产生的？我看了一下'lmer‘的代码，把它写在终端中而没有'()’，但是在那里什么都找不到。我也尝试了showMethod('l
浏览 2提问于2013-01-17得票数 1
我正在使用一个脚本，该脚本调用lme4包的lme4函数数千次(请不要担心，稍后将执行多个比较的相关更正)，并希望尽可能节省单个调用期间的时间。
我想从拟合的模型中提取t值，最快的(计算时间)方法是什么？我试过使用summary(model)，但它似乎比调用lmer本身花费的时间更长。在不使用summary()的情况下，可以从所获得的模型中得到t-值吗？
浏览 0提问于2014-05-11得票数 3
回答已采纳
1回答
我使用来自lme4的mice包和lmer进行分析。但是，pool.r.squared()在此输出上不起作用。我正在寻找关于如何在以下工作流程中包括调整后的R平方计算的建议。 require(lme4, mice)
imp <- mice(nhanes)
imp2 <- mice::complete(imp, "all") # This step is necessary in my analyses to include other variables/covariates following the multiple imputation
fit <- l
浏览 37提问于2020-07-03得票数 1
回答已采纳
我的结构如下：
'data.frame': 1041 obs. of 4 variables:
$ RST_AST: num 0.45 0.576 0.45 0.45 0.675 ...
$ subj : chr "4" "13" "24" "25" ...
$ domain : chr "perceptual" "perceptual" "perceptual" "perceptual" ...
$ ce : chr
浏览 2提问于2022-01-27得票数 -1
我正在编写一个自定义函数，它在purrr::map的帮助下为列表中的每个元素建立线性混合效果模型。代码块运行得非常好，但是当我将其转换为自定义函数时，还不清楚应该如何输入与list元素中的单个列对应的参数。
如果我让这个自定义函数正常工作，我可以用它来处理任意多的变量。否则，我将不得不对不同的变量保持复制粘贴相同的代码。
# libraries needed
library(purrr)
library(lmerTest)
data(mtcars)
# create a list of dataframes from mtcars based on a split
group_list &l
浏览 0提问于2018-01-28得票数 3
回答已采纳
尝试运行cld函数来显示显著的差异，但它不起作用。
最后，程序在列表(“cld”)中显示->错误：“cld”没有适用于“c(‘emm_UseMethod’，‘list’)”类的对象的适用方法。
所有必需的软件包都已安装。有人能帮我吗？
library(emmeans)
library(lattice)
library(pbkrtest)
library(car)
library(agricolae)
library(agriTutorial)
library(multcompView)
library(car)
View(WG)
attach(WG)
names(WG)
WG$PGPM
浏览 2提问于2020-12-21得票数 0
我有一个数据集dat2，我想在其上拟合一个线性混合效果模型。我在过去的pvals.fnc补丁中使用了lmer() (lme4包)来计算相关的p值。
但是，由于我使用新的lme4 (1.0.4)和languageR (1.4)包重新安装了R 3.0.2版本，因此我得到了一个关于lmer函数输出的错误。它说明输出不是mer对象。实际上，它的类是lmeRmod。
下面是我使用的代码：
names(dat2)<-c("auc","subj","decod","soa","vis")
attach(dat2)
mod1
浏览 4提问于2013-10-05得票数 5
我正在寻找一种方法来创建一个线性混合效应模型循环，它移除一个固定的效果自变量，运行模型，将它加回去，并从模型中删除另一个固定的独立变量。
例如：
library(lme4)
L1 <- lmer(cyl ~ disp + hp + mpg (carb|gear), mtcars)
如何创建该模型的循环，以便每个新的回归循环迭代删除单个固定效果自变量，将其添加到下一次迭代中，并删除不同的固定自变量？
这样做将产生以下模型：
原版全款：lmer(cyl ~ disp + mpg + (carb|gear), mtcars)
去除mpg的部分模型：lmer(cyl ~ disp + hp +
浏览 5提问于2022-02-09得票数 0
回答已采纳
我想使用lapply同时运行四个多级模型(使用lmer)。 使用带有一个因变量和一列自变量的lm()的一个简单示例如下： data(mtcars)
varlist <- names(mtcars)[3:6]
models <- lapply(varlist, function(x) {
lm(substitute(mpg ~ i, list(i = as.name(x))), data = mtcars)
}) 如何将其扩展为运行四个lmer()模型，每个模型都有不同的因变量和不同的自变量列表？对于所有四个模型，这两个级别将保持不变。四个(伪造的)示例模型是： data(m
浏览 18提问于2020-12-15得票数 2
回答已采纳
这是我的数据，我想做REML分析，看看随机因素的方差成分。
Plant<- rep(c("P1","P2","P3","P4"), each=9)
Leaves<- rep(rep(c("L1","L2","L3"), each=3),4)
Rep<-rep(c(1,2,3),12)
Ca<- c(3.280, 3.090, 3.185, 3.520, 3.600, 3.560, 2.880, 2.800, 2.840, 2.460, 2.440,
浏览 2提问于2021-07-30得票数 2
2回答
lmer模型的标准化系数
、、、
我曾经使用下面的代码来计算lmer模型的标准化系数。但是，随着新版本的lme，返回对象的结构发生了变化。
如何调整函数stdCoef.lmer使其与新的lme4版本一起工作？
# Install old version of lme 4
install.packages("lme4.0", type="both",
repos=c("http://lme4.r-forge.r-project.org/repos",
getOption("repos"
浏览 5提问于2014-08-05得票数 5
回答已采纳
我有一组数据，我已经用一个广义线性模型进行了分析，该模型在3向交互中有三个分类因素(factorA，factorB，factorC)，第四个连续因素(factorD)简单地添加到模型中。我试图从模型中获得一组Tukey字母组(即紧凑的字母显示)，但还没有找到一种方法来成功地包含交互。我对包含factorD不感兴趣，只对交互中的三个感兴趣。
我已经得到了Tukey调整后的成对比较：
lsmeans(my.glm, factorA*factorB*factorC)
但我不能想出如何从它产生一个紧凑的字母显示。这可以用multcomp包来完成，但我只能找到使用该包的主效果来做这件事的方法，而不是交互
浏览 4提问于2013-08-01得票数 1
我试图从lmer包中在R中执行一个lme4，但是我一直收到错误：
Paste0(“(”，vapply(findbar(F)，deparse1，"")，“)，未使用的参数(env = environment(f))”。
有什么想法吗？以虹膜为例：
library(lme4)
lmer(iris$Sepal.Length ~ iris$Sepal.Width + (1|iris$Species))
(我知道这在直觉上没有多大意义，但作为一个例子)
浏览 3提问于2022-03-11得票数 1
我试着比较两种使用多重估算建立的模型。当我试图比较这些模型时，小鼠的pool.compare()给出了一个错误:错误:对于类调用或错误的对象没有一目了然的方法：“fit1”和“fit0 0”的计算数不等，即使我使用的是相同的估算数据集。下面是一个可重复的例子：
library(mice)
library(miceadds)
library(lmerTest)
imp <- mice(nhanes, maxit = 2, m = 4)
summary(m0 <- pool(with(imp, lmerTest::lmer(bmi ~ 1 + (1 | chl)))))
summ
浏览 6提问于2020-01-16得票数 0
回答已采纳
1回答
我正在使用R中的lme4包来拟合我想在中介中使用的负二项模型(它比泊松模型更适合)。 我得到一个错误：“不支持结果模型的glmer系列” 为了巧妙起见，我简单地将类切换为lmer；例如，从 > class(model.Y)
[1] "glmerMod"
attr(,"package")
[1] "lme4" 至 > class(model.Ynew)
[1] "lmerMod"
attr(,"package")
[1] "lme4" 而且它是有效的。显然，这可能会产生问题，但我检查了残
浏览 15提问于2020-03-24得票数 0
1回答
当我使用lmer作为lme4作为na.action运行一个na.pass模型时，我会得到以下错误：
外部函数调用中的NA/NaN/Inf (第1条)
我像这样运行模型：
model1 <- lme4::lmer(agg_dv_singing ~ GMS.Musical.Training +
JAJ.ability + MDT.ability + MPT.ability + PDCT.ability +
PIAT.ability + agg_dv_long_note + demographics.age +
aggiv_entropy + aggiv_int
浏览 2提问于2021-06-06得票数 0
回答已采纳
我理解，在混合效应模型中使用连续变量或数值变量作为随机效应是没有什么意义的(例如，参见)。
但我想知道的是如果R中的lme4::lmer或nlme::lme故意阻止你这样做.
具体来说，我要问的是:如果我将lmer (或lme)任何非因素(非分类)变量作为随机效应提供，会自动将其作为一个因子来处理吗？
将factor()直接插入lmer (与使用lm时常用的方法一样)会产生以下错误：
lmer(y ~ z + (1|factor(x)), data = dat)
Error: couldn't evaluate grouping factor factor(x) within mode
浏览 1提问于2016-10-10得票数 8
回答已采纳
1回答
我正在尝试在lme4::lmer计划中加入一些drake模型，但是我得到了一个错误。
找不到“数据”，公式环境中缺少一些变量
如果我替换了一个lm模型，它就能工作。
这里是一个可重复的例子
library(drake)
library(lme4)
#> Loading required package: Matrix
#>
#> Attaching package: 'Matrix'
#> The following object is masked from 'package:drake':
#>
#>
浏览 3提问于2020-07-28得票数 2
回答已采纳
1回答
最近，我尝试在包中组织我的所有函数，以避免在每个脚本上有数百行代码。在运行了一组LME4多级模型之后，我编写了许多函数。这些函数的设计目的是生成很好的输出表(我知道有可用的包可以这样做，但是这些包通常不够灵活，无法根据我的需要定制表)。下面是一个函数的示例，它接受lmer模型的列表(例如，fm0、fm1、fm2)，并将固定效果参数及其各自的标准错误组合在一个输出表中(稍后我将使用该输出表与其他模型统计数据一起生成自定义回归表)。
#' @rdname f.regTableFix1
#' @title Fixed effects table for lmer models
#&#
浏览 0提问于2013-06-21得票数 0
回答已采纳
我一直在使用来自anova.merMod包的lme4函数，通过科学出版物的似然比测试获得固定效果的p值(大多数审查员仍然在我的领域中要求p值)。我注意到，anova.merMod函数使用REML = FALSE重新计算lmer函数(参见下面的示例)，这是一个非常好的特性，迫使不熟悉的用户正确地进行测试。但是，我一直试图阅读lme4包的大部分文档，但找不到这个特性的符号(例如，请参阅?anova.merMod，它将用户引导到?vcov.merMod)。这让我很困惑。
问：为什么在文档中没有明确提到这个特性？我是否理解错了？
Ps。在上似乎有一个关于这个问题的问题，但是答案使我更加困惑。
libr
浏览 4提问于2014-02-03得票数 2
回答已采纳
1回答
多结局二项回归R
、、、
可以使用cbind为普通lm回归指定多个结果：
set.seed(11)
df <- iris %>%
mutate(var1 = sample(c(0L,1L), 150, replace = TRUE),
var2 = sample(c(0L,1L), 150, replace = TRUE))
df %>%
lm(cbind(var1, var2) ~ Sepal.Width + Sepal.Length, data = .) %>%
coef()
# var1 var2
#
浏览 7提问于2022-10-10得票数 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
回答已采纳
我有data.frame，df
> str(df_ss)
'data.frame': 571 obs. of 4 variables:
$ final_grade : num 0.733 0.187 0.502 0.194 0.293 ...
$ time_spent : num -0.2 -0.326 -0.709 -0.168 -0.254 ...
$ gender_female: num 1 0 1 0 0 0 1 1 1 1 ...
$ course_ID : Factor w/ 26 levels "1","
浏览 6提问于2016-12-21得票数 1
回答已采纳
我想在一个lmer()循环中迭代地拟合一个for模型，并存储结果。当我遇到错误时，我不希望merMod对象(来自lmer()模型的输出)覆盖来自上一次迭代的merMod对象。示例：
# install.packages(c("lme4", "dplyr", "ggplot2"), dependencies = TRUE)
library("lme4")
library("dplyr")
library("ggplot2")
predList <- list()
j <- 1
for(i
浏览 0提问于2016-01-21得票数 0
回答已采纳
我正在使用lme4中的lmer函数分析一个庞大而复杂的数据集。我用点阵图和点阵图生成毛毛虫随机效应图。有没有办法用我的数据集中的一个因子或连续变量来对我的毛毛虫图进行颜色编码？我也使用过to数学，我可能需要帮助来理解如何使用“group”参数。这将有助于讨论要点。
library("lme4")
data(package = "lme4")
summary(grouseticks)
fit1<-lmer(TICKS~1+(1|LOCATION), grouseticks)
rr1<-ranef(fit1, condVar = TRUE)
dotp
浏览 5提问于2015-01-05得票数 4
回答已采纳
1回答
我使用R中的lmer()函数(lme4包)来分析一项纵向研究，在其中我测量了120名受试者，6次。首先，我指定了这样一个模型：
library(lme4)
model1 = lmer(DV ~ 1 + X1*X2 + (1+X1|SubjectID), REML="false")
X1是一个时变变量(级别-1)，而X2是一个主题级变量(级别-2).
因为这些主题嵌套在几个团队中，所以我被建议在团队级别(级别-3)中包含一个随机拦截。然而，我只发现如何同时包含随机拦截和斜率：
model2 = lmer(DV ~ 1 + X1*X2 + (1+X1|TeamID/Subjec
浏览 3提问于2016-11-12得票数 3
回答已采纳
1回答
我试图用而不是nlme模型重现lmer的绘图，但是predict函数出现了一个错误：
library(nlme)
library(lme4)
library(ggplot2)
fm2 <- lmer(distance ~ age + Sex + (1|Subject), data = Orthodont)
newdat <- expand.grid(Sex=unique(Orthodont$Sex),
age=c(min(Orthodont$age),
max(Orthodont$age
浏览 4提问于2022-02-23得票数 0
回答已采纳
1回答
我有以下研究，我想用混合效应模型来分析：
“受试者”分为两组(治疗A和B)。
“体重”记录在治疗前和治疗后3个月(“时间”)(重复措施)。
也需要对主语“年龄”和“性别”进行修正。
的主要问题是：两组对体重的影响是否不同？
对于混合效果，我考虑使用lmer函数的lme4包语法：
lmer(weight ~ Group*Time + age, (1|subject) + (1|gender), data=mydata)
这个语法是正确的，还是我需要使用更复杂的术语，如下面给出的术语：
(time|subject)
(time + 1|subject)
(1|subject) + (1|Group
浏览 4提问于2020-07-10得票数 3
回答已采纳
我在跟进this great answer的事。下面的foo函数获取VarCorr(fit)输出的Name列，并将它们作为summary(rePCA(fit))调用的列名。 当我们输入fm1，fm2时，它工作得很好，但我想知道为什么它在fm3上失败了？有什么解决办法吗？ library(lme4)
dat <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/sng.csv')
fm1 <- lmer(diameter ~ 1 + (1|plate) + (1|sample), Pe
浏览 33提问于2020-10-07得票数 1
回答已采纳
1回答
各位，我正在打印lmer() {nlme}和lme() {lme4}模型的{texreg}表，其中包括方差。然而，这两个模型之间的差异很大(几个数量级)。lme()方差似乎是lmer()方差的平方根。哪些是正确的？
library(plm)
library(lme4)
library(nlme)
library(texreg)
data("Grunfeld", package="plm")
reML0 <-lmer(inv ~ value + capital + (1|firm), data=Grunfeld)
reML1 <- lme(in
浏览 3提问于2015-08-18得票数 1
回答已采纳
2回答
我试图在一个循环中使用来自nlme包的lapply函数。这适用于来自lmer包的函数，但会为lme生成错误消息。我如何循环lme函数，类似于下面示例中的lmer函数？
library("nlme")
library("lme4")
set.seed(1)
dt <- data.frame(Resp1 = rnorm(100, 50, 23), Resp2 = rnorm(100, 80, 15), Pred = rnorm(100,10,2), group = factor(rep(LETTERS[1:10], each = 10)))
## Synt
浏览 1提问于2014-10-14得票数 2
回答已采纳
1回答
library("sjPlot", lib.loc="~/R/win-library/3.2")
Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) :
there is no package called ‘pbkrtest’
Error: package or namespace load failed for ‘sjPlot’
> sjp.lmer(fitted2, type = "fe")
Error: c
浏览 6提问于2016-02-20得票数 0
回答已采纳
我想绘制一个lme4::lmer模型的固定效果的估计值。目前，sjPlot::plot_model输出包括两个我不想要的额外绘图: SD (观察值)和SD (截距)。如何删除这两个额外的绘图？ library(lme4)
data(iris)
lmerMod <- lmer(Sepal.Length ~ Sepal.Width * Petal.Length + (1|Species), data=iris)
library(sjPlot)
plot1 <- plot_model(lmerMod,
type = "est",
浏览 72提问于2021-11-23得票数 1
回答已采纳
1回答
lme4模型中异常值的识别
、、、、
我想识别lmer模型中的异常值(lme4包)。我是LMERConvenienceFunctions不感兴趣移除它们(包所做的)-我只是想看看列出来的离群点。
我使用的一个模型示例：
model1<-lmer(Value~ Moisture + Planting + (day|plot), data=plants1)
浏览 2提问于2017-08-16得票数 0
我正在尝试使用mediation包来适应一个基于GLMM回归的中介模型，该模型使用了lme4包：
library(mediation)
library(lme4)
curvature_full_model <- lmer(scale(curvature_sqrt) ~ scale(edit_session_length) + scale(edit_session_duration_log) + scale(edit_session_number_sqrt) + scale(role_within_artifact_group_ratio) + scale(role_within_com
浏览 0提问于2018-06-14得票数 0
我对和线性混合模型有一个问题。当我试图用coff.glmulti函数估计模型平均系数时，我得到了以下错误：
Data.frame中的错误(.，check.names = FALSE)：参数意味着不同的行数: 1，0
我做了一些调试，发现问题从突出显示的coeff.glmulti函数行开始：
if (length(object@adi) >= 1)
for (j in 1:length(object@adi)) {
cak[[length(names(cak)) + 1]] = object@adi[[j]
浏览 2提问于2013-12-06得票数 1
回答已采纳
1回答
从lme4模型中获取预测器的类
、、、、
在对lme4模型进行拟合之后，我想知道如何从terms(fit)[[3]]获得预测器的类。
这里有一个简单的例子，但我很欣赏lme4中任何其他模型的功能答案。
注意：所有东西都必须从模型中提取出来。
library(lme4)
h <- read.csv('https://raw.githubusercontent.com/hkil/m/master/h.csv')
h$year <- as.factor(h$year)
m <- lmer(scale~ year*group + (1|stid), data = h)
terms(m)[[3]] ##
浏览 3提问于2020-05-21得票数 1
回答已采纳
1回答
我正在调整一个混合效应模型，由于观测到的异方差，有必要包含一个效应来容纳它。因此，使用lme包的nlme函数很容易解决这个问题，请参阅下面的代码：
library(nlme)
library(lme4)
Model1 <- lme(log(Var1)~log(Var2)+log(Var3)+
(Var4)+(Var5),
random = ~1|Var6, Data1, method="REML",
weights = varIdent(form=~1|
浏览 2提问于2021-04-21得票数 1
回答已采纳
假设我们假设了以下线性混合模型，我们一般称之为fit。
library(lme4)
fit <- lmer(Reaction ~ Days + (1 | Subject), data = sleepstudy)
假设我们也对从fit中提取原始公式(即固定效应和随机效应部分)感兴趣，我们想把它传递给第二个称为fit2的LMM。显然，将terms(fit)传递给lmer()的formula参数不起作用。
> fit2 <- lmer(terms(fit), data = sleepstudy)
Error: No random effects terms specified i
浏览 7提问于2022-01-11得票数 6
回答已采纳
1回答
随机模型公式对象
、、
我想把公式放在随机模型中，但我认为下面的错误是由于公式对象(?)错误造成的，但无法修复它。
set.seed(1234)
mydata <- data.frame (A = rep(1:3, each = 20), B = rep(1:2, each = 30),
C = rnorm(60, 10, 5))
mydata$A <- as.factor(mydata$A)
mydata$B <- as.factor(mydata$B)
myfunction <- function (mydata, yvars, genovar,
浏览 2提问于2011-11-15得票数 1
回答已采纳
在用lme4::lmer拟合混合效应模型时，我遇到了以下错误：
(函数(cl，name，valueClass)中的错误:类“dgTMatrix”中的对象的赋值对于“dgTMatrix”对象中的@'Dim‘无效；is (值，“整数”)不正确
我执行了谷歌搜索，但找不到任何结果，其中包含了整个错误信息。通过将我的代码恢复到一个工作状态，我终于找到了原因。我已在下文提供决议的答复。
从错误信息来看，原因很难推断，但是错误很容易发生，我在这里记录错误的原因，作为将来的参考。
错误发生在以下代码中：
library(lme4)
data(sleepstudy)
fm1 <- lme
浏览 0提问于2018-01-26得票数 1
回答已采纳
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM
腾讯云开发者
