ЗАДАНИЕ 1. anova параметрическая

library(datarium )
data("selfesteem")
library(dplyr)
## 
## Присоединяю пакет: 'dplyr'
## Следующие объекты скрыты от 'package:stats':
## 
##     filter, lag
## Следующие объекты скрыты от 'package:base':
## 
##     intersect, setdiff, setequal, union
library(tidyr)

library(rstatix)
## 
## Присоединяю пакет: 'rstatix'
## Следующий объект скрыт от 'package:stats':
## 
##     filter
dieta<-selfesteem %>% 
  pivot_longer(cols=contains("t"), values_to = "time")

library(rstatix)
dieta %>%
  group_by(name) %>%
  get_summary_stats(time, type = "mean_sd")

у людей повышается самооценка, хотя и в третьем периоде оценки будут разниться

library(ggplot2)
library(ggpubr)
bxp <- ggboxplot(dieta, x = "name", y = "time", fill="time", add = "point",  )
bxp
## Warning: The following aesthetics were dropped during statistical transformation: fill.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?

график подтверждает

dieta %>%
  group_by(name) %>%
  shapiro_test(time)

нормальное распределение

ggqqplot(dieta, "time", facet.by = "name")

res.aov <- anova_test(data = dieta, dv = time, wid=id, within = name)
get_anova_table(res.aov)

самооценки были статистически различны , и эффект большой

pwc <- dieta %>%
  pairwise_t_test(
    time ~ name, paired = TRUE,
    p.adjust.method = "bonferroni"
  )
pwc

различия в самооценках между каждыми временными периодами значимы, особенно в начале и в конце

pwc <- pwc %>% add_xy_position(x = "name")
bxp + 
  stat_pvalue_manual(pwc) +
  labs(
    subtitle = get_test_label(res.aov, detailed = TRUE),
    caption = get_pwc_label(pwc)
  )
## Warning: The following aesthetics were dropped during statistical transformation: fill.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?

ЗАДАНИЕ 1. anova непараметрический аналог дисперсионного анализа для зависимых выборок

res.fried <- dieta %>% friedman_test(time ~ name |id)
res.fried

разлтчия значимы

dieta %>% friedman_effsize(time ~ name |id)

взаимосвязь очень тесная

pwc <- dieta %>%
  wilcox_test(time ~ name, paired = TRUE, p.adjust.method = "bonferroni")
pwc

здесь тоже показывает значимую разницу между временными периодами, только 1 и 2 замеры такие же по значимости, как 1 и 3

pwc <- pwc %>% add_xy_position(x = "name")
ggboxplot(dieta, x = "name", y = "time", fill="name", add = "point") +
  stat_pvalue_manual(pwc) +
  labs(
    subtitle = get_test_label(res.fried,  detailed = TRUE),
    caption = get_pwc_label(pwc)
  )




ЗАДАНИЕ 2.

library(haven)
kl<-read_sav("C:/Users/malya/Desktop/База_КлимРиск_2023.sav")
chisq.test(kl$V1, kl$V38)
## Warning in chisq.test(kl$V1, kl$V38): аппроксимация на основе хи-квадрат может
## быть неправильной
## 
##  Pearson's Chi-squared test
## 
## data:  kl$V1 and kl$V38
## X-squared = 15.18, df = 6, p-value = 0.0189

на “Обработке” я брала эти же переменные - получились аналогичные значения.
графикик - мозаичная диаграмма

library(vcd)
## Загрузка требуемого пакета: grid
kl$V1<-sjlabelled::as_label(kl$V1)
mosaic( ~ V38 + V1, data = kl, ylab="",ylab="" )

ассоциативный график - ??? я не нашла в интернете синтаксис, а с тем, который в лекции, не получается построить