我很难理解如何使用map来应用一个函数来分析我数据中各个主题的数据library(tidyverse)这是我开始使用的模型我可以查看个人的数据：我可以获取
浏览 16提问于2017-07-13得票数 1
1回答
我使用R中的lmer()函数(lme4包)来分析一项纵向研究，在其中我测量了120名受试者，6次。首先，我指定了这样一个模型：model1 = lmer(DV ~ 1 + X1*X2 + (1+X1|SubjectID), REML="false")
X1是一个时变变量(级别-1)，而X2是一个主题级变量(级别-2).因为这些主题嵌套在几个团队中，所以我被建议在团队级别(级别-3)中包含一个随机拦截。
浏览 3提问于2016-11-12得票数 3
回答已采纳
在变量中REE ~ Time*predictor + (1 + Time | Case identifier)library(tidyverse)library(sjPlot) lmer(REE ~ Day*(x) + (1 + Day | Studynumber= lme
浏览 5提问于2022-09-23得票数 0
1回答
R多模型多函数
、、、、
我编写了一个例程，从lmer模型中提取信息来计算ICC，并从lmerTest的ranova函数中获得LRT。下面的内容可以工作，但我怀疑它可以通过以下方法得到改进：(a)将两个函数组合为一个函数并返回一个列表，但我似乎无法使用purrr的map函数访问list元素，(b)使用多个mutate/purrr行在一个位置获取所有所需的数据我的代码使用Hox (2002)中提供的"Peet“数据集，可在UCLA IDRE网站上查阅：
library(fore
浏览 0提问于2019-01-28得票数 4
1回答
我有一个长格式的数据集，包含200个变量，94个主题，每个主题对每个变量有1到3个度量。有人建议，在估算缺失值与被测对象的平均值之前，应使用重复测量方差或混合模型，以确定测量的重复性。m1 <- lme(var1 ~ measurement, random=~1|ID, data=mydata,
浏览 3提问于2021-03-19得票数 2
回答已采纳
考虑如下的数据集： library(lme4)head(sleepstudy, 10)
1 249.5600我可以这样做： library(dplyr)library(purrr)
group_by(grp = gl(2, n()/2)) %>%nest() %>%
mutate(fit = map(data, ~ lmer(Reaction
浏览 20提问于2020-01-11得票数 2
回答已采纳
我们有一个包含26个品牌的大型数据集，在93家门店销售，在399周内销售。这些品牌仍然分为子品牌(f.ex.：brand =高露洁)，但子品牌(556)仍然存在:高露洁优质白酒/高露洁额外等。我们为每个子品牌计算了每周商店级别的品牌共享价格:计算：(每个子品牌和每个商店每周的每盎司移动)除以(子品牌每盎司移动的总和，每个单独商店每周引用一个品牌)*(每周商店级别每个子品牌的日志价格)
一切都成功了！我们创建了一个包含所有详细计算的数据框架(data = tooth4)我们的最终兴趣是运行线性回归来预测价格对移动变量的影响-->现在的问题是销售变量(一个虚拟变量，表示特定商店中的特定子品牌
浏览 2提问于2019-05-15得票数 0
在afex软件包中，我们可以找到这个方差分析的例子：# estimate mixed ANOVA on the full designbetween = c("treatment", "gender"), 我的问题是，如何在lme4如果我只是写
lme
浏览 2提问于2016-05-08得票数 5
回答已采纳
我有一个通过特定的lme4::lmer()模型语法运行变量的函数，如下所示： run_my_model <- function(i)
f <- formula(paste("y_variable()一起使用时，它对数据集中的所有变量都有效： test <- map(names(my_dataset), run_my_model) 返回 Error in eval_f(x, ...) : DowndatedVtV is not positive definite 如果我<
浏览 158提问于2020-08-25得票数 0
1回答
三级部分嵌套模型
、、、
我正在用R和lme4模拟随时间变化的群体心理治疗对象。我的数据有以下结构：
lmer(outcome ~ time * treatment + (time | subject), data=data, REML=FALSE)
现在我想知道我是否应该使用一个三层部分嵌套模型纵向三级多水平模型的先验幂<e
浏览 1提问于2016-05-13得票数 5
回答已采纳
我使用ezANOVA实现了一个实验设计的分析，实验设计中包含一个主题变量和一个主题变量之间。我还用了lmer来检查残差似乎很直接。plot(fitted(model_lmer), residuals(model_lmer))
然而，由于ezANOVA也实现了球度的测试和修正，我想坚持它，并找到一种方法来检查假设，再正规的残差。
浏览 2提问于2016-07-24得票数 5
回答已采纳
为了进行功率分析，我想改变我的模型的固定效应系数，以测试不同影响大小的功率。在主SIMR ()中，在使用以下代码模板运行模型之后，我应该能够修改固定的效果：然而，当我试图用我自己的模型来做这件事时，我会收到以下错误：下面是一个可重复的例子：library(lme4)mod<-lmer(Sepal.Lengthcode work
浏览 1提问于2020-08-09得票数 4
回答已采纳
2回答
如何制作混合模型的森林样地及其置信区间。我试过这段代码random = ~1|ID, na.action = na.omit, method = "ML", plot_summs (model)
浏览 30提问于2022-06-05得票数 1
回答已采纳
2回答
多元回归模型滤波估计的点须图
、、、、
into female and male filter(sex=="F") filter(sex=="M")
#poisson regression
ab_f_LBS= lmer(LBS ~ ft + grid + (1|byear), data = subset(femaletol=1e-6,optimizer="bobyqa",optCtr
浏览 0提问于2018-10-16得票数 5
回答已采纳
1回答
我的目的是比较主题内比较(patient1比较每个变量(variable1.)在不同的时间点(00、66和01)之间。我已将数据修改为纵向数据，以便有时间作为一个因素，但我不能使用方差分析(记住，分析必须是对的！)do(tidy(t.test(.它起作用了第二进场dat %>%
select(paciente, matchesf
浏览 2提问于2022-01-17得票数 -2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
腾讯会议
云直播
对象存储活动推荐
腾讯云开发者
