r - 具有置信区间的交互图

标签 r ggplot2 interaction confidence-interval

我想为明天到期的作业绘制一个带有置信区间的交互图。这是我的数据集:

Snail0

    water   lightness   total_migration
    control dark    9.5
    control dark    20
    control dark    22
    control dark    24
    control dark    10
    control dark    18
    control dark    19
    control dark    19.5
    control dark    23
    control dark    2
    control dark    10.5
    control dark    22.5
    control dark    26
    control dark    32
    control dark    2
    control dark    12.5
    control dark    13.5
    control dark    16.5
    control dark    19
    control dark    20
    control dark    21
    control dark    22
    control dark    24
    control dark    25
    control light   11
    control light   20.5
    control light   22.5
    control light   24
    control light   21
    control light   23
    control light   23.5
    control light   28
    control light   30
    control light   18.5
    control light   19.5
    control light   19.5
    control light   23.5
    control light   24
    control light   2
    control light   20
    control light   23
    control light   28.5
    control light   29
    control light   16
    control light   17.5
    control light   26.5
    control light   28.1
    control light   30.5
    crab    light   6
    crab    light   14
    crab    light   15
    crab    light   17
    crab    light   19
    crab    light   7
    crab    light   14
    crab    light   17
    crab    light   22
    crab    light   26
    crab    light   4.5
    crab    light   10
    crab    light   11
    crab    light   21
    crab    light   24
    crab    light   14
    crab    light   20
    crab    light   22.5
    crab    light   23
    crab    light   14
    crab    light   16
    crab    light   18.5
    crab    light   22
    crab    light   23
    crab    dark    14
    crab    dark    19
    crab    dark    20
    crab    dark    23
    crab    dark    25
    crab    dark    9
    crab    dark    18
    crab    dark    18
    crab    dark    20
    crab    dark    21
    crab    dark    15
    crab    dark    15
    crab    dark    15
    crab    dark    19
    crab    dark    23
    crab    dark    12
    crab    dark    14.5
    crab    dark    21
    crab    dark    21.5
    crab    dark    22
    crab    dark    14.5
    crab    dark    16
    crab    dark    19
    crab    dark    20
    crab    dark    22.5
    snail   dark    6
    snail   dark    15
    snail   dark    20
    snail   dark    22
    snail   dark    15
    snail   dark    15
    snail   dark    16
    snail   dark    17
    snail   dark    21
    snail   dark    15
    snail   dark    20
    snail   dark    15
    snail   dark    16
    snail   dark    19.5
    snail   dark    23
    snail   dark    24
    snail   dark    4
    snail   dark    18.5
    snail   dark    20
    snail   dark    24
    snail   dark    25
    snail   light   13
    snail   light   13
    snail   light   13
    snail   light   18
    snail   light   15
    snail   light   16
    snail   light   17
    snail   light   17
    snail   light   19
    snail   light   14
    snail   light   14
    snail   light   19
    snail   light   7
    snail   light   13
    snail   light   13
    snail   light   13
    snail   light   16
    snail   light   3
    snail   light   9
    snail   light   15
    snail   light   22
    snail   light   23
interaction.plot(Snail0$water,Snail0$lightness,Snail0$total_migration) 

我一直在尝试使用 ToothGrowth data 来做一些类似于底部所做的事情。 ,但我显然做错了什么,因为我只收到错误消息。
Snail0$water <- factor(Snail0$water)


df<-with(Snail0, aggregate(total_migration, list(lightness=lightness, water=water), mean))


df$se<-with(Snail0, aggregate(total_migration, list(lightness=lightness, water=water), function(x) sd(x)/sqrt(10)))[,3]


gp<-ggplot(df, aes(x=water, y=total_migration, colour=lightness, group=lightness))

gp+geom_line(aes(linetype=lightness)) + geom_point(aes(shape=lightness)) + geom_errorbar(aes(ymax=x+se, ymin=x-se))

错误信息:
Error: Aesthetics must either be length one, or the same length as the dataProblems:lightness, water, lightness, lightness

显然,我不知道我在做什么,所以有人可以帮我更正我的代码吗?

最佳答案

创建对象时 df与第一 aggregate命令,包含因变量的列不再称为 total_migration但是 x .

因此,如果您替换 total_migration,您的代码就可以工作。与 xggplot命令:

gp <- ggplot(df, aes(x=water, y=x, colour=lightness, group=lightness))
gp + geom_line(aes(linetype=lightness)) + 
     geom_point(aes(shape=lightness)) + 
     geom_errorbar(aes(ymax=x+se, ymin=x-se))

enter image description here

关于r - 具有置信区间的交互图,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/13453602/

相关文章:

r - 从函数返回是否会创建一个额外的副本?

R:处理有时会导致 R session 崩溃的函数?

r - ggplot2 包含位图的自定义注释,具有反转的 y 比例

r - 子图中共享图例

r - 从 R 中的同一线性模型中提取特定的相互作用系数

c++ - Python 和 C++ 脚本交互

r - 如何使用 ggplot2 + directlabels 为标签使用自定义名称

r - 使用类似于 xpath 的方式访问 R 中列表的元素

r - ggplot2 条形图的多个子组

wpf - 在特定键上使用 EventTrigger