Ich habe in ggplot2 einen Interaktionsplot erstellt. Dieser vergleicht die Mortalität von Käfern welche auf Bt und nicht-Bt (Non) Pflanzen gelbt haben und an pinheads (pins) oder nicht-pinheads gegessen haben (siehe Angehängte Grafik).
Nun würde ich gerne an jedem der 4 Punkte in der Grafik ein Mass für die Varianz angeben. Am besten einen kleinen Boxplot oder einen Balken der den SE anzeigt. Das Problem ist, das meine response variable (dead.id) binär ist (0/1).
Kann mir jemand weiterhelfen?
Vielen Dank,
Michi
Untenstehend findet ihr meinen Code sowie meinen Datensatz (dput)
CODE:
Code: Alles auswählen
#INTERACTION PLOT
ggplot(data = dta,
aes(x = plant, y = dead.id, colour = pin, group=pin)) +
stat_summary(fun.y=mean, geom="point")+
stat_summary(fun.y=mean, geom="line")
dput(dta)
DATEN:
> dput(dta)
structure(list(sample = structure(c(1L, 1L, 1L, 2L, 3L, 3L, 3L,
4L, 4L, 4L, 5L, 6L, 6L, 6L, 7L, 7L, 8L, 8L, 8L, 9L, 9L, 9L, 10L,
10L, 10L, 11L, 11L, 12L, 12L, 12L, 13L, 13L, 13L, 14L, 14L, 14L,
15L, 15L, 16L, 16L, 16L, 17L, 17L, 17L, 18L, 18L, 18L, 19L, 19L,
19L, 20L, 20L, 21L, 22L, 22L, 22L, 23L, 23L, 23L, 24L, 24L, 24L,
25L, 26L, 26L, 26L, 27L, 27L, 27L, 28L, 28L, 29L, 29L, 29L, 30L,
30L, 30L, 31L, 31L, 31L, 32L, 32L, 32L, 33L, 33L, 34L, 34L, 34L,
35L, 35L, 35L, 36L, 36L, 36L, 37L, 37L, 37L, 38L, 38L, 38L, 39L,
39L, 39L, 40L, 40L, 40L, 41L, 41L, 42L, 42L, 42L, 43L, 44L, 44L,
44L, 45L, 45L, 45L, 46L, 46L, 47L, 48L, 48L, 48L, 49L, 49L, 50L,
50L, 50L, 51L, 51L, 51L, 52L, 52L, 53L, 53L, 53L, 54L, 54L, 54L,
55L, 55L, 55L, 56L, 56L, 56L, 57L, 57L, 57L, 58L, 58L, 58L, 59L,
59L, 59L, 60L, 60L, 60L, 61L, 61L, 61L, 62L, 62L, 62L, 63L, 63L,
63L, 64L, 64L, 64L, 65L, 65L, 65L, 66L, 66L, 67L, 67L, 67L, 68L,
68L, 68L, 69L, 69L, 69L, 70L, 70L, 70L, 71L, 71L, 71L, 72L, 72L,
72L, 73L, 73L, 73L, 74L, 74L, 74L, 75L, 75L, 75L, 76L, 76L, 76L,
77L, 77L, 77L, 78L, 78L, 78L, 79L, 79L, 79L, 80L, 80L, 81L, 81L,
81L, 82L, 82L, 82L, 83L, 83L, 84L, 84L, 84L, 85L, 85L, 85L, 86L,
86L, 86L, 87L, 87L, 87L, 88L, 88L, 89L, 89L, 89L, 90L, 90L, 90L,
91L, 91L, 91L, 92L, 92L, 93L, 93L, 93L, 94L, 94L, 94L, 95L, 95L,
95L, 96L, 96L, 97L, 97L, 97L, 98L, 98L, 99L, 99L, 99L, 100L,
101L, 101L, 101L, 102L, 102L, 102L, 103L, 103L, 103L, 104L, 104L,
104L, 105L, 105L, 105L, 106L, 106L, 106L, 107L, 107L, 107L, 108L,
108L, 108L, 109L, 109L, 109L, 110L, 110L, 110L, 111L, 112L, 112L,
112L, 113L, 113L, 113L, 114L, 114L, 114L, 115L, 115L, 115L, 116L,
116L, 116L, 117L, 117L, 118L, 118L, 118L, 119L, 119L, 119L, 120L,
120L, 120L, 121L, 121L, 121L, 122L, 122L, 122L, 123L, 123L, 124L,
124L, 125L, 125L, 126L, 126L), .Label = c("10BtC", "10BtCat",
"10BtJa", "10nonC", "10nonCat", "10nonJa", "11BtC", "11BtCat",
"11nonC", "11nonCat", "11nonJa", "12BtCat", "12BtJa", "12nonC",
"12nonCat", "12nonJa", "13BtCat", "13nonC", "14BtC", "14BtCat",
"14nonJa", "15BtC", "15nonC", "15nonCat", "15nonJa", "16BtC",
"16BtCat", "16BtJa", "16nonC", "16nonCat", "16nonJa", "17BtC",
"17BtJa", "17nonC", "17nonCat", "18BtC", "18BtCat", "18nonCat",
"18nonJa", "19BtC", "19BtCat", "19BtJa", "19nonC", "19nonCat",
"19nonJa", "1BtC", "1BtJa", "1nonC", "1nonCat", "1nonJa", "20BtC",
"20BtJa", "20nonCat", "20nonJa", "21BtCat", "21BtJa", "21nonC",
"21nonCat", "21nonJa", "22BtC", "22BtCat", "22BtJa", "22nonC",
"22nonCat", "22nonJa", "23BtC", "23BtCat", "23BtJa", "23nonC",
"23nonCat", "23nonJa", "24BtC", "24BtJa", "24nonC", "24nonCat",
"24nonJa", "25BtC", "25BtCat", "25BtJa", "25nonC", "25nonCat",
"25nonJa", "26nonC", "2BtC", "2BtCat", "2BtJa", "2nonCat", "2nonJa",
"3BtC", "3BtCat", "3BtJa", "3nonCat", "3nonJa", "4BtC", "4BtCat",
"4BtJa", "4nonCat", "4nonJa", "5BtCat", "5BtJa", "5nonC", "5nonCat",
"5nonJa", "6BtC", "6BtCat", "6BtJa", "6nonC", "6nonCat", "6nonJa",
"7BtCat", "7BtJa", "7nonC", "7nonCat", "7nonJa", "8BtC", "8BtCat",
"8BtJa", "8nonC", "8nonCat", "8nonJa", "9BtC", "9BtCat", "9BtJa",
"9nonC", "9nonCat", "9nonJa"), class = "factor"), id = structure(c(1L,
1L, 1L, 2L, 3L, 3L, 3L, 4L, 4L, 4L, 5L, 6L, 6L, 6L, 1L, 1L, 2L,
2L, 2L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 2L, 2L, 2L, 3L, 3L, 3L,
4L, 4L, 4L, 5L, 5L, 6L, 6L, 6L, 2L, 2L, 2L, 4L, 4L, 4L, 1L, 1L,
1L, 2L, 2L, 6L, 1L, 1L, 1L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 1L, 1L,
1L, 2L, 2L, 2L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 1L,
1L, 1L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 1L, 1L, 1L, 2L, 2L, 2L,
5L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 4L, 5L,
5L, 5L, 6L, 6L, 6L, 1L, 1L, 3L, 4L, 4L, 4L, 5L, 5L, 6L, 6L, 6L,
1L, 1L, 1L, 3L, 3L, 5L, 5L, 5L, 6L, 6L, 6L, 2L, 2L, 2L, 3L, 3L,
3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L, 2L,
3L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 2L, 2L,
2L, 3L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L,
3L, 3L, 3L, 4L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L,
2L, 2L, 3L, 3L, 3L, 4L, 4L, 5L, 5L, 5L, 6L, 6L, 6L, 4L, 4L, 1L,
1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 5L, 5L, 5L, 6L, 6L, 1L, 1L, 1L,
2L, 2L, 2L, 3L, 3L, 3L, 5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L,
2L, 3L, 3L, 5L, 5L, 5L, 6L, 6L, 2L, 2L, 2L, 3L, 4L, 4L, 4L, 5L,
5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L,
4L, 5L, 5L, 5L, 6L, 6L, 6L, 2L, 2L, 2L, 3L, 4L, 4L, 4L, 5L, 5L,
5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 4L, 4L, 4L, 5L,
5L, 5L, 6L, 6L, 6L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 4L, 4L, 5L,
5L, 6L, 6L), .Label = c("BtC", "BtCat", "BtJa", "nonC", "nonCat",
"nonJa"), class = "factor"), plant = structure(c(1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
2L), .Label = c("Bt", "Non"), class = "factor"), pin = c(0L,
1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L,
1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L,
1L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L,
1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L,
0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L,
1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L,
0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L,
0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,
1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 1L,
0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L), dead.id = c(1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L,
1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L,
1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L,
0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L,
1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L,
1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L,
1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L,
1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L,
0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L,
0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 0L,
1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L,
0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L,
0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,
0L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
0L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 1L,
1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L)), .Names = c("sample",
"id", "plant", "pin", "dead.id"), class = "data.frame", row.names = c(NA,
-340L))