Создание гистограммы нескольких лет в R
У меня есть набор данных с данными за несколько лет. Я взял каждый год и получил среднее значение, которое показано ниже.
## structure(list(X = 1:10, Sector = c("B01", "B02", "B03", "B04",
"B05", "B06", "B07", "B08", "B09", "B10"), Octant = c("NW", "N",
"NE", "NW", "NW", "N", "N", "NE", "NE", "W"), Distance_m = c(1600L,
1600L, 1600L, 3200L, 3200L, 3200L, 3200L, 3200L, 3200L, 4800L
), Yr1951 = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), Yr1952 = c(4.16722556361885,
8.07009038922937, 19.9170847727104, 3.14049476463692, 3.14049476463692,
6.08176260934446, 6.08176260934446, 15.0098667568185, 15.0098667568185,
4.65320683317187), Yr1953 = c(15.7901552659006, 30.578613590843,
75.4684036666576, 11.8997398120925, 11.8997398120925, 23.0445831227103,
23.0445831227103, 56.8743215341634, 56.8743215341634, 17.6316009916985
), Yr1954 = c(23.2458594045635, 45.0170464032273, 111.102637787894,
17.5184900949112, 17.5184900949112, 33.9256410267316, 33.9256410267316,
83.7289095546898, 83.7289095546898, 25.9567883170532), Yr1955 = c(23.2458594045635,
45.0170464032273, 111.102637787894, 17.5184900949112, 17.5184900949112,
33.9256410267316, 33.9256410267316, 83.7289095546898, 83.7289095546898,
25.9567883170532), Yr1956 = c(23.2458594045635, 45.0170464032273,
111.102637787894, 17.5184900949112, 17.5184900949112, 33.9256410267316,
33.9256410267316, 83.7289095546898, 83.7289095546898, 25.9567883170532
), Yr1957 = c(23.2458594045635, 45.0170464032273, 111.102637787894,
17.5184900949112, 17.5184900949112, 33.9256410267316, 33.9256410267316,
83.7289095546898, 83.7289095546898, 25.9567883170532), Yr1958 = c(24.6439039639345,
47.724446276369, 117.78453479105, 18.5720811598549, 18.5720811598549,
35.9659853751658, 35.9659853751658, 88.7645051172278, 88.7645051172278,
27.517872639807), Yr1959 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1960 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1961 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1962 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1963 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1964 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1965 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1966 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1967 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1968 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1969 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1970 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1971 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1972 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1973 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1974 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1975 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1976 = c(28.8380376420477,
55.8466458957943, 137.830225800518, 21.7328543546861, 21.7328543546861,
42.0870184204684, 42.0870184204684, 103.871291804842, 103.871291804842,
32.2011256080683), Yr1977 = c(28.8380376420477, 55.8466458957943,
137.830225800518, 21.7328543546861, 21.7328543546861, 42.0870184204684,
42.0870184204684, 103.871291804842, 103.871291804842, 32.2011256080683
), Yr1978 = c(28.8380376420477, 55.8466458957943, 137.830225800518,
21.7328543546861, 21.7328543546861, 42.0870184204684, 42.0870184204684,
103.871291804842, 103.871291804842, 32.2011256080683), Yr1979 = c(16.6426332897682,
32.2294900798683, 79.5427876444474, 12.5421823029183, 12.5421823029183,
24.2887127940458, 24.2887127940458, 59.9448492404336, 59.9448492404336,
18.5834948849451), Yr1980 = c(4.44722893748883, 8.61233426394242,
21.2553494883767, 3.35151025115061, 3.35151025115061, 6.49040716762317,
6.49040716762317, 16.0184066760251, 16.0184066760251, 4.96586416182185
), Yr1981 = c(4.44722893748883, 8.61233426394242, 21.2553494883767,
3.35151025115061, 3.35151025115061, 6.49040716762317, 6.49040716762317,
16.0184066760251, 16.0184066760251, 4.96586416182185), Yr1982 = c(4.44722893748883,
8.61233426394242, 21.2553494883767, 3.35151025115061, 3.35151025115061,
6.49040716762317, 6.49040716762317, 16.0184066760251, 16.0184066760251,
4.96586416182185), Yr1983 = c(4.44722893748883, 8.61233426394242,
21.2553494883767, 3.35151025115061, 3.35151025115061, 6.49040716762317,
6.49040716762317, 16.0184066760251, 16.0184066760251, 4.96586416182185
), Yr1984 = c(4.44722893748883, 8.61233426394242, 21.2553494883767,
3.35151025115061, 3.35151025115061, 6.49040716762317, 6.49040716762317,
16.0184066760251, 16.0184066760251, 4.96586416182185), Yr1985 = c(4.44722893748883,
8.61233426394242, 21.2553494883767, 3.35151025115061, 3.35151025115061,
6.49040716762317, 6.49040716762317, 16.0184066760251, 16.0184066760251,
4.96586416182185), Yr1986 = c(4.44722893748883, 8.61233426394242,
21.2553494883767, 3.35151025115061, 3.35151025115061, 6.49040716762317,
6.49040716762317, 16.0184066760251, 16.0184066760251, 4.96586416182185
), Yr1987 = c(4.44722893748883, 8.61233426394242, 21.2553494883767,
3.35151025115061, 3.35151025115061, 6.49040716762317, 6.49040716762317,
16.0184066760251, 16.0184066760251, 4.96586416182185), Yr1988 = c(1.1104088237638,
2.15037545723583, 5.30715372557456, 0.836823695861693, 0.836823695861693,
1.62056091333541, 1.62056091333541, 3.99956475497726, 3.99956475497726,
1.23990454739508)), class = "data.frame", row.names = c(NA, -10L
))
Мне нужно создать гистограмму, содержащую отдельные столбцы для каждого года. Как бы я это сделал? Ось y можно назвать экспозицией, а x — годом. У меня не было опыта, кроме очень простых гистограмм. Мне также было интересно, как я могу объединить все эти годы вместе и получить общее среднее значение и медиану.
@jdobres я отредактировал, чтобы добавить исходный код. Я не был уверен, как получить средства каждого года для создания гистограммы, поэтому изначально я их разделил.
—
19.03.2022 16:15
Стоит ли изучать PHP в 2023-2024 годах?
Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.
Приемы CSS-макетирования - floats и Flexbox
Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...
Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest
В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...
Концепция локализации и ее применение в приложениях React ⚡️
Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...
Пользовательский скаляр GraphQL
Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...
1
2
39
1
Перейти к ответу
Данный вопрос помечен как решенный
Ответы 1
Ответ принят как подходящий
Для таких операций очень удобен набор пакетов tidyverse. tidyverse также включает ggplot, отличную библиотеку для создания визуализаций. Существует множество отличных ресурсов для получения дополнительной информации о tidyverse, таких как Вот этот.
Здесь я назвал ваши данные df. Во-первых, давайте преобразуем данные в «длинный» формат, в котором каждая строка содержит одно интересующее значение данных, а затем рассчитаем среднее значение за год:
library(tidyverse)
df_means <- df %>%
pivot_longer(starts_with('Yr'), names_to = 'year') %>%
mutate(year_numeric = as.numeric(gsub('Yr', '', year))) %>%
group_by(year_numeric) %>%
summarize(value = mean(value, na.rm = T))
year_numeric value
<dbl> <dbl>
1 1951 0
2 1952 8.53
3 1953 32.3
4 1954 47.6
5 1955 47.6
6 1956 47.6
7 1957 47.6
8 1958 50.4
9 1959 59.0
10 1960 59.0
# … with 28 more rows
Мы можем выполнить аналогичную операцию без команды group_by, чтобы получить среднее значение и медиану:
df_overall <- df %>%
pivot_longer(starts_with('Yr'), names_to = 'year') %>%
summarize(
grand_mean = mean(value, na.rm = T),
grand_median = median(value, na.rm = T)
) %>%
pivot_longer(everything(), names_to = 'measure')
measure value
<chr> <dbl>
1 grand_mean 41.3
2 grand_median 28.8
И, наконец, мы можем создать график с ggplot. Обратите внимание, как мало графического кода требуется для создания чего-то относительно сложного:
plot_means <- df_means %>%
ggplot(data = ., aes(x = year_numeric, y = value)) +
geom_col() +
geom_hline(data = df_overall, aes(color = measure, yintercept = value)) +
labs(x = 'Year', y = 'Exposure')
print(plot_means)
Благодарю вас! Это работает хорошо, за исключением того, что в моем полном наборе данных есть NA, поэтому фактическая гистограмма не может быть сформирована. Где в коде я могу вставить na.rm=TRUE, чтобы получить желаемый результат?
—
19.03.2022 17:30
Внутри вызовов функций mean и median (см. редактирование). Если мой ответ работает для вас, отметьте его как правильный и, при желании, проголосуйте за него.
—
19.03.2022 17:30
Другие вопросы по теме
Вторичная ось в процентах ggplot2
Проблема с шириной столбцов в морской гистограмме
Оси x в визуализации графика в matplotlib
Как центрировать столбцы в R barplot() вокруг тиков
Могу ли я сделать горизонтальную гистограмму с двумя осями категорий?
Как установить цвет полосы в barChart в JasperStudio на основе значения в наборе данных?
Гистограмма с накоплением (с процентами) только одной переменной в R?
Гистограмма из Chart.js, не отображающая второй набор данных
Построение гистограмм рядом друг с другом
Есть ли способ изменить плитки элементов легенды C3 на значки?
Похожие вопросы
Как получить эксклюзивные диапазоны дат в R?
Правильно размещайте метки geom_text на лесном участке с параметром «fill» в ggplot2
Работа со списками в R — объединение двух фреймов данных и вычисление значений с несовпадающими размерами
Почему strsplit не работает с моей строкой символов?
Чтение CSV в R создает новые столбцы с X
R: Проблема с вводом `mutate()` `index_date`. Строка символов x не имеет стандартного однозначного формата. i Ввод `index_date` равен `as.Date(...)
Рисование трехосевого графика в ggplot2 с помощью r
Как изменить цвета в гистограммах с несколькими стопками в R?
Команда R для извлечения текста между двумя строками, содержащими фигурные скобки
График в цикле
Похоже, что среднее значение за каждый год хранится как отдельная переменная в вашей рабочей области. Это не очень хороший способ управления данными, и это очень затрудняет создание графиков. Можете ли вы поделиться образцом исходных данных, на которых это основано?