rm(list = ls())
date()
## [1] "Sat Jan 18 13:14:10 2025"
sessionInfo()
## R version 4.4.2 (2024-10-31)
## Platform: x86_64-apple-darwin20
## Running under: macOS Ventura 13.7.2
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/lib/libRlapack.dylib; LAPACK version 3.12.0
##
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
##
## time zone: Europe/Moscow
## tzcode source: internal
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## loaded via a namespace (and not attached):
## [1] digest_0.6.36 R6_2.5.1 fastmap_1.2.0 xfun_0.46
## [5] cachem_1.1.0 knitr_1.48 htmltools_0.5.8.1 rmarkdown_2.27
## [9] lifecycle_1.0.4 cli_3.6.3 sass_0.4.9 jquerylib_0.1.4
## [13] compiler_4.4.2 rstudioapi_0.16.0 tools_4.4.2 evaluate_0.24.0
## [17] bslib_0.7.0 yaml_2.3.9 rlang_1.1.4 jsonlite_1.8.8
options(scipen = 999) # Убирает научную запись чисел
library(ggplot2)
library(tidyr)
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(psych)
##
## Attaching package: 'psych'
## The following objects are masked from 'package:ggplot2':
##
## %+%, alpha
load(file = "Data.RData")
Data2 <- select(Data, !(Age:ResidenceAdult))
По стране
Data %>%
#select(1:13) %>% #Берем только первые 13 колонок
pivot_longer(cols = starts_with("W_"), names_to = "Variable", values_to = "Value") %>% # Колонки, которые начинаются с "W_"
ggplot(., aes(x = Value, y = EnvIdentity, col = Nationality)) +
geom_point() +
geom_smooth(method = 'lm', se = TRUE) + # se - убирает стандартную ошибку st. error
facet_grid(cols = vars(Variable)) +
labs(title = "Well-being", x = "Well-Being", y = "EnvIdentity") +
theme_light()
## `geom_smooth()` using formula = 'y ~ x'
Благополучие по полу
Data %>%
#select(1:13) %>% #Берем только первые 13 колонок
pivot_longer(cols = starts_with("W_"), names_to = "Variable", values_to = "Value") %>% # Колонки, которые начинаются с "W_"
ggplot(., aes(x = Value, y = EnvIdentity, col = Sex)) +
geom_point() +
geom_smooth(method = 'lm', se = FALSE) + # se - убирает стандартную ошибку st. error se = FALSE
facet_grid(cols = vars(Variable)) +
labs(title = "Well-being", x = "Well-Being", y = "EnvIdentity") +
theme_light()
## `geom_smooth()` using formula = 'y ~ x'
По стране
Data %>%
#select(1:13) %>% #Берем только первые 13 колонок
pivot_longer(cols = starts_with("V_"), names_to = "Variable", values_to = "Value") %>% # Колонки, которые начинаются с "V_"
ggplot(., aes(x = Value, y = EnvIdentity, col = Nationality)) +
geom_point() +
geom_smooth(method = 'lm', se = FALSE) + # se - убирает стандартную ошибку st. error
facet_grid(cols = vars(Variable)) +
labs(title = "Values", x = "Values", y = "EnvIdentity") +
theme_light()
## `geom_smooth()` using formula = 'y ~ x'
По полу
Data %>%
#select(1:13) %>% #Берем только первые 13 колонок
pivot_longer(cols = starts_with("V_"), names_to = "Variable", values_to = "Value") %>% # Колонки, которые начинаются с "V_"
ggplot(., aes(x = Value, y = EnvIdentity, col = Sex)) +
geom_point() +
geom_smooth(method = 'lm', se = FALSE) + # se - убирает стандартную ошибку st. error
facet_grid(cols = vars(Variable)) +
labs(title = "Values", x = "Values", y = "EnvIdentity") +
theme_light()
## `geom_smooth()` using formula = 'y ~ x'
corPlot(Data[,6:23], stars = TRUE, MAR = 1, xlas = 2)
EnvIdentity ~ Nationality + Sex + W_PositiveRel + W_Autonomy + W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance + WellBeing + V_Conformity + V_tradition + V_benevolence + V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism + V_Achievement + V_Power + V_Security
##Все факторы
lm(EnvIdentity ~ ., Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ Nationality + W_Autonomy + V_benevolence +
## V_universalism + V_Selfdirection, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -27.4914 -5.9294 0.2515 5.8118 22.1345
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 46.4738 9.6637 4.809 0.00000519 ***
## NationalityRussia 8.3964 2.8441 2.952 0.00391 **
## W_Autonomy -0.3478 0.1574 -2.209 0.02935 *
## V_benevolence 2.6353 1.3870 1.900 0.06022 .
## V_universalism 5.3193 1.2979 4.098 0.00008313 ***
## V_Selfdirection -2.2027 1.4232 -1.548 0.12476
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11 on 103 degrees of freedom
## Multiple R-squared: 0.2771, Adjusted R-squared: 0.242
## F-statistic: 7.896 on 5 and 103 DF, p-value: 0.000002502
##Все без пола
lm(EnvIdentity ~ Nationality + W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ Nationality + W_Autonomy + V_benevolence +
## V_universalism + V_Selfdirection, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -27.4914 -5.9294 0.2515 5.8118 22.1345
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 46.4738 9.6637 4.809 0.00000519 ***
## NationalityRussia 8.3964 2.8441 2.952 0.00391 **
## W_Autonomy -0.3478 0.1574 -2.209 0.02935 *
## V_benevolence 2.6353 1.3870 1.900 0.06022 .
## V_universalism 5.3193 1.2979 4.098 0.00008313 ***
## V_Selfdirection -2.2027 1.4232 -1.548 0.12476
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11 on 103 degrees of freedom
## Multiple R-squared: 0.2771, Adjusted R-squared: 0.242
## F-statistic: 7.896 on 5 and 103 DF, p-value: 0.000002502
##Все без страны
lm(EnvIdentity ~ Sex + W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ W_PositiveRel + W_Autonomy + W_PersGrowth +
## V_benevolence + V_universalism + V_Selfdirection + V_Power +
## V_Security, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -26.279 -5.403 -1.255 6.639 22.688
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 34.7222 10.9217 3.179 0.00197 **
## W_PositiveRel -0.2382 0.1634 -1.458 0.14797
## W_Autonomy -0.2588 0.1582 -1.636 0.10502
## W_PersGrowth 0.4556 0.1529 2.980 0.00361 **
## V_benevolence 2.2434 1.3979 1.605 0.11167
## V_universalism 4.6865 1.3396 3.499 0.00070 ***
## V_Selfdirection -2.2173 1.4855 -1.493 0.13869
## V_Power -1.0600 0.7710 -1.375 0.17221
## V_Security 1.5156 1.0659 1.422 0.15818
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.02 on 100 degrees of freedom
## Multiple R-squared: 0.2945, Adjusted R-squared: 0.238
## F-statistic: 5.217 on 8 and 100 DF, p-value: 0.00001843
lm(EnvIdentity ~ Nationality * Sex * (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security)
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ Nationality + Sex + W_PositiveRel +
## W_Autonomy + W_EnvMastery + W_PersGrowth + W_PurposeLife +
## V_Conformity + V_benevolence + V_universalism + V_Selfdirection +
## V_Hedonism + V_Power + V_Security + Nationality:Sex + Nationality:W_PositiveRel +
## Nationality:W_Autonomy + Nationality:W_PurposeLife + Nationality:V_Conformity +
## Nationality:V_benevolence + Nationality:V_Selfdirection +
## Nationality:V_Hedonism + Sex:W_PositiveRel + Sex:W_Autonomy +
## Sex:W_PurposeLife + Sex:V_Selfdirection + Sex:V_Hedonism +
## Sex:V_Security + Nationality:Sex:W_PositiveRel + Nationality:Sex:V_Selfdirection +
## Nationality:Sex:V_Hedonism, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -29.4577 -3.6079 0.1787 5.4925 19.1874
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 9.8743 31.5307 0.313 0.755003
## NationalityRussia 11.2410 36.1360 0.311 0.756584
## SexMale 94.4939 44.2530 2.135 0.035915
## W_PositiveRel 0.3922 0.4054 0.967 0.336327
## W_Autonomy 0.4216 0.3504 1.203 0.232562
## W_EnvMastery -0.4751 0.2439 -1.948 0.055036
## W_PersGrowth 0.3816 0.2622 1.455 0.149638
## W_PurposeLife -0.2901 0.2283 -1.271 0.207721
## V_Conformity 1.8037 2.1641 0.833 0.407154
## V_benevolence 4.6220 2.1866 2.114 0.037773
## V_universalism 5.1132 1.3949 3.666 0.000452
## V_Selfdirection -6.5751 2.5995 -2.529 0.013471
## V_Hedonism 3.3061 1.7649 1.873 0.064832
## V_Power -2.3155 0.8912 -2.598 0.011224
## V_Security -0.8080 1.5027 -0.538 0.592325
## NationalityRussia:SexMale -18.7538 50.6171 -0.371 0.712024
## NationalityRussia:W_PositiveRel -0.4592 0.4769 -0.963 0.338628
## NationalityRussia:W_Autonomy -0.6518 0.4089 -1.594 0.114994
## NationalityRussia:W_PurposeLife 0.7801 0.3166 2.464 0.015976
## NationalityRussia:V_Conformity 5.4521 2.9455 1.851 0.068008
## NationalityRussia:V_benevolence -7.5792 3.2128 -2.359 0.020856
## NationalityRussia:V_Selfdirection 5.0932 3.7664 1.352 0.180248
## NationalityRussia:V_Hedonism 0.1551 2.6270 0.059 0.953060
## SexMale:W_PositiveRel -1.6438 0.7834 -2.098 0.039161
## SexMale:W_Autonomy -0.8274 0.4083 -2.026 0.046194
## SexMale:W_PurposeLife 0.4968 0.3105 1.600 0.113651
## SexMale:V_Selfdirection 5.3484 3.6414 1.469 0.145962
## SexMale:V_Hedonism -7.7637 3.0858 -2.516 0.013955
## SexMale:V_Security 3.2621 2.1162 1.541 0.127308
## NationalityRussia:SexMale:W_PositiveRel 1.0067 0.8188 1.229 0.222629
## NationalityRussia:SexMale:V_Selfdirection -9.9038 5.5641 -1.780 0.079030
## NationalityRussia:SexMale:V_Hedonism 6.4499 4.2479 1.518 0.133015
##
## (Intercept)
## NationalityRussia
## SexMale *
## W_PositiveRel
## W_Autonomy
## W_EnvMastery .
## W_PersGrowth
## W_PurposeLife
## V_Conformity
## V_benevolence *
## V_universalism ***
## V_Selfdirection *
## V_Hedonism .
## V_Power *
## V_Security
## NationalityRussia:SexMale
## NationalityRussia:W_PositiveRel
## NationalityRussia:W_Autonomy
## NationalityRussia:W_PurposeLife *
## NationalityRussia:V_Conformity .
## NationalityRussia:V_benevolence *
## NationalityRussia:V_Selfdirection
## NationalityRussia:V_Hedonism
## SexMale:W_PositiveRel *
## SexMale:W_Autonomy *
## SexMale:W_PurposeLife
## SexMale:V_Selfdirection
## SexMale:V_Hedonism *
## SexMale:V_Security
## NationalityRussia:SexMale:W_PositiveRel
## NationalityRussia:SexMale:V_Selfdirection .
## NationalityRussia:SexMale:V_Hedonism
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.12 on 77 degrees of freedom
## Multiple R-squared: 0.5418, Adjusted R-squared: 0.3573
## F-statistic: 2.937 on 31 and 77 DF, p-value: 0.00006898
lm(EnvIdentity ~ Nationality : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
Sex : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ W_PositiveRel + W_Autonomy + W_EnvMastery +
## W_PersGrowth + W_PurposeLife + V_Conformity + V_benevolence +
## V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
## V_Achievement + V_Power + W_EnvMastery:Nationality + V_benevolence:Nationality +
## V_Stimulation:Nationality + V_Achievement:Nationality + W_PositiveRel:Sex +
## W_Autonomy:Sex + W_PersGrowth:Sex + W_PurposeLife:Sex + V_Hedonism:Sex +
## V_Power:Sex, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -28.0945 -6.0592 0.5671 5.6963 19.3593
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 58.012039 17.625649 3.291 0.00145 **
## W_PositiveRel 0.006831 0.228332 0.030 0.97620
## W_Autonomy 0.027647 0.254295 0.109 0.91368
## W_EnvMastery -0.700392 0.255200 -2.744 0.00739 **
## W_PersGrowth 0.098038 0.294254 0.333 0.73982
## W_PurposeLife -0.166563 0.186134 -0.895 0.37339
## V_Conformity 3.818401 1.591056 2.400 0.01858 *
## V_benevolence 5.858784 2.059554 2.845 0.00557 **
## V_universalism 4.048200 1.406086 2.879 0.00504 **
## V_Selfdirection -3.972379 1.679431 -2.365 0.02029 *
## V_Stimulation -1.784736 1.246287 -1.432 0.15580
## V_Hedonism 3.339069 1.343102 2.486 0.01487 *
## V_Achievement -1.045728 1.592607 -0.657 0.51320
## V_Power -4.020971 1.286262 -3.126 0.00242 **
## W_EnvMastery:NationalityRussia 0.432833 0.204034 2.121 0.03680 *
## V_benevolence:NationalityRussia -7.477614 2.514506 -2.974 0.00383 **
## V_Stimulation:NationalityRussia 2.330988 1.662958 1.402 0.16464
## V_Achievement:NationalityRussia 3.902936 2.109337 1.850 0.06774 .
## W_PositiveRel:SexMale -0.489528 0.375394 -1.304 0.19574
## W_Autonomy:SexMale -0.627866 0.343485 -1.828 0.07107 .
## W_PersGrowth:SexMale 0.547982 0.320280 1.711 0.09074 .
## W_PurposeLife:SexMale 0.569176 0.244139 2.331 0.02210 *
## V_Hedonism:SexMale -2.745043 1.915181 -1.433 0.15544
## V_Power:SexMale 3.317726 1.566000 2.119 0.03704 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.19 on 85 degrees of freedom
## Multiple R-squared: 0.4874, Adjusted R-squared: 0.3487
## F-statistic: 3.514 on 23 and 85 DF, p-value: 0.00001265
lm(EnvIdentity ~ Nationality * (Sex + W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security)
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ Nationality + Sex + W_Autonomy + W_PurposeLife +
## V_benevolence + V_universalism + V_Power + V_Security + Nationality:Sex +
## Nationality:W_Autonomy + Nationality:W_PurposeLife + Nationality:V_benevolence +
## Nationality:V_Power, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -30.1457 -6.1668 -0.3201 6.4459 22.8636
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 32.3213 18.7704 1.722 0.08834 .
## NationalityRussia 13.9902 23.1531 0.604 0.54712
## SexMale -3.8839 2.9100 -1.335 0.18517
## W_Autonomy 0.1252 0.3080 0.407 0.68524
## W_PurposeLife -0.2787 0.1474 -1.891 0.06167 .
## V_benevolence 3.0434 1.7133 1.776 0.07887 .
## V_universalism 3.9805 1.2557 3.170 0.00205 **
## V_Power -2.9982 1.0605 -2.827 0.00573 **
## V_Security 2.2732 1.1119 2.044 0.04367 *
## NationalityRussia:SexMale 8.4953 4.7396 1.792 0.07625 .
## NationalityRussia:W_Autonomy -0.7192 0.3791 -1.897 0.06085 .
## NationalityRussia:W_PurposeLife 0.5842 0.2490 2.346 0.02103 *
## NationalityRussia:V_benevolence -3.0114 2.1643 -1.391 0.16736
## NationalityRussia:V_Power 2.9401 1.4154 2.077 0.04049 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.57 on 95 degrees of freedom
## Multiple R-squared: 0.3842, Adjusted R-squared: 0.2999
## F-statistic: 4.56 on 13 and 95 DF, p-value: 0.000005199
lm(EnvIdentity ~ Sex * (Nationality + W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security)
, Data2) %>%
step(direction = "both", steps = 1000000, trace = 0) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ Sex + Nationality + W_PositiveRel +
## W_Autonomy + W_PurposeLife + V_Conformity + V_universalism +
## V_Power + Sex:Nationality + Sex:W_PositiveRel + Sex:W_Autonomy +
## Sex:W_PurposeLife, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -26.7975 -5.9185 0.7032 7.1128 20.4659
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 46.66680 12.89172 3.620 0.000473 ***
## SexMale 38.06551 22.84752 1.666 0.098960 .
## NationalityRussia 6.71984 3.39313 1.980 0.050518 .
## W_PositiveRel 0.10338 0.22083 0.468 0.640746
## W_Autonomy -0.09425 0.23840 -0.395 0.693466
## W_PurposeLife -0.31613 0.16883 -1.872 0.064181 .
## V_Conformity 3.62655 1.32331 2.741 0.007316 **
## V_universalism 4.75743 1.12230 4.239 0.0000516 ***
## V_Power -2.24012 0.81389 -2.752 0.007075 **
## SexMale:NationalityRussia 13.31710 5.77925 2.304 0.023362 *
## SexMale:W_PositiveRel -0.80955 0.34870 -2.322 0.022366 *
## SexMale:W_Autonomy -0.72919 0.35306 -2.065 0.041586 *
## SexMale:W_PurposeLife 0.64801 0.25057 2.586 0.011207 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.59 on 96 degrees of freedom
## Multiple R-squared: 0.3752, Adjusted R-squared: 0.2971
## F-statistic: 4.804 on 12 and 96 DF, p-value: 0.000004006
Задаем полную и минимальную модель
start_model <- lm(EnvIdentity ~ 1, data = Data2)
start_model_1 <- lm(EnvIdentity ~ V_Power:Nationality, data = Data2)
start_model_2 <- lm(EnvIdentity ~ Nationality : (W_PersGrowth + W_PurposeLife + V_Power)
, data = Data2)
full_model <- lm(EnvIdentity ~ Nationality : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
Sex : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security, data = Data2)
full_model_1 <- lm(EnvIdentity ~ Sex + Nationality + Nationality : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
Sex : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security, data = Data2)
full_model_1N <- lm(EnvIdentity ~ Nationality + Nationality : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
Sex : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security, data = Data2)
full_model_1S <- lm(EnvIdentity ~ Sex + Nationality : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
Sex : (W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security) +
W_PositiveRel + W_Autonomy +
W_EnvMastery + W_PersGrowth + W_PurposeLife + W_SelfAcceptance +
WellBeing + V_Conformity + V_tradition + V_benevolence +
V_universalism + V_Selfdirection + V_Stimulation + V_Hedonism +
V_Achievement + V_Power + V_Security, data = Data2)
Поиск моделей в направлении “forward”
step(start_model, scope = list(lower = start_model, upper = full_model), direction = "forward") %>%
summary()
## Start: AIC=553.85
## EnvIdentity ~ 1
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3359.8 13866 532.20
## + V_benevolence 1 2063.5 15162 541.94
## + V_Conformity 1 1136.8 16089 548.40
## + V_tradition 1 1129.1 16097 548.46
## + V_Security 1 1068.7 16157 548.86
## + V_Selfdirection 1 742.5 16483 551.04
## <none> 17226 553.85
## + V_Achievement 1 217.8 17008 554.46
## + V_Hedonism 1 210.1 17016 554.51
## + W_Autonomy 1 189.4 17036 554.64
## + V_Stimulation 1 173.4 17052 554.74
## + W_SelfAcceptance 1 145.6 17080 554.92
## + W_EnvMastery 1 88.4 17137 555.28
## + W_PersGrowth 1 82.7 17143 555.32
## + V_Power 1 63.4 17162 555.44
## + W_PositiveRel 1 63.3 17162 555.44
## + WellBeing 1 7.0 17219 555.80
## + W_PurposeLife 1 4.9 17221 555.81
##
## Step: AIC=532.2
## EnvIdentity ~ V_universalism
##
## Df Sum of Sq RSS AIC
## + V_universalism:Sex 1 299.676 13566 531.81
## + W_PersGrowth 1 294.195 13572 531.86
## + V_universalism:Nationality 1 265.546 13600 532.09
## <none> 13866 532.20
## + V_benevolence 1 144.606 13721 533.05
## + W_PurposeLife 1 95.261 13771 533.44
## + V_Stimulation 1 80.701 13785 533.56
## + V_Security 1 79.682 13786 533.57
## + V_Selfdirection 1 76.249 13790 533.59
## + V_Power 1 71.691 13794 533.63
## + W_Autonomy 1 60.408 13805 533.72
## + W_SelfAcceptance 1 51.447 13814 533.79
## + V_tradition 1 46.258 13820 533.83
## + V_Conformity 1 44.186 13822 533.85
## + W_EnvMastery 1 31.016 13835 533.95
## + W_PositiveRel 1 25.092 13841 534.00
## + V_Hedonism 1 8.175 13858 534.13
## + WellBeing 1 3.096 13863 534.17
## + V_Achievement 1 0.042 13866 534.20
##
## Step: AIC=531.81
## EnvIdentity ~ V_universalism + V_universalism:Sex
##
## Df Sum of Sq RSS AIC
## + V_universalism:Nationality 1 257.476 13309 531.73
## + W_PersGrowth 1 248.933 13317 531.80
## <none> 13566 531.81
## + V_benevolence 1 149.748 13416 532.60
## + V_Security 1 111.153 13455 532.92
## + V_Power 1 100.335 13466 533.01
## + W_PurposeLife 1 70.605 13496 533.25
## + V_Stimulation 1 57.711 13508 533.35
## + V_tradition 1 46.276 13520 533.44
## + V_Conformity 1 40.499 13526 533.49
## + V_Selfdirection 1 31.915 13534 533.56
## + W_Autonomy 1 20.663 13546 533.65
## + W_SelfAcceptance 1 17.853 13548 533.67
## + W_EnvMastery 1 15.190 13551 533.69
## + W_PositiveRel 1 14.980 13551 533.69
## + V_Hedonism 1 9.895 13556 533.73
## + V_Achievement 1 0.189 13566 533.81
## + WellBeing 1 0.117 13566 533.81
##
## Step: AIC=531.73
## EnvIdentity ~ V_universalism + V_universalism:Sex + V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## + W_Autonomy 1 338.14 12970 530.92
## <none> 13309 531.73
## + V_benevolence 1 225.39 13083 531.86
## + WellBeing 1 172.71 13136 532.30
## + V_tradition 1 170.17 13138 532.32
## + V_Security 1 130.05 13179 532.66
## + W_PositiveRel 1 102.54 13206 532.88
## + W_PurposeLife 1 98.67 13210 532.91
## + V_Power 1 98.42 13210 532.92
## + V_Conformity 1 91.65 13217 532.97
## + W_EnvMastery 1 81.22 13228 533.06
## + W_SelfAcceptance 1 46.31 13262 533.35
## + V_Selfdirection 1 41.64 13267 533.38
## + V_Hedonism 1 19.94 13289 533.56
## + W_PersGrowth 1 18.30 13290 533.58
## + V_Stimulation 1 15.54 13293 533.60
## + V_Achievement 1 0.01 13309 533.73
##
## Step: AIC=530.92
## EnvIdentity ~ V_universalism + W_Autonomy + V_universalism:Sex +
## V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## + V_benevolence 1 237.902 12733 530.90
## <none> 12970 530.92
## + V_tradition 1 212.262 12758 531.12
## + V_Conformity 1 155.364 12815 531.61
## + V_Power 1 106.662 12864 532.02
## + W_PersGrowth 1 78.361 12892 532.26
## + V_Security 1 71.825 12899 532.32
## + W_PurposeLife 1 67.856 12903 532.35
## + V_Selfdirection 1 53.469 12917 532.47
## + W_PositiveRel 1 49.032 12922 532.51
## + V_Stimulation 1 31.987 12938 532.65
## + WellBeing 1 29.451 12941 532.67
## + W_Autonomy:Nationality 1 20.949 12950 532.74
## + W_SelfAcceptance 1 17.079 12954 532.78
## + W_EnvMastery 1 10.354 12960 532.83
## + W_Autonomy:Sex 1 1.913 12969 532.90
## + V_Achievement 1 0.455 12970 532.92
## + V_Hedonism 1 0.388 12970 532.92
##
## Step: AIC=530.9
## EnvIdentity ~ V_universalism + W_Autonomy + V_benevolence + V_universalism:Sex +
## V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## <none> 12733 530.90
## + V_Selfdirection 1 179.025 12554 531.36
## + V_Power 1 159.910 12573 531.53
## + V_tradition 1 109.967 12623 531.96
## + W_PurposeLife 1 101.367 12631 532.03
## + W_PositiveRel 1 91.607 12641 532.12
## + WellBeing 1 83.689 12649 532.18
## + V_Stimulation 1 75.014 12658 532.26
## + W_SelfAcceptance 1 60.264 12672 532.39
## + W_PersGrowth 1 50.222 12682 532.47
## + V_Achievement 1 49.578 12683 532.48
## + V_benevolence:Sex 1 42.491 12690 532.54
## + V_Conformity 1 40.374 12692 532.56
## + W_EnvMastery 1 34.482 12698 532.61
## + V_benevolence:Nationality 1 31.074 12702 532.64
## + V_Security 1 25.887 12707 532.68
## + W_Autonomy:Sex 1 23.884 12709 532.70
## + W_Autonomy:Nationality 1 20.894 12712 532.72
## + V_Hedonism 1 5.654 12727 532.85
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + W_Autonomy + V_benevolence +
## V_universalism:Sex + V_universalism:Nationality, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -28.4955 -6.1408 0.0724 6.3700 20.7577
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 44.1388 9.7930 4.507 0.0000174 ***
## V_universalism 3.8811 1.2019 3.229 0.00167 **
## W_Autonomy -0.2771 0.1645 -1.684 0.09517 .
## V_benevolence 1.8246 1.3153 1.387 0.16836
## V_universalism:SexMale -0.4890 0.4538 -1.078 0.28375
## V_universalism:NationalityRussia 1.3442 0.5769 2.330 0.02175 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.12 on 103 degrees of freedom
## Multiple R-squared: 0.2608, Adjusted R-squared: 0.225
## F-statistic: 7.269 on 5 and 103 DF, p-value: 0.00000723
Поиск моделей в направлении “both”
step(start_model, scope = list(lower = start_model, upper = full_model), direction = "both") %>%
summary()
## Start: AIC=553.85
## EnvIdentity ~ 1
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3359.8 13866 532.20
## + V_benevolence 1 2063.5 15162 541.94
## + V_Conformity 1 1136.8 16089 548.40
## + V_tradition 1 1129.1 16097 548.46
## + V_Security 1 1068.7 16157 548.86
## + V_Selfdirection 1 742.5 16483 551.04
## <none> 17226 553.85
## + V_Achievement 1 217.8 17008 554.46
## + V_Hedonism 1 210.1 17016 554.51
## + W_Autonomy 1 189.4 17036 554.64
## + V_Stimulation 1 173.4 17052 554.74
## + W_SelfAcceptance 1 145.6 17080 554.92
## + W_EnvMastery 1 88.4 17137 555.28
## + W_PersGrowth 1 82.7 17143 555.32
## + V_Power 1 63.4 17162 555.44
## + W_PositiveRel 1 63.3 17162 555.44
## + WellBeing 1 7.0 17219 555.80
## + W_PurposeLife 1 4.9 17221 555.81
##
## Step: AIC=532.2
## EnvIdentity ~ V_universalism
##
## Df Sum of Sq RSS AIC
## + V_universalism:Sex 1 299.7 13566 531.81
## + W_PersGrowth 1 294.2 13572 531.86
## + V_universalism:Nationality 1 265.5 13600 532.09
## <none> 13866 532.20
## + V_benevolence 1 144.6 13721 533.05
## + W_PurposeLife 1 95.3 13771 533.44
## + V_Stimulation 1 80.7 13785 533.56
## + V_Security 1 79.7 13786 533.57
## + V_Selfdirection 1 76.2 13790 533.59
## + V_Power 1 71.7 13794 533.63
## + W_Autonomy 1 60.4 13805 533.72
## + W_SelfAcceptance 1 51.4 13814 533.79
## + V_tradition 1 46.3 13820 533.83
## + V_Conformity 1 44.2 13822 533.85
## + W_EnvMastery 1 31.0 13835 533.95
## + W_PositiveRel 1 25.1 13841 534.00
## + V_Hedonism 1 8.2 13858 534.13
## + WellBeing 1 3.1 13863 534.17
## + V_Achievement 1 0.0 13866 534.20
## - V_universalism 1 3359.8 17226 553.85
##
## Step: AIC=531.81
## EnvIdentity ~ V_universalism + V_universalism:Sex
##
## Df Sum of Sq RSS AIC
## + V_universalism:Nationality 1 257.476 13309 531.73
## + W_PersGrowth 1 248.933 13317 531.80
## <none> 13566 531.81
## - V_universalism:Sex 1 299.676 13866 532.20
## + V_benevolence 1 149.748 13416 532.60
## + V_Security 1 111.153 13455 532.92
## + V_Power 1 100.335 13466 533.01
## + W_PurposeLife 1 70.605 13496 533.25
## + V_Stimulation 1 57.711 13508 533.35
## + V_tradition 1 46.276 13520 533.44
## + V_Conformity 1 40.499 13526 533.49
## + V_Selfdirection 1 31.915 13534 533.56
## + W_Autonomy 1 20.663 13546 533.65
## + W_SelfAcceptance 1 17.853 13548 533.67
## + W_EnvMastery 1 15.190 13551 533.69
## + W_PositiveRel 1 14.980 13551 533.69
## + V_Hedonism 1 9.895 13556 533.73
## + V_Achievement 1 0.189 13566 533.81
## + WellBeing 1 0.117 13566 533.81
##
## Step: AIC=531.73
## EnvIdentity ~ V_universalism + V_universalism:Sex + V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## + W_Autonomy 1 338.14 12970 530.92
## <none> 13309 531.73
## - V_universalism:Nationality 1 257.48 13566 531.81
## + V_benevolence 1 225.39 13083 531.86
## - V_universalism:Sex 1 291.61 13600 532.09
## + WellBeing 1 172.71 13136 532.30
## + V_tradition 1 170.17 13138 532.32
## + V_Security 1 130.05 13179 532.66
## + W_PositiveRel 1 102.54 13206 532.88
## + W_PurposeLife 1 98.67 13210 532.91
## + V_Power 1 98.42 13210 532.92
## + V_Conformity 1 91.65 13217 532.97
## + W_EnvMastery 1 81.22 13228 533.06
## + W_SelfAcceptance 1 46.31 13262 533.35
## + V_Selfdirection 1 41.64 13267 533.38
## + V_Hedonism 1 19.94 13289 533.56
## + W_PersGrowth 1 18.30 13290 533.58
## + V_Stimulation 1 15.54 13293 533.60
## + V_Achievement 1 0.01 13309 533.73
##
## Step: AIC=530.92
## EnvIdentity ~ V_universalism + W_Autonomy + V_universalism:Sex +
## V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## - V_universalism:Sex 1 142.02 13113 530.11
## + V_benevolence 1 237.90 12733 530.90
## <none> 12970 530.92
## + V_tradition 1 212.26 12758 531.12
## + V_Conformity 1 155.36 12815 531.61
## - W_Autonomy 1 338.14 13309 531.73
## + V_Power 1 106.66 12864 532.02
## + W_PersGrowth 1 78.36 12892 532.26
## + V_Security 1 71.83 12899 532.32
## + W_PurposeLife 1 67.86 12903 532.35
## + V_Selfdirection 1 53.47 12917 532.47
## + W_PositiveRel 1 49.03 12922 532.51
## + V_Stimulation 1 31.99 12938 532.65
## + WellBeing 1 29.45 12941 532.67
## + W_Autonomy:Nationality 1 20.95 12950 532.74
## + W_SelfAcceptance 1 17.08 12954 532.78
## + W_EnvMastery 1 10.35 12960 532.83
## + W_Autonomy:Sex 1 1.91 12969 532.90
## + V_Achievement 1 0.45 12970 532.92
## + V_Hedonism 1 0.39 12970 532.92
## - V_universalism:Nationality 1 574.96 13546 533.65
##
## Step: AIC=530.11
## EnvIdentity ~ V_universalism + W_Autonomy + V_universalism:Nationality
##
## Df Sum of Sq RSS AIC
## <none> 13113 530.11
## + V_benevolence 1 236.39 12876 530.12
## + V_tradition 1 220.98 12892 530.25
## + V_Conformity 1 172.22 12940 530.67
## + V_universalism:Sex 1 142.02 12970 530.92
## + W_Autonomy:Sex 1 121.20 12991 531.10
## + W_PersGrowth 1 119.17 12993 531.11
## + V_Selfdirection 1 89.18 13023 531.36
## + V_Power 1 87.46 13025 531.38
## + W_PurposeLife 1 79.03 13034 531.45
## + W_PositiveRel 1 54.16 13058 531.66
## + V_Security 1 46.98 13066 531.72
## + V_Stimulation 1 46.66 13066 531.72
## + W_SelfAcceptance 1 32.98 13080 531.83
## + WellBeing 1 29.37 13083 531.86
## + W_Autonomy:Nationality 1 21.93 13091 531.93
## - W_Autonomy 1 487.72 13600 532.09
## + W_EnvMastery 1 0.54 13112 532.10
## + V_Achievement 1 0.38 13112 532.10
## + V_Hedonism 1 0.04 13112 532.11
## - V_universalism:Nationality 1 692.86 13805 533.72
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + W_Autonomy + V_universalism:Nationality,
## data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -29.5124 -7.4646 0.1875 6.4359 21.6983
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 49.3011 9.3499 5.273 0.000000724 ***
## V_universalism 4.8944 0.9497 5.154 0.000001207 ***
## W_Autonomy -0.3163 0.1600 -1.976 0.0508 .
## V_universalism:NationalityRussia 1.3335 0.5661 2.355 0.0204 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.18 on 105 degrees of freedom
## Multiple R-squared: 0.2388, Adjusted R-squared: 0.217
## F-statistic: 10.98 on 3 and 105 DF, p-value: 0.000002491
Поиск моделей в направлении “forward”
step(start_model_1, scope = list(lower = start_model_1, upper = full_model), direction = "forward") %>%
summary()
## Start: AIC=556.97
## EnvIdentity ~ V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3694.2 13511 531.37
## + V_benevolence 1 2256.1 14949 542.39
## + V_tradition 1 1264.8 15940 549.39
## + V_Conformity 1 1237.4 15968 549.58
## + V_Security 1 1102.1 16103 550.50
## + V_Selfdirection 1 752.8 16452 552.84
## <none> 17205 555.71
## + W_Autonomy 1 280.4 16924 555.92
## + V_Hedonism 1 210.5 16994 556.37
## + V_Achievement 1 203.5 17002 556.42
## + W_EnvMastery 1 198.6 17006 556.45
## + V_Stimulation 1 178.6 17026 556.58
## + W_SelfAcceptance 1 141.0 17064 556.82
## + V_Power 1 116.7 17088 556.97
## + W_PersGrowth 1 64.3 17141 557.31
## + W_PositiveRel 1 58.5 17146 557.34
## + W_PurposeLife 1 4.0 17201 557.69
## + WellBeing 1 0.8 17204 557.71
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_Power 1 333.49 13177 530.64
## + W_Autonomy 1 249.77 13261 531.33
## + V_universalism:Sex 1 246.20 13265 531.36
## <none> 13511 531.37
## + V_benevolence 1 188.38 13322 531.84
## + W_PurposeLife 1 136.57 13374 532.26
## + V_Selfdirection 1 127.97 13383 532.33
## + WellBeing 1 120.96 13390 532.39
## + V_Stimulation 1 99.65 13411 532.56
## + W_SelfAcceptance 1 95.23 13416 532.60
## + V_tradition 1 88.36 13422 532.65
## + V_Security 1 73.80 13437 532.77
## + V_Conformity 1 68.34 13442 532.82
## + W_PositiveRel 1 57.07 13454 532.91
## + W_PersGrowth 1 48.77 13462 532.97
## + W_EnvMastery 1 22.18 13489 533.19
## + V_Achievement 1 14.94 13496 533.25
## + V_universalism:Nationality 1 11.02 13500 533.28
## + V_Hedonism 1 3.90 13507 533.34
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_Conformity 1 439.73 12738 528.94
## + V_tradition 1 428.14 12749 529.04
## + W_Autonomy 1 423.39 12754 529.08
## + V_Security 1 408.94 12768 529.21
## + V_benevolence 1 311.14 12866 530.04
## + V_universalism:Sex 1 279.31 12898 530.31
## + WellBeing 1 271.90 12905 530.37
## <none> 13177 530.64
## + W_PositiveRel 1 161.70 13016 531.30
## + W_PurposeLife 1 145.00 13032 531.44
## + W_EnvMastery 1 108.44 13069 531.74
## + W_SelfAcceptance 1 100.12 13077 531.81
## + V_Hedonism 1 76.25 13101 532.01
## + V_Power:Sex 1 61.40 13116 532.14
## + V_universalism:Nationality 1 60.27 13117 532.14
## + V_Selfdirection 1 54.72 13122 532.19
## + V_Achievement 1 53.49 13124 532.20
## + V_Stimulation 1 4.08 13173 532.61
## + W_PersGrowth 1 0.43 13177 532.64
##
## Step: AIC=528.94
## EnvIdentity ~ V_universalism + V_Power + V_Conformity + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + W_Autonomy 1 613.83 12124 525.56
## + WellBeing 1 439.22 12298 527.12
## + W_EnvMastery 1 310.63 12427 528.25
## + V_Conformity:Sex 1 301.28 12436 528.34
## + V_universalism:Sex 1 283.05 12454 528.50
## + V_Conformity:Nationality 1 244.95 12493 528.83
## + W_PositiveRel 1 235.56 12502 528.91
## <none> 12738 528.94
## + V_Security 1 208.67 12529 529.14
## + W_SelfAcceptance 1 156.57 12581 529.60
## + W_PurposeLife 1 148.11 12589 529.67
## + V_tradition 1 142.23 12595 529.72
## + V_Selfdirection 1 124.62 12613 529.87
## + V_universalism:Nationality 1 111.80 12626 529.98
## + V_Hedonism 1 103.45 12634 530.06
## + V_benevolence 1 96.69 12641 530.11
## + V_Power:Sex 1 26.25 12711 530.72
## + W_PersGrowth 1 9.83 12728 530.86
## + V_Stimulation 1 0.52 12737 530.94
## + V_Achievement 1 0.00 12738 530.94
##
## Step: AIC=525.56
## EnvIdentity ~ V_universalism + V_Power + V_Conformity + W_Autonomy +
## V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## <none> 12124 525.56
## + V_Conformity:Sex 1 136.150 11988 526.33
## + V_Security 1 135.031 11989 526.34
## + V_universalism:Sex 1 126.075 11998 526.42
## + V_Selfdirection 1 123.382 12000 526.45
## + W_Autonomy:Sex 1 104.128 12020 526.62
## + V_tradition 1 90.481 12033 526.74
## + W_PurposeLife 1 87.987 12036 526.77
## + W_PositiveRel 1 80.598 12043 526.83
## + W_PersGrowth 1 75.011 12049 526.88
## + WellBeing 1 73.144 12051 526.90
## + W_EnvMastery 1 66.918 12057 526.96
## + W_SelfAcceptance 1 63.055 12061 526.99
## + V_benevolence 1 58.584 12065 527.03
## + V_Hedonism 1 36.750 12087 527.23
## + V_Conformity:Nationality 1 30.596 12093 527.29
## + V_Stimulation 1 5.836 12118 527.51
## + V_Power:Sex 1 3.305 12120 527.53
## + V_universalism:Nationality 1 2.866 12121 527.54
## + V_Achievement 1 1.317 12122 527.55
## + W_Autonomy:Nationality 1 0.540 12123 527.56
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + V_Power + V_Conformity +
## W_Autonomy + V_Power:Nationality, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -30.490 -6.941 0.090 6.870 22.556
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 49.0541 8.7884 5.582 0.000000194 ***
## V_universalism 4.3303 1.0770 4.021 0.000111 ***
## V_Power -2.6694 0.8890 -3.003 0.003359 **
## V_Conformity 3.0307 1.3098 2.314 0.022663 *
## W_Autonomy -0.3179 0.1392 -2.284 0.024446 *
## V_Power:NationalityRussia 2.9129 0.8242 3.534 0.000614 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.85 on 103 degrees of freedom
## Multiple R-squared: 0.2962, Adjusted R-squared: 0.262
## F-statistic: 8.669 on 5 and 103 DF, p-value: 0.0000006923
Поиск моделей в направлении “forward”
step(start_model_2, scope = list(lower = start_model_2, upper = full_model), direction = "forward") %>%
summary()
## Start: AIC=562.88
## EnvIdentity ~ Nationality:(W_PersGrowth + W_PurposeLife + V_Power)
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3359.8 13866 532.20
## + V_benevolence 1 2063.5 15162 541.94
## + V_Conformity 1 1136.8 16089 548.40
## + V_tradition 1 1129.1 16097 548.46
## + V_Security 1 1068.7 16157 548.86
## + V_Selfdirection 1 742.5 16483 551.04
## <none> 17226 553.85
## + V_Achievement 1 217.8 17008 554.46
## + V_Hedonism 1 210.1 17016 554.51
## + W_Autonomy 1 189.4 17036 554.64
## + V_Stimulation 1 173.4 17052 554.74
## + W_SelfAcceptance 1 145.6 17080 554.92
## + W_EnvMastery 1 88.4 17137 555.28
## + W_PersGrowth 1 82.7 17143 555.32
## + V_Power 1 63.4 17162 555.44
## + W_PositiveRel 1 63.3 17162 555.44
## + WellBeing 1 7.0 17219 555.80
## + W_PurposeLife 1 4.9 17221 555.81
##
## Step: AIC=536.18
## EnvIdentity ~ V_universalism + Nationality:W_PersGrowth + Nationality:W_PurposeLife +
## Nationality:V_Power
##
## Df Sum of Sq RSS AIC
## + V_universalism:Sex 1 299.676 13566 531.81
## + W_PersGrowth 1 294.195 13572 531.86
## + V_universalism:Nationality 1 265.546 13600 532.09
## <none> 13866 532.20
## + V_benevolence 1 144.606 13721 533.05
## + W_PurposeLife 1 95.261 13771 533.44
## + V_Stimulation 1 80.701 13785 533.56
## + V_Security 1 79.682 13786 533.57
## + V_Selfdirection 1 76.249 13790 533.59
## + V_Power 1 71.691 13794 533.63
## + W_Autonomy 1 60.408 13805 533.72
## + W_SelfAcceptance 1 51.447 13814 533.79
## + V_tradition 1 46.258 13820 533.83
## + V_Conformity 1 44.186 13822 533.85
## + W_EnvMastery 1 31.016 13835 533.95
## + W_PositiveRel 1 25.092 13841 534.00
## + V_Hedonism 1 8.175 13858 534.13
## + WellBeing 1 3.096 13863 534.17
## + V_Achievement 1 0.042 13866 534.20
##
## Step: AIC=536.73
## EnvIdentity ~ V_universalism + Nationality:W_PersGrowth + Nationality:W_PurposeLife +
## Nationality:V_Power + V_universalism:Sex
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + Nationality:W_PersGrowth +
## Nationality:W_PurposeLife + Nationality:V_Power + V_universalism:Sex,
## data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -26.872 -6.727 0.670 7.381 21.938
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 38.194836 13.361843 2.859 0.00518 **
## V_universalism 5.748294 1.053516 5.456 0.000000353 ***
## NationalityFrance:W_PersGrowth 0.114392 0.257532 0.444 0.65787
## NationalityRussia:W_PersGrowth -0.079496 0.276445 -0.288 0.77427
## NationalityFrance:W_PurposeLife -0.160089 0.142122 -1.126 0.26268
## NationalityRussia:W_PurposeLife 0.004213 0.219521 0.019 0.98473
## NationalityFrance:V_Power -1.417363 0.909928 -1.558 0.12247
## NationalityRussia:V_Power 0.486133 1.021542 0.476 0.63520
## V_universalism:SexMale -0.538818 0.465338 -1.158 0.24966
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.27 on 100 degrees of freedom
## Multiple R-squared: 0.262, Adjusted R-squared: 0.203
## F-statistic: 4.438 on 8 and 100 DF, p-value: 0.000126
Поиск моделей в направлении “forward”
step(start_model_1, scope = list(lower = start_model_1, upper = full_model_1), direction = "forward") %>%
summary()
## Start: AIC=556.97
## EnvIdentity ~ V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3359.8 13866 532.20
## + V_benevolence 1 2063.5 15162 541.94
## + V_Conformity 1 1136.8 16089 548.40
## + V_tradition 1 1129.1 16097 548.46
## + V_Security 1 1068.7 16157 548.86
## + V_Selfdirection 1 742.5 16483 551.04
## + Sex 1 704.1 16522 551.30
## <none> 17226 553.85
## + V_Achievement 1 217.8 17008 554.46
## + V_Hedonism 1 210.1 17016 554.51
## + W_Autonomy 1 189.4 17036 554.64
## + V_Stimulation 1 173.4 17052 554.74
## + W_SelfAcceptance 1 145.6 17080 554.92
## + W_EnvMastery 1 88.4 17137 555.28
## + W_PersGrowth 1 82.7 17143 555.32
## + V_Power 1 63.4 17162 555.44
## + W_PositiveRel 1 63.3 17162 555.44
## + WellBeing 1 7.0 17219 555.80
## + W_PurposeLife 1 4.9 17221 555.81
## + Nationality 1 0.0 17226 555.85
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + Nationality 1 323.51 13542 531.62
## + W_PersGrowth 1 294.20 13572 531.86
## + Sex 1 280.59 13585 531.97
## <none> 13866 532.20
## + V_benevolence 1 144.61 13721 533.05
## + W_PurposeLife 1 95.26 13771 533.44
## + V_Stimulation 1 80.70 13785 533.56
## + V_Security 1 79.68 13786 533.57
## + V_Selfdirection 1 76.25 13790 533.59
## + V_Power 1 71.69 13794 533.63
## + W_Autonomy 1 60.41 13805 533.72
## + W_SelfAcceptance 1 51.45 13814 533.79
## + V_tradition 1 46.26 13820 533.83
## + V_Conformity 1 44.19 13822 533.85
## + W_EnvMastery 1 31.02 13835 533.95
## + W_PositiveRel 1 25.09 13841 534.00
## + V_Hedonism 1 8.18 13858 534.13
## + WellBeing 1 3.10 13863 534.17
## + V_Achievement 1 0.04 13866 534.20
##
## Step: AIC=532.35
## EnvIdentity ~ V_universalism + Nationality + V_Power:Nationality
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + Nationality + V_Power:Nationality,
## data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -25.7661 -7.4689 -0.2818 6.8209 23.5597
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 35.4962 5.9030 6.013 0.0000000273 ***
## V_universalism 5.4079 0.9938 5.441 0.0000003535 ***
## NationalityRussia -2.3080 4.3311 -0.533 0.595
## NationalityFrance:V_Power -1.4959 0.9124 -1.639 0.104
## NationalityRussia:V_Power 0.7081 1.0161 0.697 0.487
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.24 on 104 degrees of freedom
## Multiple R-squared: 0.2371, Adjusted R-squared: 0.2078
## F-statistic: 8.081 on 4 and 104 DF, p-value: 0.0000103
Поиск моделей в направлении “forward”
step(start_model_1, scope = list(lower = start_model_1, upper = full_model_1N), direction = "forward") %>%
summary()
## Start: AIC=556.97
## EnvIdentity ~ V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3359.8 13866 532.20
## + V_benevolence 1 2063.5 15162 541.94
## + V_Conformity 1 1136.8 16089 548.40
## + V_tradition 1 1129.1 16097 548.46
## + V_Security 1 1068.7 16157 548.86
## + V_Selfdirection 1 742.5 16483 551.04
## <none> 17226 553.85
## + V_Achievement 1 217.8 17008 554.46
## + V_Hedonism 1 210.1 17016 554.51
## + W_Autonomy 1 189.4 17036 554.64
## + V_Stimulation 1 173.4 17052 554.74
## + W_SelfAcceptance 1 145.6 17080 554.92
## + W_EnvMastery 1 88.4 17137 555.28
## + W_PersGrowth 1 82.7 17143 555.32
## + V_Power 1 63.4 17162 555.44
## + W_PositiveRel 1 63.3 17162 555.44
## + WellBeing 1 7.0 17219 555.80
## + W_PurposeLife 1 4.9 17221 555.81
## + Nationality 1 0.0 17226 555.85
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + Nationality 1 323.51 13542 531.62
## + V_universalism:Sex 1 299.68 13566 531.81
## + W_PersGrowth 1 294.20 13572 531.86
## <none> 13866 532.20
## + V_benevolence 1 144.61 13721 533.05
## + W_PurposeLife 1 95.26 13771 533.44
## + V_Stimulation 1 80.70 13785 533.56
## + V_Security 1 79.68 13786 533.57
## + V_Selfdirection 1 76.25 13790 533.59
## + V_Power 1 71.69 13794 533.63
## + W_Autonomy 1 60.41 13805 533.72
## + W_SelfAcceptance 1 51.45 13814 533.79
## + V_tradition 1 46.26 13820 533.83
## + V_Conformity 1 44.19 13822 533.85
## + W_EnvMastery 1 31.02 13835 533.95
## + W_PositiveRel 1 25.09 13841 534.00
## + V_Hedonism 1 8.18 13858 534.13
## + WellBeing 1 3.10 13863 534.17
## + V_Achievement 1 0.04 13866 534.20
##
## Step: AIC=532.35
## EnvIdentity ~ V_universalism + Nationality + V_Power:Nationality
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + Nationality + V_Power:Nationality,
## data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -25.7661 -7.4689 -0.2818 6.8209 23.5597
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 35.4962 5.9030 6.013 0.0000000273 ***
## V_universalism 5.4079 0.9938 5.441 0.0000003535 ***
## NationalityRussia -2.3080 4.3311 -0.533 0.595
## NationalityFrance:V_Power -1.4959 0.9124 -1.639 0.104
## NationalityRussia:V_Power 0.7081 1.0161 0.697 0.487
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 11.24 on 104 degrees of freedom
## Multiple R-squared: 0.2371, Adjusted R-squared: 0.2078
## F-statistic: 8.081 on 4 and 104 DF, p-value: 0.0000103
Поиск моделей в направлении “forward”
step(start_model_1, scope = list(lower = start_model_1, upper = full_model_1S), direction = "forward") %>%
summary()
## Start: AIC=556.97
## EnvIdentity ~ V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_universalism 1 3694.2 13511 531.37
## + V_benevolence 1 2256.1 14949 542.39
## + V_tradition 1 1264.8 15940 549.39
## + V_Conformity 1 1237.4 15968 549.58
## + V_Security 1 1102.1 16103 550.50
## + V_Selfdirection 1 752.8 16452 552.84
## + Sex 1 692.1 16513 553.24
## <none> 17205 555.71
## + W_Autonomy 1 280.4 16924 555.92
## + V_Hedonism 1 210.5 16994 556.37
## + V_Achievement 1 203.5 17002 556.42
## + W_EnvMastery 1 198.6 17006 556.45
## + V_Stimulation 1 178.6 17026 556.58
## + W_SelfAcceptance 1 141.0 17064 556.82
## + V_Power 1 116.7 17088 556.97
## + W_PersGrowth 1 64.3 17141 557.31
## + W_PositiveRel 1 58.5 17146 557.34
## + W_PurposeLife 1 4.0 17201 557.69
## + WellBeing 1 0.8 17204 557.71
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_Power 1 333.49 13177 530.64
## + W_Autonomy 1 249.77 13261 531.33
## <none> 13511 531.37
## + Sex 1 220.13 13291 531.58
## + V_benevolence 1 188.38 13322 531.84
## + W_PurposeLife 1 136.57 13374 532.26
## + V_Selfdirection 1 127.97 13383 532.33
## + WellBeing 1 120.96 13390 532.39
## + V_Stimulation 1 99.65 13411 532.56
## + W_SelfAcceptance 1 95.23 13416 532.60
## + V_tradition 1 88.36 13422 532.65
## + V_Security 1 73.80 13437 532.77
## + V_Conformity 1 68.34 13442 532.82
## + W_PositiveRel 1 57.07 13454 532.91
## + W_PersGrowth 1 48.77 13462 532.97
## + W_EnvMastery 1 22.18 13489 533.19
## + V_Achievement 1 14.94 13496 533.25
## + V_universalism:Nationality 1 11.02 13500 533.28
## + V_Hedonism 1 3.90 13507 533.34
##
## Step: AIC=530.64
## EnvIdentity ~ V_universalism + V_Power + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + V_Conformity 1 439.73 12738 528.94
## + V_tradition 1 428.14 12749 529.04
## + W_Autonomy 1 423.39 12754 529.08
## + V_Security 1 408.94 12768 529.21
## + V_benevolence 1 311.14 12866 530.04
## + WellBeing 1 271.90 12905 530.37
## + Sex 1 251.10 12926 530.55
## <none> 13177 530.64
## + W_PositiveRel 1 161.70 13016 531.30
## + W_PurposeLife 1 145.00 13032 531.44
## + W_EnvMastery 1 108.44 13069 531.74
## + W_SelfAcceptance 1 100.12 13077 531.81
## + V_Hedonism 1 76.25 13101 532.01
## + V_universalism:Nationality 1 60.27 13117 532.14
## + V_Selfdirection 1 54.72 13122 532.19
## + V_Achievement 1 53.49 13124 532.20
## + V_Stimulation 1 4.08 13173 532.61
## + W_PersGrowth 1 0.43 13177 532.64
##
## Step: AIC=528.94
## EnvIdentity ~ V_universalism + V_Power + V_Conformity + V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## + W_Autonomy 1 613.83 12124 525.56
## + WellBeing 1 439.22 12298 527.12
## + W_EnvMastery 1 310.63 12427 528.25
## + V_Conformity:Nationality 1 244.95 12493 528.83
## + Sex 1 240.28 12497 528.87
## + W_PositiveRel 1 235.56 12502 528.91
## <none> 12738 528.94
## + V_Security 1 208.67 12529 529.14
## + W_SelfAcceptance 1 156.57 12581 529.60
## + W_PurposeLife 1 148.11 12589 529.67
## + V_tradition 1 142.23 12595 529.72
## + V_Selfdirection 1 124.62 12613 529.87
## + V_universalism:Nationality 1 111.80 12626 529.98
## + V_Hedonism 1 103.45 12634 530.06
## + V_benevolence 1 96.69 12641 530.11
## + W_PersGrowth 1 9.83 12728 530.86
## + V_Stimulation 1 0.52 12737 530.94
## + V_Achievement 1 0.00 12738 530.94
##
## Step: AIC=525.56
## EnvIdentity ~ V_universalism + V_Power + V_Conformity + W_Autonomy +
## V_Power:Nationality
##
## Df Sum of Sq RSS AIC
## <none> 12124 525.56
## + V_Security 1 135.031 11989 526.34
## + V_Selfdirection 1 123.382 12000 526.45
## + Sex 1 91.935 12032 526.73
## + V_tradition 1 90.481 12033 526.74
## + W_PurposeLife 1 87.987 12036 526.77
## + W_PositiveRel 1 80.598 12043 526.83
## + W_PersGrowth 1 75.011 12049 526.88
## + WellBeing 1 73.144 12051 526.90
## + W_EnvMastery 1 66.918 12057 526.96
## + W_SelfAcceptance 1 63.055 12061 526.99
## + V_benevolence 1 58.584 12065 527.03
## + V_Hedonism 1 36.750 12087 527.23
## + V_Conformity:Nationality 1 30.596 12093 527.29
## + V_Stimulation 1 5.836 12118 527.51
## + V_universalism:Nationality 1 2.866 12121 527.54
## + V_Achievement 1 1.317 12122 527.55
## + W_Autonomy:Nationality 1 0.540 12123 527.56
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + V_Power + V_Conformity +
## W_Autonomy + V_Power:Nationality, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -30.490 -6.941 0.090 6.870 22.556
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 49.0541 8.7884 5.582 0.000000194 ***
## V_universalism 4.3303 1.0770 4.021 0.000111 ***
## V_Power -2.6694 0.8890 -3.003 0.003359 **
## V_Conformity 3.0307 1.3098 2.314 0.022663 *
## W_Autonomy -0.3179 0.1392 -2.284 0.024446 *
## V_Power:NationalityRussia 2.9129 0.8242 3.534 0.000614 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.85 on 103 degrees of freedom
## Multiple R-squared: 0.2962, Adjusted R-squared: 0.262
## F-statistic: 8.669 on 5 and 103 DF, p-value: 0.0000006923
lm(formula = EnvIdentity ~ V_universalism + V_Power + V_Conformity +
W_Autonomy + V_Power:Nationality, data = Data2) %>%
summary()
##
## Call:
## lm(formula = EnvIdentity ~ V_universalism + V_Power + V_Conformity +
## W_Autonomy + V_Power:Nationality, data = Data2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -30.490 -6.941 0.090 6.870 22.556
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 49.0541 8.7884 5.582 0.000000194 ***
## V_universalism 4.3303 1.0770 4.021 0.000111 ***
## V_Power -2.6694 0.8890 -3.003 0.003359 **
## V_Conformity 3.0307 1.3098 2.314 0.022663 *
## W_Autonomy -0.3179 0.1392 -2.284 0.024446 *
## V_Power:NationalityRussia 2.9129 0.8242 3.534 0.000614 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 10.85 on 103 degrees of freedom
## Multiple R-squared: 0.2962, Adjusted R-squared: 0.262
## F-statistic: 8.669 on 5 and 103 DF, p-value: 0.0000006923
Data %>%
select(Nationality, EnvIdentity, V_universalism, V_Power, V_Conformity, W_Autonomy) %>% # Выбираем только нужные переменные
pivot_longer(
cols = c(V_universalism, V_Power, V_Conformity, W_Autonomy),
names_to = "Variable",
values_to = "Value"
) %>%
# Добавляем новые метки для переменных
mutate(Variable = recode(
Variable,
V_universalism = "Универсализм (ЦО)",
V_Power = "Власть (ЦО)",
V_Conformity = "Конформность (ЦО)",
W_Autonomy = "Автономность (ПБ)"
)) %>%
ggplot(aes(x = Value, y = EnvIdentity, col = Nationality)) +
geom_point(alpha = 0.7, size = 2) + # Добавление прозрачности и размера точек
geom_smooth(method = 'lm', se = FALSE, linewidth = 0.8) + # Линия тренда
facet_grid(cols = vars(Variable), scales = "free_x") + # Свободная шкала для X
labs(
#title = "Факторы регрессионной модели и их связь с идентификацией с природой",
x = "Значение переменной",
y = "Идентификация с природой",
col = "Страна",
caption = "Примечание: (ЦО) — Ценностные ориентации, (ПБ) — Психологическое благополучие"
) +
theme_light(base_size = 12) +
theme(
strip.text = element_text(size = 11, face = "bold"), # Оформление заголовков фасетов
legend.position = "bottom", # Перемещение легенды вниз
legend.title = element_text(size = 11, face = "bold"), # Форматирование заголовка легенды
legend.text = element_text(size = 10), # Форматирование текста легенды
plot.title = element_text(size = 14, face = "bold", hjust = 0.5), # Центровка и форматирование заголовка
plot.caption = element_text(size = 10, hjust = 0), # Оформление подписи
panel.grid.minor = element_blank() # Удаление мелкой сетки
)
## `geom_smooth()` using formula = 'y ~ x'