程序员经验分享-hwhale

在 lmer 输出上运行汽车方差分析后出现列表类型错误

2023-11-29

我正在尝试在 R 中做一个随机效应模型。我想运行一个anova在输出中，我浏览了一些教程并复制了他们的示例，但我的不起作用。我在 Anova 上不断收到以下错误：

(1:length(names))[-which.term][sapply(names[-which.term], function(term2) is.relative(term, : 无效的下标类型“list”

我正在运行的代码如下，任何人对问题所在有任何想法将不胜感激！

library(lme4)
library(car)

div1<-rep(1,5)
div2<-rep(2,3)
div3<-rep(3,3)
div6<-rep(6,6)
div<-c(div1,div2,div3,div6)
res<-c(3.8082479,7.7819745,3.3792467,7.2288647,3.4564646,1.8043898,5.1443293,3.9467614,2.5922306,1.9996585,4.2004104,0.7290807,2.1854365,3.4118980,3.2464388,2.9607496,1.9993038)
df<-data.frame(div,res)

randeff<-lmer(res~1+(1|div),data=df,REML=FALSE)
summary(randeff)
Anova(randeff)

car::Anova专为分析而设计固定效应仅有的。您试图分析一个没有随机效应的模型，因此它不高兴（不幸的是，错误消息难以理解）。要检验此假设，请向数据集添加随机列，并向模型添加固定效应：

> df$y <- rnorm(nrow(df))
> randeff2 <- lmer(res~y+(1|div),data=df,REML=FALSE)
> Anova(randeff2)
Analysis of Deviance Table (Type II Wald chisquare tests)

Response: res
   Chisq Df Pr(>Chisq)
y 0.0418  1      0.838
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)

r

lme4

ANOVA

在 lmer 输出上运行汽车方差分析后出现列表类型错误 的相关文章

  • 绘图不适合 R 窗口

    我正在使用metaR 中的包对比例进行荟萃分析 我的分析中有大约一百项研究 当我打字时forest result 创建我的森林图 它可以工作 但图形在顶部和底部被切掉 无论我做什么 它都不适合绘图窗口 我对如何解决这个问题感到困惑 有任何想

  • 保存/加载 data.table 的最快方法

    我想做的实际上是使用最快的可用方法来存储data table以便进一步处理 大致如下 从 CSV RDS 读取原始数据 将其转换为data table 将其保存为针对重新读取而优化的格式 RDS 似乎不适用于data table 是对的吗

  • 展平 ggplot2 中的 alpha 通道

    ggplot2 的某些对象 函数 支持 alpha 通道 http docs ggplot2 org current geom point html http docs ggplot2 org current geom point html

  • lme4：如何指定 2 个与随机截距的相关性，而不添加随机斜率之间的相关性

    重新发布自stats stackexchange com https stats stackexchange com q 195385 33560 我试图在 R 的 lme4 包中指定一个模型 其中随机截距和随机斜率之间有 2 个相关性 但

  • 在 R/ggplot2 中将字符串转换为函数参数的最佳方法？ [复制]

    这个问题在这里已经有答案了 我正在开发一个闪亮的应用程序 用户可以选择可以使用 ggplot2 绘制哪些变量 但是我完全不确定将字符串 即要绘制的变量的名称 转换为合适的函数参数的最佳方法 考虑以下非常人为的 有效的示例 df lt dat

  • R：使用 dcast 时包含没有条目的因子

    我在数据帧上使用 reshape2 函数 dcast 其中一个变量是某些级别未出现在数据框中的因素 但我会将所有值包含在创建的新列中 例如 假设我运行以下命令 library reshape2 dataDF lt data frame id

  • 抑制 r markdown 中的控制台输出，但保留绘图

    嗨 我有以下降价块 r echo FALSE warning FALSE message FALSE error FALSE lapply obj function x plot x main some plot box axis 1 at

  • 使用梯度下降（最速下降）估计线性回归

    示例数据 X lt matrix c rep 1 97 runif 97 nrow 97 ncol 2 y lt matrix runif 97 nrow 97 ncol 1 我已经成功创建了成本函数 COST lt function th

  • 如何使用 ggplot 绘制矩阵图

    我想可视化一个矩阵 MAT lt matrix c 100 7 0 0 49 0 0 0 49 nrow 3 ncol 3 gt MAT 1 2 3 1 100 7 0 2 0 49 0 3 0 0 49 然而 标准方法不能正确地对小数字进

  • 没有垂直线的直方图

    当我创建直方图时 它看起来很像这样 set seed 1 x lt 1 100 y lt x rnorm 50 y round y hist y Is there a way to make a histogram look a bit l

  • 两个闪亮的操作按钮

    我正在编写一个闪亮的函数 其中包含两个操作按钮 这两个按钮是左右按钮 单击时可以帮助绘图移动 第一次单击时 这两个按钮都工作正常 当我重新单击它们时 就会出现问题 无论我单击哪个按钮 它都会在两个范围之间来回弹跳 我猜可能是 右 和 左 按

  • 如何在闪亮的应用程序中显示矩阵，并用条件指定颜色？

    我有一个具有正值和负值的矩阵 M 我正在尝试使用 DT 包在闪亮的应用程序中显示为表格 我想用不同的颜色显示矩阵 红色的正数和负数 例如 到目前为止 我只能以一对一的方式添加颜色 但我想以这种方式添加颜色 如果matrix values g

  • Sweave 缓存包

    我正在尝试编写一份报告 我的问题是每次我编译 R 时都会加载我在报告中使用的包 如 ggplot2 MASS cubature 这是非常耗时的 有没有办法查包裹 I found 缓存编织但它不起作用 这是我在 sweave 文件中添加的块

  • 插值时间序列

    我有两组具有不同时间戳的数据 一组数据包含校准数据 另一组包含样本数据 校准的频率比样品少得多 我想做的是将校准数据 低频 插值到采样时间序列 高频 上 sam lt textConnection time value 01 00 52 2

  • 在 R 中使用科学记数法和 xtable

    我将 data frame 传递给 xtable dat table lt xtable dat 1 20 digits 10 我宁愿使用科学记数法 而不是像那样显示数字 我该怎么做呢 看过了 但我发现的只是R 格式化xtable中的数字

  • 如何在 ggplot2 中向 x 轴添加特定值？

    我正在尝试在 ggplot2 中绘制图表 我希望 x 轴显示 2 84 以及下面键入的序列 除了在 Breaks 中输入所有精确值之外 还有其他方法吗 我尝试了谷歌 但它没有解决我的问题 scale x continuous limits

  • ggplot2：图例中的斜体

    我正在尝试编辑图例中的标签 以便第一个标签 WT 为纯文本 而后续 7 个标签为斜体 我一直在使用element text face c plain rep italic 7 但这导致没有任何标签被转换为斜体 我有点困惑为什么它不起作用 因

  • 如何对plot_ly()图表进行分面？

    Using ggplot2 and plotly制作交互式散点图facet wrap library ggplot2 library plotly g lt iris gt ggplot aes x Sepal Length y Sepal

  • R 中的频率加权，与 Stata 的结果比较

    我正在尝试分析明尼苏达大学 IPUMS 数据集中的数据1990 年美国人口普查 http usa ipums org usa sampdesc shtml us1990a in R 我正在使用survey http faculty wash

  • 如何将带有几行代码的字符数组转换为 data.frame？

    我有以下数组 my list lt c Jan 01 Dec 31 00 00 24 00 Jan 01 Jun 30 12 00 18 00 Jul 06 Dec 31 09 00 19 00 导致以下结果的最短代码是什么 x1 x2 x

随机推荐

热门标签