summary(cars)
## speed dist
## Min. : 4.0 Min. : 2.00
## 1st Qu.:12.0 1st Qu.: 26.00
## Median :15.0 Median : 36.00
## Mean :15.4 Mean : 42.98
## 3rd Qu.:19.0 3rd Qu.: 56.00
## Max. :25.0 Max. :120.00
file_path <- "~/Downloads/База_КлимРиск_2023.sav"
data <- read_sav(file_path)
# Извлекаю переменную для анализа
variable <- iris$Sepal.Length
# Провожу тест Шапиро-Уилка
shapiro_test <- shapiro.test(variable)
shapiro_test
##
## Shapiro-Wilk normality test
##
## data: variable
## W = 0.97609, p-value = 0.01018
# Провожу тест Колмогорова-Смирнова
ks_test <- ks.test(variable, "pnorm", mean = mean(variable), sd = sd(variable))
## Warning in ks.test.default(variable, "pnorm", mean = mean(variable), sd =
## sd(variable)): ties should not be present for the one-sample Kolmogorov-Smirnov
## test
ks_test
##
## Asymptotic one-sample Kolmogorov-Smirnov test
##
## data: variable
## D = 0.088654, p-value = 0.1891
## alternative hypothesis: two-sided
# Провожу тест Андерсона-Дарлинга
ad_test <- ad.test(variable)
ad_test
##
## Anderson-Darling normality test
##
## data: variable
## A = 0.8892, p-value = 0.02251
# Одномерный анализ
# Таблицы частот
table(data$V12)
##
## 1 2 3 88 99
## 190 437 189 19 64
table(data$V13)
##
## 1 2 3
## 575 149 166
table(data$V15)
##
## 1 2 3
## 586 153 154
# Столбчатые диаграммы
ggplot(data, aes(x = V12)) +
geom_bar(color = "red") +
ggtitle("Распределение V12") +
xlab("V12") +
ylab("Частота") +
theme_minimal()
## Warning: Removed 14 rows containing non-finite outside the scale range
## (`stat_count()`).
ggplot(data, aes(x = V13)) +
geom_bar(color = "red") +
ggtitle("Распределение V13") +
xlab("V13") +
ylab("Частота") +
theme_minimal()
## Warning: Removed 23 rows containing non-finite outside the scale range
## (`stat_count()`).
ggplot(data, aes(x = V15)) +
geom_bar(color = "red") +
ggtitle("Распределение V15") +
xlab("V15") +
ylab("Частота") +
theme_minimal()
## Warning: Removed 20 rows containing non-finite outside the scale range
## (`stat_count()`).
# Таблицы сопряженности
table(data$Region, data$V12)
##
## 1 2 3 88 99
## 1 67 91 74 1 32
## 2 91 252 73 13 15
## 3 32 94 42 5 17
## 4 0 0 0 0 0
table(data$Region, data$V13)
##
## 1 2 3
## 1 106 84 76
## 2 327 47 61
## 3 142 18 29
## 4 0 0 0
table(data$Region, data$V15)
##
## 1 2 3
## 1 133 57 76
## 2 310 76 51
## 3 143 20 27
## 4 0 0 0
# Создание факторной переменной для регионов с метками
data$Region <- factor(data$Region, levels = c(1, 2, 3, 4),
labels = c("Алтайский край", "Республика Алтай", "Республика Тыва", "Монголия"))
# Столбчатая диаграмма для V12 по регионам
ggplot(data, aes(x = Region, fill = as.factor(V12))) +
geom_bar(position = "dodge") +
ggtitle("Распределение V12 по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "V12") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V13 по регионам
ggplot(data, aes(x = Region, fill = as.factor(V13))) +
geom_bar(position = "dodge") +
ggtitle("Распределение V13 по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "V13") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V15 по регионам
ggplot(data, aes(x = Region, fill = as.factor(V15))) +
geom_bar(position = "dodge") +
ggtitle("Распределение V15 по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "V15") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Преобразование переменных V14_1, V14_2 и т.д. в факторы с метками
data$V14_1 <- factor(data$V14_1, levels = c(0, 1), labels = c("нет", "да"))
data$V14_2 <- factor(data$V14_2, levels = c(0, 1), labels = c("нет", "да"))
data$V14_3 <- factor(data$V14_3, levels = c(0, 1), labels = c("нет", "да"))
data$V14_4 <- factor(data$V14_4, levels = c(0, 1), labels = c("нет", "да"))
data$V14_5 <- factor(data$V14_5, levels = c(0, 1), labels = c("нет", "да"))
data$V14_6 <- factor(data$V14_6, levels = c(0, 1), labels = c("нет", "да"))
data$V14_7 <- factor(data$V14_7, levels = c(0, 1), labels = c("нет", "да"))
data$V14_8 <- factor(data$V14_8, levels = c(0, 1), labels = c("нет", "да"))
data$V14_9 <- factor(data$V14_9, levels = c(0, 1), labels = c("нет", "да"))
data$V14_10 <- factor(data$V14_10, levels = c(0, 1), labels = c("нет", "да"))
data$V14_88 <- factor(data$V14_88, levels = c(0, 1), labels = c("нет", "да"))
# Столбчатая диаграмма для V14_1 по регионам
ggplot(data, aes(x = Region, fill = V14_1)) +
geom_bar(position = "dodge") +
ggtitle("Распределение резких перепадов температуры по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Резкие перепады температуры") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_2 по регионам
ggplot(data, aes(x = Region, fill = V14_2)) +
geom_bar(position = "dodge") +
ggtitle("Распределение затяжных морозов по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Затяжные морозы") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_3 по регионам
ggplot(data, aes(x = Region, fill = V14_3)) +
geom_bar(position = "dodge") +
ggtitle("Распределение частых оттепелей по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Частые оттепели") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_4 по регионам
ggplot(data, aes(x = Region, fill = V14_4)) +
geom_bar(position = "dodge") +
ggtitle("Распределение увеличения гололеда на дорогах по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Увеличение гололеда на дорогах") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_5 по регионам
ggplot(data, aes(x = Region, fill = V14_5)) +
geom_bar(position = "dodge") +
ggtitle("Распределение обильных снегопадов по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Обильные снегопады") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_6 по регионам
ggplot(data, aes(x = Region, fill = V14_6)) +
geom_bar(position = "dodge") +
ggtitle("Распределение схода снежных лавин по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Сход снежных лавин") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_7 по регионам
ggplot(data, aes(x = Region, fill = V14_7)) +
geom_bar(position = "dodge") +
ggtitle("Распределение увеличения количества пасмурных дней по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Увеличение количества пасмурных дней") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_8 по регионам
ggplot(data, aes(x = Region, fill = V14_8)) +
geom_bar(position = "dodge") +
ggtitle("Распределение раннего таяния и позднего образования речного льда по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Раннее таяние и позднее образование речного льда") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_9 по регионам
ggplot(data, aes(x = Region, fill = V14_9)) +
geom_bar(position = "dodge") +
ggtitle("Распределение участившихся ледовых зажоров по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Участившиеся ледовые зажоры") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_10 по регионам
ggplot(data, aes(x = Region, fill = V14_10)) +
geom_bar(position = "dodge") +
ggtitle("Распределение усиления ветров и метелей по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Усиление ветров и метелей") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V14_88 по регионам
ggplot(data, aes(x = Region, fill = V14_88)) +
geom_bar(position = "dodge") +
ggtitle("Распределение других факторов по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Другие факторы") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Преобразование переменных V16_1, V16_2 и т.д. в факторы с метками
data$V16_1 <- factor(data$V16_1, levels = c(0, 1), labels = c("нет", "да"))
data$V16_2 <- factor(data$V16_2, levels = c(0, 1), labels = c("нет", "да"))
data$V16_3 <- factor(data$V16_3, levels = c(0, 1), labels = c("нет", "да"))
data$V16_4 <- factor(data$V16_4, levels = c(0, 1), labels = c("нет", "да"))
data$V16_5 <- factor(data$V16_5, levels = c(0, 1), labels = c("нет", "да"))
data$V16_6 <- factor(data$V16_6, levels = c(0, 1), labels = c("нет", "да"))
data$V16_7 <- factor(data$V16_7, levels = c(0, 1), labels = c("нет", "да"))
data$V16_8 <- factor(data$V16_8, levels = c(0, 1), labels = c("нет", "да"))
data$V16_9 <- factor(data$V16_9, levels = c(0, 1), labels = c("нет", "да"))
data$V16_88 <- factor(data$V16_88, levels = c(0, 1), labels = c("нет", "да"))
# Столбчатая диаграмма для V16_1 по регионам
ggplot(data, aes(x = Region, fill = V16_1)) +
geom_bar(position = "dodge") +
ggtitle("Распределение увеличения количества засушливых дней по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Увеличение количества засушливых дней") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_2 по регионам
ggplot(data, aes(x = Region, fill = V16_2)) +
geom_bar(position = "dodge") +
ggtitle("Распределение увеличения периодов аномальной жары по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Увеличение периодов аномальной жары") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_3 по регионам
ggplot(data, aes(x = Region, fill = V16_3)) +
geom_bar(position = "dodge") +
ggtitle("Распределение большого количества осадков по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Большое количество осадков") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_4 по регионам
ggplot(data, aes(x = Region, fill = V16_4)) +
geom_bar(position = "dodge") +
ggtitle("Распределение увеличения количества камнепадов и оползней по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Увеличение количества камнепадов и оползней") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_5 по регионам
ggplot(data, aes(x = Region, fill = V16_5)) +
geom_bar(position = "dodge") +
ggtitle("Распределение усиления паводков по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Усиление паводков") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_6 по регионам
ggplot(data, aes(x = Region, fill = V16_6)) +
geom_bar(position = "dodge") +
ggtitle("Распределение таяния вечной мерзлоты по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Таяние вечной мерзлоты") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_7 по регионам
ggplot(data, aes(x = Region, fill = V16_7)) +
geom_bar(position = "dodge") +
ggtitle("Распределение сильных ветров и штормов по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Сильные ветры и штормы") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_8 по регионам
ggplot(data, aes(x = Region, fill = V16_8)) +
geom_bar(position = "dodge") +
ggtitle("Распределение роста количества насекомых по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Рост количества насекомых") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_9 по регионам
ggplot(data, aes(x = Region, fill = V16_9)) +
geom_bar(position = "dodge") +
ggtitle("Распределение роста количества насекомых-вредителей по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Рост количества насекомых-вредителей") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
# Столбчатая диаграмма для V16_88 по регионам
ggplot(data, aes(x = Region, fill = V16_88)) +
geom_bar(position = "dodge") +
ggtitle("Распределение других факторов по регионам") +
xlab("Регион") +
ylab("Частота") +
scale_fill_discrete(name = "Другие факторы") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1))
