#Source: https://t.me/rstudioprogr
library(AER)
## Warning: пакет 'AER' был собран под R версии 4.3.3
## Загрузка требуемого пакета: car
## Загрузка требуемого пакета: carData
## Загрузка требуемого пакета: lmtest
## Загрузка требуемого пакета: zoo
##
## Присоединяю пакет: 'zoo'
## Следующие объекты скрыты от 'package:base':
##
## as.Date, as.Date.numeric
## Загрузка требуемого пакета: sandwich
## Загрузка требуемого пакета: survival
library(ggstatsplot)
## You can cite this package as:
## Patil, I. (2021). Visualizations with statistical details: The 'ggstatsplot' approach.
## Journal of Open Source Software, 6(61), 3167, doi:10.21105/joss.03167
library(gridExtra)
library(dplyr)
##
## Присоединяю пакет: 'dplyr'
## Следующий объект скрыт от 'package:gridExtra':
##
## combine
## Следующий объект скрыт от 'package:car':
##
## recode
## Следующие объекты скрыты от 'package:stats':
##
## filter, lag
## Следующие объекты скрыты от 'package:base':
##
## intersect, setdiff, setequal, union
data("Affairs")
head(Affairs)
## affairs gender age yearsmarried children religiousness education occupation
## 4 0 male 37 10.00 no 3 18 7
## 5 0 female 27 4.00 no 4 14 6
## 11 0 female 32 15.00 yes 1 12 1
## 16 0 male 57 15.00 yes 5 18 6
## 23 0 male 22 0.75 no 2 17 6
## 29 0 female 32 1.50 no 2 17 5
## rating
## 4 4
## 5 4
## 11 4
## 16 5
## 23 3
## 29 5
Affairs1 <- Affairs
Affairs1$children <- ifelse(Affairs$children == "no", 0, 1)
g1 <- ggcorrmat(data = Affairs1 %>% dplyr::filter(gender == "male") %>% dplyr::select(-gender),
colors = c("#187ABB", "white", "#3EB561"), title = "Корреляция внутри подвыборки мужчин")
g2 <- ggcorrmat(data = Affairs1 %>% dplyr::filter(gender == "female") %>% dplyr::select(-gender),
colors = c("#187ABB", "white", "#3EB561"), title = "Корреляция внутри подвыборки женщин")
grid.arrange(g1, g2, ncol = 2)
