[英]Interaction plot with confidence intervals
我想繪制一個帶有置信區間的交互圖,用於明天到期的分配。 這是我的數據集:
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數據執行類似底部的操作 ,但我顯然做錯了,因為我只收到錯誤消息。
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
顯然,我不知道我在做什么,所以有人可以幫我糾正我的代碼嗎?
使用第一個aggregate
命令創建對象df
時,包含因變量的列不再稱為total_migration
而是x
。
因此,如果在ggplot
命令中用x
替換total_migration
,則代碼可以正常工作:
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))
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.