#r #ggplot2 #grouping #line
#r #ggplot2 #группировка #строка
Вопрос:
Я строю значения на основе взаимной пересадки (origin -> new, origin -> origin) для 2 сайтов. Я хочу отформатировать тип линии для этих значений на основе истории, где origin -> origin — сплошная линия, а origin -> new — пунктирная линия. Для 2 сайтов значения для «origin» и «new» разные. Ниже приведено подмножество данных:
structure(list(FragID = c("1004", "1006", "1015", "1038", "1087",
"1089", "1107", "1116"), ParentID = c("166", "166", "166", "166",
"144", "144", "144", "144"), ParentSite.x = c("Inner Bay", "Inner
Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer
Bay", "Outer Bay"), FragSite.x = c("Inner Bay", "Outer Bay", "Inner
Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer
Bay"), TotalSA = c(0.0171970755726674, 0.0338197513082082,
0.0215722601402604, 0.030712272182997, 0.027529366126288,
0.029650482611575, 0.0316984120058258, 0.0356299659679559), T03SA =
c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665,
0.000966259734683879, 0.000701528253168926, 0.000828107993427705,
0.00079488114602085, 0.000904344998291552)), .Names = c("FragID",
"ParentID", "ParentSite.x", "FragSite.x", "TotalSA", "T03SA"),
class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L, 22L, 28L,
29L))
Я хочу, чтобы внутренний отсек -> Внутренний отсек и внешний отсек -> Внешний отсек были сплошными линиями, а внутренний отсек -> Внешний отсек и внешний отсек -> Внутренний отсек были пунктирными. То, что у меня есть до сих пор, приведено ниже:
GPBWM03 = ggplot(M3, aes(x = TotalSA, y = T03SA,
group = interaction(FragSite.x, ParentSite.x), color = FragSite.x,
linetype = ParentSite.x))
geom_point(alpha = 0.5, stroke = 0)
geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange =
TRUE)
scale_color_manual(values = c("Inner Bay" = "coral2",
"Outer Bay" = "skyblue4"))
Комментарии:
1. Не совсем уверен, что вы имеете в виду, но я думаю, что вы хотите просто добавить переменную, которая указывает, что либо это наблюдение находится внутри групп (от внутреннего к внутреннему или от внешнего к внешнему), либо между группами (от внутреннего к внешнему или от внешнего к внутреннему). Вы можете добавить это, а затем назначить ему тип строки
Ответ №1:
Другим способом было бы создать дополнительную переменную, которая кодирует linetype . В результате вы получаете легенду, которая показывает два фактора в отдельности:
M3 <- structure(list(FragID = c("1004", "1006", "1015", "1038", "1087", "1089", "1107", "1116"),
ParentID = c("166", "166", "166", "166", "144", "144", "144", "144"),
ParentSite.x = c("Inner Bay", "Inner Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay", "Outer Bay"),
FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"),
TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604, 0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258, 0.0356299659679559),
T03SA = c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665, 0.000966259734683879, 0.000701528253168926, 0.000828107993427705, 0.00079488114602085, 0.000904344998291552)),
.Names = c("FragID", "ParentID", "ParentSite.x", "FragSite.x", "TotalSA", "T03SA"),
class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L, 22L, 28L, 29L))
M3$in.out <- interaction(M3$ParentSite.x, M3$FragSite.x)
M3$comb <- 'same'
M3$comb[which(M3$ParentSite.x!=M3$FragSite.x)] <- 'different'
GPBWM03 = ggplot(M3, aes(x = TotalSA, y = T03SA,
group = interaction(FragSite.x, ParentSite.x), color = FragSite.x,
linetype = comb))
geom_point(alpha = 0.5, stroke = 0)
geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange =
TRUE)
scale_color_manual(values = c("Inner Bay" = "coral2",
"Outer Bay" = "skyblue4"))
GPBWM03
Ответ №2:
Это то, что вам нужно?
ggplot(M3,
aes(
x = TotalSA,
y = T03SA,
color = FragSite.x,
linetype = interaction(ParentSite.x, FragSite.x)
))
geom_smooth(
method = "lm",
formula = y ~ x,
se = FALSE,
fullrange =
TRUE
)
scale_linetype_manual(values = c(1, 2, 2, 1))
scale_color_manual(values =
c("Inner Bay" = "coral2", "Outer Bay" = "skyblue4"))
Чтобы получить вывод, подобный:
Это было сделано с использованием приведенного ниже набора данных, поскольку в некоторых значениях в вашем исходном dput выводе выше были разрывы строк:
structure(list(FragID = c("1004", "1006", "1015", "1038", "1087",
"1089", "1107", "1116"), ParentID = c("166", "166", "166", "166",
"144", "144", "144", "144"), ParentSite.x = c("Inner Bay", "Inner Bay",
"Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay",
"Outer Bay"), FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay",
"Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"
), TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604,
0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258,
0.0356299659679559), T03SA = c(0.000709604810935872, 0.00148788795124323,
0.00109901406229665, 0.000966259734683879, 0.000701528253168926,
0.000828107993427705, 0.00079488114602085, 0.000904344998291552
)), class = "data.frame", row.names = c(2L, 3L, 5L, 9L, 21L,
22L, 28L, 29L))
Редактировать: для комбинированной легенды может быть проще просто вручную назначить цвета в дополнение к типам линий:
ggplot(M3,
aes(
x = TotalSA,
y = T03SA,
color = interaction(ParentSite.x, FragSite.x),
linetype = interaction(ParentSite.x, FragSite.x)
))
geom_smooth(
method = "lm",
formula = y ~ x,
se = FALSE,
fullrange =
TRUE
)
scale_linetype_manual(name = "Combination", values = c(1, 2, 2, 1))
scale_color_manual(name = "Combination", values = c("coral2", "coral2", "skyblue4", "skyblue4"))
Вывод:
В качестве альтернативы, вам может быть проще создать переменную взаимодействия, как предложено некоторыми другими ответами здесь.
Комментарии:
1. Как мне отформатировать легенду, чтобы это была одна легенда, отображающая как цвет, так и тип линии? Таким образом, объединяя 2 легенды в одну.
Ответ №3:
У вас dput был неправильный формат. Я создал его вручную, используя следующее:
df <- data.frame(FragID = c("1004", "1006", "1015", "1038", "1087", "1089", "1107", "1116"),
ParentID = c("166", "166", "166", "166", "144", "144", "144", "144"),
ParentSite.x = c("Inner Bay", "Inner Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay", "Outer Bay", "Outer Bay"),
FragSite.x = c("Inner Bay", "Outer Bay", "Inner Bay", "Outer Bay", "Inner Bay", "Inner Bay", "Outer Bay", "Outer Bay"),
TotalSA = c(0.0171970755726674, 0.0338197513082082, 0.0215722601402604, 0.030712272182997, 0.027529366126288, 0.029650482611575, 0.0316984120058258, 0.0356299659679559),
T03SA = c(0.000709604810935872, 0.00148788795124323, 0.00109901406229665, 0.000966259734683879, 0.000701528253168926, 0.000828107993427705, 0.00079488114602085, 0.000904344998291552))
Затем я добавил новый столбец к вашим данным, который определяет тип строки на основе взаимодействия ParentSite.x и FragSite.x использования tidyverse пакета:
df <- df %>%
mutate(linetype = ifelse(ParentSite.x=='Inner Bay'amp;FragSite.x=='Inner Bay' | ParentSite.x=='Outer Bay'amp;FragSite.x=='Outer Bay'
, 'Matching', 'Not Matching'))
Затем мы строим:
ggplot(df, aes(x = TotalSA, y = T03SA, color = FragSite.x, linetype = linetype))
geom_point(alpha = 0.5, stroke = 0)
geom_smooth(method = "lm", formula = y~x, se = FALSE, fullrange = TRUE)
scale_color_manual(values = c("Inner Bay" = "coral2",
"Outer Bay" = "skyblue4"))
labs(linetype='Line Type')
