date()
## [1] "Sun Apr 21 00:33:59 2024"
options(scipen = 999)

Обозначения

Психологическое благополучие

  • positive relationship - позитивные отношения

  • autonomy - автономность

  • environmental mastery - управление средой

  • personal growth - личностный рост

  • purpose in life - цель в жизни

  • self-acceptance - самопринятие

  • pwb - psychologycal well-being - психологическое благополучие

ei - environmental identity - природная идентичность

Ценностные ориентации

  • conformity - конформность

  • tradition - традиции

  • benevolence - доброта

  • universalism - универсализм

  • selfdirection - самостоятельность

  • stimulation - стимуляция

  • hedonism - гедонизм

  • achievement - достижения

  • power - власть

  • security - безопасность

  • pwb - psychologycal well-being

возраст, пол, в каком населенном пункте выросли (child), в каком живут сейчас (adult) # Libraries

library(kableExtra)
library(reshape2)
library(ggplot2)
library(dplyr)
## 
## Присоединяю пакет: 'dplyr'
## Следующий объект скрыт от 'package:kableExtra':
## 
##     group_rows
## Следующие объекты скрыты от 'package:stats':
## 
##     filter, lag
## Следующие объекты скрыты от 'package:base':
## 
##     intersect, setdiff, setequal, union
library(stats)
library(tidyr)
## 
## Присоединяю пакет: 'tidyr'
## Следующий объект скрыт от 'package:reshape2':
## 
##     smiths
library(readxl)
thes <- read_excel("thes.xlsx")

dividing samples into rus and fr

histograms and boxplots

Гистограммы для общей выборки

# Определение границ для осей x и y
x_min <- min(min(thes$positiverelat), min(thes$autonomy), min(thes$envmattery), min(thes$persgrowth), min(thes$purposeinlife), min(thes$selfacceptance))
x_max <- max(max(thes$positiverelat), max(thes$autonomy), max(thes$envmattery), max(thes$persgrowth), max(thes$purposeinlife), max(thes$selfacceptance))
y_max <- max(hist(thes$positiverelat)$counts, hist(thes$autonomy)$counts, hist(thes$envmattery)$counts, hist(thes$persgrowth)$counts, hist(thes$purposeinlife)$counts, hist(thes$selfacceptance)$counts)

# Разделение графиков на 2 строки и 3 столбца
par(mfrow = c(2, 3))

# Создание гистограмм для каждой переменной
hist(thes$positiverelat, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "positiverelat", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$autonomy, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "autonomy", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$envmattery, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "envmattery", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$persgrowth, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "persgrowth", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$purposeinlife, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "purposeinlife", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$selfacceptance, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfacceptance", xlab = "Значения", ylab = "Частота", breaks = 15)

# Возврат к одиночному графику
par(mfrow = c(1, 1))


x_min <- min(min(thes$conformity), min(thes$tradition), min(thes$benevolence), min(thes$universalism), min(thes$selfdirection), min(thes$stimulation), min(thes$hedonism), min(thes$achievement), min(thes$power), min(thes$security))
x_max <- max(max(thes$conformity), max(thes$tradition), max(thes$benevolence), max(thes$universalism), max(thes$selfdirection), max(thes$stimulation), max(thes$hedonism), max(thes$achievement), max(thes$power), max(thes$security))
y_max <- max(hist(thes$conformity)$counts, hist(thes$tradition)$counts, hist(thes$benevolence)$counts, hist(thes$universalism)$counts, hist(thes$selfdirection)$counts, hist(thes$stimulation)$counts, hist(thes$hedonism)$counts, hist(thes$achievement)$counts, hist(thes$power)$counts, hist(thes$security)$counts)

par(mfrow = c(2, 5))

hist(thes$conformity, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "conformity", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$tradition, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "tradition", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$benevolence, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "benevolence", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$universalism, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "universalism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$selfdirection, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfdirection", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$stimulation, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "stimulation", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$hedonism, col = "cyan", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "hedonism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$achievement, col = "magenta", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "achievement", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$power, col = "brown", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "power", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(thes$security, col = "gray", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "security", xlab = "Значения", ylab = "Частота", breaks = 15)

par(mfrow = c(1, 1))

Боксплоты для общей выборки

par(mfrow = c(1, 6))
max_val <- max(thes$positiverelat, thes$autonomy, thes$envmattery, thes$persgrowth, thes$purposeinlife, thes$selfacceptance, na.rm = TRUE)
ylim <- c(0, max_val)
boxplot(thes$positiverelat, main = "positiverelat", ylim = ylim)
boxplot(thes$autonomy, main = "autonomym", ylim = ylim)
boxplot(thes$envmattery, main = "envmattery", ylim = ylim)
boxplot(thes$persgrowth, main = "persgrowt", ylim = ylim)
boxplot(thes$purposeinlife, main = "purposeinlife", ylim = ylim)
boxplot(thes$selfacceptance, main = "selfacceptanc", ylim = ylim)

par(mfrow = c(1, 5))
max_val <- max(c(thes$conformity, thes$tradition, thes$benevolence, thes$universalism, thes$selfdirection, thes$stimulation, thes$hedonism, thes$achievement, thes$power, thes$security), na.rm = TRUE)
ylim <- c(-1, max_val)
boxplot(thes$conformity, main = "conformity", ylim = ylim)
boxplot(thes$tradition, main = "tradition", ylim = ylim)
boxplot(thes$benevolence, main = "benevolence", ylim = ylim)
boxplot(thes$universalism, main = "universalism", ylim = ylim)
boxplot(thes$selfdirection, main = "selfdirection", ylim = ylim)

boxplot(thes$stimulation, main = "stimulation", ylim = ylim)
boxplot(thes$hedonism, main = "hedonism", ylim = ylim)
boxplot(thes$achievement, main = "achievement", ylim = ylim)
boxplot(thes$power, main = "power", ylim = ylim)
boxplot(thes$security, main = "security", ylim = ylim)

Гистограммы для фр выборки

x_min <- min(min(french$positiverelat), min(french$autonomy), min(french$envmattery), min(french$persgrowth), min(french$purposeinlife), min(french$selfacceptance))
x_max <- max(max(french$positiverelat), max(french$autonomy), max(french$envmattery), max(french$persgrowth), max(french$purposeinlife), max(french$selfacceptance))
y_max <- max(hist(french$positiverelat)$counts, hist(french$autonomy)$counts, hist(french$envmattery)$counts, hist(french$persgrowth)$counts, hist(french$purposeinlife)$counts, hist(french$selfacceptance)$counts)

par(mfrow = c(2, 3))

hist(french$positiverelat, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "positiverelat", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$autonomy, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "autonomy", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$envmattery, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "envmattery", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$persgrowth, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "persgrowth", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$purposeinlife, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "purposeinlife", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$selfacceptance, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfacceptance", xlab = "Значения", ylab = "Частота", breaks = 15)

par(mfrow = c(1, 1))


x_min <- min(min(french$conformity), min(french$tradition), min(french$benevolence), min(french$universalism), min(french$selfdirection), min(french$stimulation), min(french$hedonism), min(french$achievement), min(french$power), min(french$security))
x_max <- max(max(french$conformity), max(french$tradition), max(french$benevolence), max(french$universalism), max(french$selfdirection), max(french$stimulation), max(french$hedonism), max(french$achievement), max(french$power), max(french$security))
y_max <- max(hist(french$conformity)$counts, hist(french$tradition)$counts, hist(french$benevolence)$counts, hist(french$universalism)$counts, hist(french$selfdirection)$counts, hist(french$stimulation)$counts, hist(french$hedonism)$counts, hist(french$achievement)$counts, hist(french$power)$counts, hist(french$security)$counts)

par(mfrow = c(2, 5))

hist(french$conformity, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "conformity", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$tradition, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "tradition", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$benevolence, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "benevolence", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$universalism, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "universalism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$selfdirection, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfdirection", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$stimulation, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "stimulation", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$hedonism, col = "cyan", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "hedonism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$achievement, col = "magenta", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "achievement", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$power, col = "brown", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "power", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(french$security, col = "gray", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "security", xlab = "Значения", ylab = "Частота", breaks = 15)

par(mfrow = c(1, 1))

Боксплоты для фр выборки

par(mfrow = c(1, 6))
max_val <- max(french$positiverelat, french$autonomy, french$envmattery, french$persgrowth, french$purposeinlife, french$selfacceptance)
ylim <- c(0, max_val)
boxplot(french$positiverelat, main = "positiverelat", ylim = ylim)
boxplot(french$autonomy, main = "autonomym", ylim = ylim)
boxplot(french$envmattery, main = "envmattery", ylim = ylim)
boxplot(french$persgrowth, main = "persgrowt", ylim = ylim)
boxplot(french$purposeinlife, main = "purposeinlife", ylim = ylim)
boxplot(french$selfacceptance, main = "selfacceptanc", ylim = ylim)

par(mfrow = c(1, 5))
max_val <- max(c(french$conformity, french$tradition, french$benevolence, french$universalism, french$selfdirection, french$stimulation, french$hedonism, french$achievement, french$power, french$security))
ylim <- c(-1, max_val)
boxplot(french$conformity, main = "frenchconformity", ylim = ylim)
boxplot(french$tradition, main = "frenchtradition", ylim = ylim)
boxplot(french$benevolence, main = "frenchbenevolence", ylim = ylim)
boxplot(french$universalism, main = "frenchuniversalism", ylim = ylim)
boxplot(french$selfdirection, main = "frenchselfdirection", ylim = ylim)

boxplot(french$stimulation, main = "frenchstimulation", ylim = ylim)
boxplot(french$hedonism, main = "frenchhedonism", ylim = ylim)
boxplot(french$achievement, main = "frenchachievement", ylim = ylim)
boxplot(french$power, main = "frenchpower", ylim = ylim)
boxplot(french$security, main = "frenchsecurity", ylim = ylim)

Гистограммы для рус выборки

x_min <- min(min(russian$positiverelat), min(russian$autonomy), min(russian$envmattery), min(russian$persgrowth), min(russian$purposeinlife), min(russian$selfacceptance))
x_max <- max(max(russian$positiverelat), max(russian$autonomy), max(russian$envmattery), max(russian$persgrowth), max(russian$purposeinlife), max(russian$selfacceptance))
y_max <- max(hist(russian$positiverelat)$counts, hist(russian$autonomy)$counts, hist(russian$envmattery)$counts, hist(russian$persgrowth)$counts, hist(russian$purposeinlife)$counts, hist(russian$selfacceptance)$counts)

par(mfrow = c(2, 3))

hist(russian$positiverelat, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "positiverelat", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$autonomy, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "autonomy", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$envmattery, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "envmattery", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$persgrowth, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "persgrowth", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$purposeinlife, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "purposeinlife", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$selfacceptance, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfacceptance", xlab = "Значения", ylab = "Частота", breaks = 15)

par(mfrow = c(1, 1))


x_min <- min(min(russian$conformity), min(russian$tradition), min(russian$benevolence), min(russian$universalism), min(russian$selfdirection), min(russian$stimulation), min(russian$hedonism), min(russian$achievement), min(russian$power), min(russian$security))
x_max <- max(max(russian$conformity), max(russian$tradition), max(russian$benevolence), max(russian$universalism), max(russian$selfdirection), max(russian$stimulation), max(russian$hedonism), max(russian$achievement), max(russian$power), max(russian$security))
y_max <- max(hist(russian$conformity)$counts, hist(russian$tradition)$counts, hist(russian$benevolence)$counts, hist(russian$universalism)$counts, hist(russian$selfdirection)$counts, hist(russian$stimulation)$counts, hist(russian$hedonism)$counts, hist(russian$achievement)$counts, hist(russian$power)$counts, hist(russian$security)$counts)

par(mfrow = c(2, 5))

hist(russian$conformity, col = "blue", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "conformity", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$tradition, col = "red", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "tradition", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$benevolence, col = "green", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "benevolence", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$universalism, col = "yellow", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "universalism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$selfdirection, col = "orange", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "selfdirection", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$stimulation, col = "purple", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "stimulation", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$hedonism, col = "cyan", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "hedonism", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$achievement, col = "magenta", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "achievement", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$power, col = "brown", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "power", xlab = "Значения", ylab = "Частота", breaks = 15)
hist(russian$security, col = "gray", xlim = c(x_min, x_max), ylim = c(0, y_max), main = "security", xlab = "Значения", ylab = "Частота", breaks = 15)

par(mfrow = c(1, 1))

Боксплоты для рус выборки

par(mfrow = c(1, 6))
max_val <- max(russian$positiverelat, russian$autonomy, russian$envmattery, russian$persgrowth, russian$purposeinlife, russian$selfacceptance)
ylim <- c(0, max_val)
boxplot(russian$positiverelat, main = "russianpositiverelat", ylim = ylim)
boxplot(russian$autonomy, main = "russianautonomym", ylim = ylim)
boxplot(russian$envmattery, main = "russianenvmattery", ylim = ylim)
boxplot(russian$persgrowth, main = "russianpersgrowt", ylim = ylim)
boxplot(russian$purposeinlife, main = "russianpurposeinlife", ylim = ylim)
boxplot(russian$selfacceptance, main = "russianselfacceptanc", ylim = ylim)

par(mfrow = c(1, 5))
max_val <- max(c(russian$conformity, russian$tradition, russian$benevolence, russian$universalism, russian$selfdirection, russian$stimulation, russian$hedonism, russian$achievement, russian$power, russian$security))
ylim <- c(-1, max_val)
boxplot(russian$conformity, main = "russianconformity", ylim = ylim)
boxplot(russian$tradition, main = "russiantradition", ylim = ylim)
boxplot(russian$benevolence, main = "russianbenevolence", ylim = ylim)
boxplot(russian$universalism, main = "russianuniversalism", ylim = ylim)
boxplot(russian$selfdirection, main = "russianselfdirection", ylim = ylim)

boxplot(russian$stimulation, main = "russianstimulation", ylim = ylim)
boxplot(russian$hedonism, main = "russianhedonism", ylim = ylim)
boxplot(russian$achievement, main = "russianhachievement", ylim = ylim)
boxplot(russian$power, main = "russianpower", ylim = ylim)
boxplot(russian$security, main = "russiansecurity", ylim = ylim)

more histograms

for (var in c("autonomy", "positiverelat", "envmattery", "persgrowth", "purposeinlife", "selfacceptance", "pwb")) {
  p <- thes %>%
    ggplot(aes_string(x = var)) +
    geom_histogram(binwidth = 5) +
    facet_grid(rows = vars(sex), cols = vars(nationality)) +
    labs(title = " ", x = var, y = "Count") +
    theme_light()
  print(p)
  
}
## Warning: `aes_string()` was deprecated in ggplot2 3.0.0.
## ℹ Please use tidy evaluation idioms with `aes()`.
## ℹ See also `vignette("ggplot2-in-packages")` for more information.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

for (var in c("conformity", "tradition", "benevolence", "universalism", "selfdirection", "stimulation", "hedonism", "power", "achievement", "security")) {
  p <- thes %>%
    ggplot(aes_string(x = var)) +
    geom_histogram(bins = 20) +
    facet_grid(rows = vars(sex), cols = vars(nationality)) +
    labs(title = " ", x = var, y = "Count") +
    theme_light()
  print(p)

}

сравнение боксплотов для фр и рус

ggplot(gather(thes, key = "variable", value = "value", positiverelat:pwb), 
       aes(x = nationality, y = value, fill = nationality)) +
  geom_boxplot() +
  facet_wrap(~ variable, scales = "free_y") +  # Добавление разбиения по переменным
  scale_x_discrete(labels = c("", ""))+
  theme_minimal() +
  theme(axis.title.x = element_blank(),  
        axis.title.y = element_blank()) 

ggplot(gather(thes, key = "variable", value = "value", conformity:security), 
       aes(x = nationality, y = value, fill = nationality)) +
  geom_boxplot() +
  facet_wrap(~ variable, scales = "free_y") +  
  scale_x_discrete(labels = c("", "")) +
  theme_minimal() +
  theme(axis.title.x = element_blank(),  
        axis.title.y = element_blank()) 

table with descriptive statistics

generate_descriptive_table <- function(data) {
  data <- data[, -(1:5)]
  means <- sapply(data, mean)
  sds <- sapply(data, sd)

  descriptive_table <- data.frame(
    Mean = means,
    SD = sds
    
  )
  
  return(descriptive_table)
}

# Generate descriptive table for 'thes', 'french', and 'russian'
descriptive_table_thes <- generate_descriptive_table(thes)
descriptive_table_french <- generate_descriptive_table(french)
descriptive_table_russian <- generate_descriptive_table(russian)
descriptive_table_frenchmen <- generate_descriptive_table(french_men)
descriptive_table_russianmen <- generate_descriptive_table(russian_men)
descriptive_table_frenchwomen <- generate_descriptive_table(french_women)
descriptive_table_russianwomen <- generate_descriptive_table(russian_women)
descriptive_table_women <- generate_descriptive_table(women)
descriptive_table_men <- generate_descriptive_table(men)

print(descriptive_table_thes)
##                      Mean         SD
## positiverelat   56.293578  7.6995704
## autonomy        51.752294  8.4734399
## envmattery      49.541284  8.7416967
## persgrowth      55.880734  9.1941253
## purposeinlife   58.614679  9.6243405
## selfacceptance  56.642202 11.0147401
## pwb            328.724771 37.9021331
## ei              59.853211 12.6292031
## conformity       3.166972  1.0898569
## tradition        3.095413  1.5032302
## benevolence      4.985321  1.0518789
## universalism     4.962110  1.1389226
## selfdirection    5.344954  0.9842732
## stimulation      4.339174  1.5484671
## hedonism         4.975413  1.2970505
## achievement      4.266055  1.2893237
## power            2.711009  1.5935977
## security         4.460550  1.2879799
print(descriptive_table_french)
##                      Mean         SD
## positiverelat   54.596774  5.2274800
## autonomy        47.322581  4.7105539
## envmattery      44.451613  5.2344302
## persgrowth      49.387097  4.0218862
## purposeinlife   58.258065  9.9208928
## selfacceptance  56.258065 11.4389835
## pwb            310.274194 23.1347934
## ei              59.838710 12.8116156
## conformity       3.454839  1.0895446
## tradition        3.664516  1.3807929
## benevolence      5.267742  0.9560462
## universalism     5.251935  1.0799998
## selfdirection    5.445161  1.0659680
## stimulation      4.714516  1.5557674
## hedonism         5.086129  1.2739114
## achievement      4.310484  1.3832172
## power            2.713710  1.5777960
## security         4.661290  1.4708484
print(descriptive_table_russian)
##                      Mean         SD
## positiverelat   58.531915  9.6927832
## autonomy        57.595745  8.8210706
## envmattery      56.255319  7.8808190
## persgrowth      64.446809  6.7237687
## purposeinlife   59.085106  9.3036454
## selfacceptance  57.148936 10.5294178
## pwb            353.063830 39.9871742
## ei              59.872340 12.5223574
## conformity       2.787234  0.9774895
## tradition        2.344681  1.3284199
## benevolence      4.612766  1.0659996
## universalism     4.579787  1.1118267
## selfdirection    5.212766  0.8581114
## stimulation      3.844043  1.4069311
## hedonism         4.829362  1.3264556
## achievement      4.207447  1.1660030
## power            2.707447  1.6313164
## security         4.195745  0.9486735
print(descriptive_table_frenchmen)
##                    Mean         SD
## positiverelat   55.6400  4.3193364
## autonomy        47.4000  4.6992907
## envmattery      43.8000  5.5226805
## persgrowth      48.1200  3.8656608
## purposeinlife   60.1600  9.7625133
## selfacceptance  58.3200 12.0163777
## pwb            313.4400 20.4818619
## ei              55.7200 12.4616746
## conformity       3.1840  1.0565983
## tradition        3.3200  1.3625956
## benevolence      5.0480  0.9544283
## universalism     5.0028  0.8974246
## selfdirection    5.5200  0.9165151
## stimulation      4.8000  1.3678542
## hedonism         5.2000  0.9954563
## achievement      4.3100  1.1732256
## power            2.6800  1.4974284
## security         4.7840  1.3927431
print(descriptive_table_russianmen)
##                    Mean         SD
## positiverelat   57.2500  9.4081432
## autonomy        61.8500  8.4683622
## envmattery      54.8500  9.0220490
## persgrowth      63.9000  8.1233905
## purposeinlife   57.1500 10.1373463
## selfacceptance  56.7500 11.5752686
## pwb            351.7500 42.9784218
## ei              58.2000 10.3902789
## conformity       2.8200  1.0620734
## tradition        2.4200  1.4095912
## benevolence      4.5500  1.1199154
## universalism     4.3845  1.1816112
## selfdirection    5.2300  0.9342039
## stimulation      3.9000  1.1996622
## hedonism         4.7160  1.5278378
## achievement      4.1000  1.2576084
## power            2.3500  1.8872843
## security         4.0000  0.9267941
print(descriptive_table_frenchwomen)
##                      Mean        SD
## positiverelat   53.891892  5.709562
## autonomy        47.270270  4.782193
## envmattery      44.891892  5.059555
## persgrowth      50.243243  3.946907
## purposeinlife   56.972973  9.951232
## selfacceptance  54.864865 10.977660
## pwb            308.135135 24.810574
## ei              62.621622 12.443989
## conformity       3.637838  1.087135
## tradition        3.897297  1.361879
## benevolence      5.416216  0.940896
## universalism     5.420270  1.169298
## selfdirection    5.394595  1.165701
## stimulation      4.656757  1.686831
## hedonism         5.009189  1.440180
## achievement      4.310811  1.524573
## power            2.736486  1.649859
## security         4.578378  1.534625
print(descriptive_table_russianwomen)
##                      Mean         SD
## positiverelat   59.481481  9.9667539
## autonomy        54.444444  7.8118910
## envmattery      57.296296  6.9104982
## persgrowth      64.851852  5.6001628
## purposeinlife   60.518519  8.5456708
## selfacceptance  57.444444  9.9007899
## pwb            354.037037 38.4292371
## ei              61.111111 13.9568933
## conformity       2.762963  0.9299082
## tradition        2.288889  1.2894046
## benevolence      4.659259  1.0433892
## universalism     4.724444  1.0562573
## selfdirection    5.200000  0.8152395
## stimulation      3.802593  1.5640347
## hedonism         4.913333  1.1789402
## achievement      4.287037  1.1109241
## power            2.972222  1.3907893
## security         4.340741  0.9556549
print(descriptive_table_women)
##                      Mean        SD
## positiverelat   56.250000  8.207623
## autonomy        50.296875  7.141828
## envmattery      50.125000  8.512363
## persgrowth      56.406250  8.644144
## purposeinlife   58.468750  9.478411
## selfacceptance  55.953125 10.534041
## pwb            327.500000 38.513242
## ei              61.984375 13.017073
## conformity       3.268750  1.105380
## tradition        3.218750  1.544974
## benevolence      5.096875  1.047441
## universalism     5.126719  1.166916
## selfdirection    5.312500  1.029640
## stimulation      4.296406  1.678171
## hedonism         4.968750  1.327063
## achievement      4.300781  1.355602
## power            2.835938  1.538667
## security         4.478125  1.317823
print(descriptive_table_men)
##                      Mean         SD
## positiverelat   56.355556  7.0037508
## autonomy        53.822222  9.7847026
## envmattery      48.711111  9.0895100
## persgrowth      55.133333  9.9763356
## purposeinlife   58.822222  9.9322452
## selfacceptance  57.622222 11.7149650
## pwb            330.466667 37.3774361
## ei              56.822222 11.5272186
## conformity       3.022222  1.0627770
## tradition        2.920000  1.4405807
## benevolence      4.826667  1.0493288
## universalism     4.728000  1.0671172
## selfdirection    5.391111  0.9253555
## stimulation      4.400000  1.3588866
## hedonism         4.984889  1.2679325
## achievement      4.216667  1.2020343
## power            2.533333  1.6699211
## security         4.435556  1.2586348

u-test tables

table with u-test results french vs russians

numeric_data_fr <- russian[sapply(russian, is.numeric)]
numeric_data_rus <- french[sapply(french, is.numeric)]
# Create an empty data frame to store the results
test_results <- data.frame(Variable = character(), 
                            Test = character(),
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

# Loop through each column in the Russian and French data frames
for (col in names(numeric_data_rus)) {
  # Perform Shapiro-Wilk test for normality
  shapiro_rus <- shapiro.test(numeric_data_rus[[col]])
  shapiro_fr <- shapiro.test(numeric_data_fr[[col]])
  
  # Check if both distributions are normal
  if (shapiro_rus$p.value > 0.05 & shapiro_fr$p.value > 0.05) {
    # If both distributions are normal, perform t-test
    test_result <- t.test(numeric_data_rus[[col]], numeric_data_fr[[col]])
    test_type <- "t-test"
  } else {
    # If one or both distributions are not normal, perform Wilcoxon rank sum test
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    test_type <- "Wilcoxon rank-sum test"
  }
  
  # Store the results in the data frame
  test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Test = test_type,
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable                   Test     Statistic
## W              age Wilcoxon rank-sum test 1407.00000000
## W1           child Wilcoxon rank-sum test 1819.00000000
## W2           adult Wilcoxon rank-sum test 1836.00000000
## t    positiverelat                 t-test   -2.51937681
## t1        autonomy                 t-test   -7.23991782
## t2      envmattery                 t-test   -8.88890856
## t3      persgrowth                 t-test  -13.61885175
## t4   purposeinlife                 t-test   -0.44661654
## t5  selfacceptance                 t-test   -0.42139627
## t6             pwb                 t-test   -6.55182598
## t7              ei                 t-test   -0.01374833
## t8      conformity                 t-test    3.36008862
## t9       tradition                 t-test    5.05027471
## W3     benevolence Wilcoxon rank-sum test 1917.50000000
## t10   universalism                 t-test    3.16452637
## W4   selfdirection Wilcoxon rank-sum test 1697.00000000
## W5     stimulation Wilcoxon rank-sum test 1927.50000000
## W6        hedonism Wilcoxon rank-sum test 1590.00000000
## t11    achievement                 t-test    0.42139785
## t12          power                 t-test    0.02013239
## t13       security                 t-test    2.00261062
##                           p_value
## W   0.760917326822712625222777660
## W1  0.018871827866932731204707707
## W2  0.001375120237398201215525684
## t   0.014181663399671364314547084
## t1  0.000000000618217525161084634
## t2  0.000000000000228653832081123
## t3  0.000000000000000000002053784
## t4  0.656096788679796594223603279
## t5  0.674344061676286798956425628
## t6  0.000000008566101849043720266
## t7  0.989058071110399561476356212
## t8  0.001090612095625236597010166
## t9  0.000001960664467735884105778
## W3  0.004769742780453545098751089
## t10 0.002071457713388009328620010
## W4  0.141849162633408409917024073
## W5  0.003942249985367435316518936
## W6  0.415735553171173499720225664
## t11 0.674319929112946248039861530
## t12 0.983978935595853476669958582
## t13 0.047807575468086400960476823
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

u-test only, both

numeric_data_fr <- russian[sapply(russian, is.numeric)]
numeric_data_rus <- french[sapply(french, is.numeric)]

test_results <- data.frame(Variable = character(), 
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

for (col in names(numeric_data_rus)) {
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    
  
    test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable Statistic                   p_value
## W              age    1407.0 0.76091732682271262522278
## W1           child    1819.0 0.01887182786693273120471
## W2           adult    1836.0 0.00137512023739820121553
## W3   positiverelat    1111.5 0.03455242911204069178321
## W4        autonomy     425.0 0.00000000025839638949183
## W5      envmattery     322.0 0.00000000000370075517420
## W6      persgrowth      80.0 0.00000000000000003443398
## W7   purposeinlife    1405.5 0.75484660166107520051781
## W8  selfacceptance    1381.0 0.64388326379770788321366
## W9             pwb     511.5 0.00000000734352084766878
## W10             ei    1474.0 0.91953571231547903686021
## W11     conformity    1988.5 0.00112604054899900080100
## W12      tradition    2214.5 0.00000352218617667844349
## W13    benevolence    1917.5 0.00476974278045354509875
## W14   universalism    1952.5 0.00242865988142633004027
## W15  selfdirection    1697.0 0.14184916263340840991702
## W16    stimulation    1927.5 0.00394224998536743531652
## W17       hedonism    1590.0 0.41573555317117349972023
## W18    achievement    1532.0 0.64793831620696029638395
## W19          power    1463.0 0.97310907449596739660080
## W20       security    1767.0 0.05782807333616002409826
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

table with u-test results french men vs russian men

numeric_data_fr <- russian_men[sapply(russian_men, is.numeric)]
numeric_data_rus <- french_men[sapply(french_men, is.numeric)]
test_results <- data.frame(Variable = character(), 
                            Test = character(),
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

for (col in names(numeric_data_rus)) {
  shapiro_rus <- shapiro.test(numeric_data_rus[[col]])
  shapiro_fr <- shapiro.test(numeric_data_fr[[col]])
  
    if (shapiro_rus$p.value > 0.05 & shapiro_fr$p.value > 0.05) {
    test_result <- t.test(numeric_data_rus[[col]], numeric_data_fr[[col]])
    test_type <- "t-test"
  } else {
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    test_type <- "Wilcoxon rank-sum test"
  }
  
  test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Test = test_type,
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable                   Test   Statistic          p_value
## t              age                 t-test  -0.8812136 0.38312329109061
## W            child Wilcoxon rank-sum test 301.0000000 0.21810130141972
## W1           adult Wilcoxon rank-sum test 291.5000000 0.17946527230264
## t1   positiverelat                 t-test  -0.7079454 0.48543123423379
## t2        autonomy                 t-test  -6.8353868 0.00000019356305
## t3      envmattery                 t-test  -4.8044076 0.00004055472391
## t4      persgrowth                 t-test  -7.9933763 0.00000001880539
## t5   purposeinlife                 t-test   1.0061008 0.32039178362566
## t6  selfacceptance                 t-test   0.4445064 0.65898731654626
## t7             pwb                 t-test  -3.6671202 0.00111397163532
## t8              ei                 t-test  -0.7278511 0.47065821703397
## t9      conformity                 t-test   1.1450391 0.25887103036891
## t10      tradition                 t-test   2.1599818 0.03678785471831
## W2     benevolence Wilcoxon rank-sum test 302.5000000 0.23321016693708
## t11   universalism                 t-test   1.9357343 0.06107502543070
## t12  selfdirection                 t-test   1.0434839 0.30291148048278
## t13    stimulation                 t-test   2.3489771 0.02353538138754
## W3        hedonism Wilcoxon rank-sum test 289.0000000 0.37649764133648
## t14    achievement                 t-test   0.5733826 0.56963541226276
## t15          power                 t-test   0.6377079 0.52772478319487
## t16       security                 t-test   2.2581688 0.02922151519819
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

u-test only, men

numeric_data_fr <- russian_men[sapply(russian_men, is.numeric)]
numeric_data_rus <- french_men[sapply(french_men, is.numeric)]

test_results <- data.frame(Variable = character(), 
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

for (col in names(numeric_data_rus)) {
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    
  
    test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable Statistic          p_value
## W              age     206.0 0.31829951075971
## W1           child     301.0 0.21810130141972
## W2           adult     291.5 0.17946527230264
## W3   positiverelat     226.5 0.59832092707405
## W4        autonomy      25.5 0.00000029227946
## W5      envmattery      85.5 0.00017415881637
## W6      persgrowth      15.0 0.00000008133706
## W7   purposeinlife     302.5 0.23447906509972
## W8  selfacceptance     279.5 0.50722062849332
## W9             pwb     106.5 0.00108612613947
## W10             ei     221.5 0.52209881857258
## W11     conformity     311.0 0.16548455409799
## W12      tradition     341.5 0.03724657321170
## W13    benevolence     302.5 0.23321016693708
## W14   universalism     339.0 0.04289960976727
## W15  selfdirection     295.0 0.30787976121890
## W16    stimulation     335.0 0.05254672867911
## W17       hedonism     289.0 0.37649764133648
## W18    achievement     279.0 0.51397237369999
## W19          power     292.0 0.34168360400989
## W20       security     334.0 0.05600564256254
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

table with u-test results french women vs russian women

numeric_data_fr <- russian_women[sapply(russian_women, is.numeric)]
numeric_data_rus <- french_women[sapply(french_women, is.numeric)]
test_results <- data.frame(Variable = character(), 
                            Test = character(),
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

for (col in names(numeric_data_rus)) {
  shapiro_rus <- shapiro.test(numeric_data_rus[[col]])
  shapiro_fr <- shapiro.test(numeric_data_fr[[col]])
  
  if (shapiro_rus$p.value > 0.05 & shapiro_fr$p.value > 0.05) {
    test_result <- t.test(numeric_data_rus[[col]], numeric_data_fr[[col]])
    test_type <- "t-test"
  } else {
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    test_type <- "Wilcoxon rank-sum test"
  }
  
  test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Test = test_type,
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable                   Test    Statistic                p_value
## W              age Wilcoxon rank-sum test 515.50000000 0.83143687975764191123
## W1           child Wilcoxon rank-sum test 637.50000000 0.04738706015072723099
## W2           adult Wilcoxon rank-sum test 663.50000000 0.00269920160280503682
## t    positiverelat                 t-test  -2.61751504 0.01260821199794503154
## t1        autonomy                 t-test  -4.22867091 0.00013269000398823542
## t2      envmattery                 t-test  -7.90784213 0.00000000044804961332
## t3      persgrowth                 t-test -11.61251454 0.00000000000000531399
## t4   purposeinlife                 t-test  -1.52842775 0.13163693453807800982
## W3  selfacceptance Wilcoxon rank-sum test 417.00000000 0.26447304053143699809
## t5             pwb                 t-test  -5.43480783 0.00000266262816130139
## t6              ei                 t-test   0.44737780 0.65645249885344125751
## t7      conformity                 t-test   3.45907835 0.00100018763035765847
## t8       tradition                 t-test   4.81240555 0.00001107387434929198
## W4     benevolence Wilcoxon rank-sum test 689.00000000 0.00999873610256550582
## W5    universalism Wilcoxon rank-sum test 687.00000000 0.01088984146276868528
## W6   selfdirection Wilcoxon rank-sum test 583.00000000 0.25795849910366508428
## t9     stimulation                 t-test   2.08703729 0.04125492106869686759
## W7        hedonism Wilcoxon rank-sum test 527.50000000 0.70754536665384359928
## t10    achievement                 t-test   0.07216479 0.94270306173292839347
## t11          power                 t-test  -0.61862658 0.53847947357061609619
## t12       security                 t-test   0.76114427 0.44951807371780816425
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

u-test only, women

numeric_data_fr <- russian_women[sapply(russian_women, is.numeric)]
numeric_data_rus <- french_women[sapply(french_women, is.numeric)]

test_results <- data.frame(Variable = character(), 
                            Statistic = numeric(),
                            p_value = numeric(),
                            stringsAsFactors = FALSE)

for (col in names(numeric_data_rus)) {
    test_result <- wilcox.test(numeric_data_rus[[col]], numeric_data_fr[[col]], exact = FALSE)
    
  
    test_results <- rbind(test_results, 
                        data.frame(Variable = col, 
                                   Statistic = test_result$statistic,
                                   p_value = test_result$p.value))
}

print(test_results)
##           Variable Statistic             p_value
## W              age     515.5 0.83143687975764191
## W1           child     637.5 0.04738706015072723
## W2           adult     663.5 0.00269920160280504
## W3   positiverelat     338.0 0.02844723981303018
## W4        autonomy     210.0 0.00008196269587791
## W5      envmattery      68.0 0.00000000447555841
## W6      persgrowth      19.5 0.00000000006781702
## W7   purposeinlife     396.5 0.16309444926110833
## W8  selfacceptance     417.0 0.26447304053143700
## W9             pwb     151.0 0.00000222911842201
## W10             ei     537.0 0.61465228518332737
## W11     conformity     730.5 0.00167456882710692
## W12      tradition     816.5 0.00001645773506943
## W13    benevolence     689.0 0.00999873610256551
## W14   universalism     687.0 0.01088984146276869
## W15  selfdirection     583.0 0.25795849910366508
## W16    stimulation     654.5 0.03530430363219323
## W17       hedonism     527.5 0.70754536665384360
## W18    achievement     511.5 0.87554045605105346
## W19          power     455.0 0.54890268657708696
## W20       security     555.5 0.44942598435587572
test_results$p_value <- format(test_results$p_value, scientific = FALSE)

#regresions ## both

thes %>% 
  pivot_longer(cols = c("autonomy", "positiverelat", "envmattery", "persgrowth", "purposeinlife", "selfacceptance"),                names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

 ggplot(thes, aes(x = ei, y = pwb , col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        labs(title = "Well-being", x = "EnvIdentity", y = "Well-being" ) +
        theme_light()
## `geom_smooth()` using formula = 'y ~ x'

thes %>% 
  pivot_longer(cols = c("conformity", "tradition", "benevolence", "universalism", "selfdirection", "stimulation", "hedonism", "power", "achievement", "security"),
               names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

men

men %>% 
  pivot_longer(cols = c("autonomy", "positiverelat", "envmattery", "persgrowth", "purposeinlife", "selfacceptance"),                names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

 ggplot(men, aes(x = ei, y = pwb , col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        labs(title = "Well-being", x = "EnvIdentity", y = "Well-being" ) +
        theme_light()
## `geom_smooth()` using formula = 'y ~ x'

men %>% 
  pivot_longer(cols = c("conformity", "tradition", "benevolence", "universalism", "selfdirection", "stimulation", "hedonism", "power", "achievement", "security"),
               names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

women

women %>% 
  pivot_longer(cols = c("autonomy", "positiverelat", "envmattery", "persgrowth", "purposeinlife", "selfacceptance"),                names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

 ggplot(women, aes(x = ei, y = pwb , col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        labs(title = "Well-being", x = "EnvIdentity", y = "Well-being" ) +
        theme_light()
## `geom_smooth()` using formula = 'y ~ x'

women %>% 
  pivot_longer(cols = c("conformity", "tradition", "benevolence", "universalism", "selfdirection", "stimulation", "hedonism", "power", "achievement", "security"),
               names_to = "Variable", 
               values_to = "Value") %>%
  ggplot(., aes(x = ei, y = Value, col = nationality)) +
        geom_point() +
        geom_smooth(method = 'lm') +
        facet_grid(cols = vars(Variable)) +
        labs(title = "", x = "Well-Being", y = "EnvIdentity") +
        theme_light() 
## `geom_smooth()` using formula = 'y ~ x'

correlations

correlation map with p-value for both samples

numeric_data <- thes[sapply(thes, is.numeric)]

correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 0.17 -0.07 -0.27 -0.2 -0.15 -0.07 -0.15 -0.22 0.25 0.12 0.14 0.14 0.11 0.07 0.14 0.06 0.12 0.01 0.13
adult 0.17 1 -0.09 -0.21 -0.33 -0.21 0 -0.12 -0.23 0.12 0.03 0.18 0.06 0.15 0.15 0.15 -0.01 -0.07 -0.04 0.02
positiverelat -0.07 -0.09 1 0.32 0.33 0.39 0.48 0.48 0.71 0.06 0.07 -0.12 0.21 0.22 0.16 0.06 0.02 0.14 -0.17 0.03
autonomy -0.27* -0.21* 0.32* 1 0.57 0.59 0.12 0.19 0.65 -0.1 -0.05 -0.18 -0.13 -0.1 -0.04 -0.22 -0.22 0.01 -0.04 -0.17
envmattery -0.2* -0.33* 0.33* 0.57* 1 0.71 0.26 0.24 0.73 -0.07 -0.09 -0.22 -0.15 -0.26 -0.18 -0.36 -0.22 0.05 -0.01 -0.21
persgrowth -0.15 -0.21* 0.39* 0.59* 0.71* 1 0.16 0.1 0.69 0.07 -0.2 -0.33 -0.14 -0.14 -0.01 -0.18 -0.21 0.03 -0.03 -0.16
purposeinlife -0.07 0 0.48* 0.12 0.26* 0.16 1 0.71 0.68 0.02 0.09 0.1 0.2 0.2 0.37 0.21 -0.02 0.33 0.03 0.19
selfacceptance -0.15 -0.12 0.48* 0.19* 0.24* 0.1 0.71* 1 0.69 0.09 0.22 0.18 0.34 0.33 0.36 0.2 0.05 0.22 0.04 0.25
pwb -0.22* -0.23* 0.71* 0.65* 0.73* 0.69* 0.68* 0.69* 1 0.02 0.02 -0.12 0.1 0.08 0.18 -0.05 -0.14 0.2 -0.03 0
ei 0.25* 0.12 0.06 -0.1 -0.07 0.07 0.02 0.09 0.02 1 0.26 0.26 0.35 0.44 0.21 0.1 0.11 0.11 -0.06 0.25
conformity 0.12 0.03 0.07 -0.05 -0.09 -0.2* 0.09 0.22* 0.02 0.26* 1 0.71 0.61 0.48 0.46 0.28 0.2 0.55 0.39 0.59
tradition 0.14 0.18 -0.12 -0.18 -0.22* -0.33* 0.1 0.18 -0.12 0.26* 0.71* 1 0.52 0.48 0.4 0.33 0.1 0.36 0.33 0.44
benevolence 0.14 0.06 0.21* -0.13 -0.15 -0.14 0.2* 0.34* 0.1 0.35* 0.61* 0.52* 1 0.62 0.57 0.39 0.27 0.49 0.11 0.42
universalism 0.11 0.15 0.22* -0.1 -0.26* -0.14 0.2* 0.33* 0.08 0.44* 0.48* 0.48* 0.62* 1 0.59 0.37 0.2 0.25 0.01 0.42
selfdirection 0.07 0.15 0.16 -0.04 -0.18 -0.01 0.37* 0.36* 0.18 0.21* 0.46* 0.4* 0.57* 0.59* 1 0.48 0.35 0.5 0.21 0.47
stimulation 0.14 0.15 0.06 -0.22* -0.36* -0.18 0.21* 0.2* -0.05 0.1 0.28* 0.33* 0.39* 0.37* 0.48* 1 0.41 0.44 0.4 0.43
hedonism 0.06 -0.01 0.02 -0.22* -0.22* -0.21* -0.02 0.05 -0.14 0.11 0.2* 0.1 0.27* 0.2* 0.35* 0.41* 1 0.42 0.32 0.43
achievement 0.12 -0.07 0.14 0.01 0.05 0.03 0.33* 0.22* 0.2* 0.11 0.55* 0.36* 0.49* 0.25* 0.5* 0.44* 0.42* 1 0.56 0.45
power 0.01 -0.04 -0.17 -0.04 -0.01 -0.03 0.03 0.04 -0.03 -0.06 0.39* 0.33* 0.11 0.01 0.21* 0.4* 0.32* 0.56* 1 0.41
security 0.13 0.02 0.03 -0.17 -0.21* -0.16 0.19* 0.25* 0 0.25* 0.59* 0.44* 0.42* 0.42* 0.47* 0.43* 0.43* 0.45* 0.41* 1

correlation matrix for russian

numeric_data <- russian[sapply(russian, is.numeric)]

correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 0.53 0 -0.12 -0.24 0.03 -0.13 -0.31 -0.18 0.18 0.03 0.06 -0.04 0.15 0.07 0.17 0.12 0.11 0.04 0.26
adult 0.53* 1 -0.13 -0.18 -0.32 -0.27 -0.16 -0.3 -0.3 0.01 -0.23 -0.07 -0.31 -0.08 0.09 0 0.18 -0.01 -0.01 -0.02
positiverelat 0 -0.13 1 0.3 0.43 0.52 0.4 0.46 0.7 0.09 0.25 -0.06 0.45 0.31 0.1 0.1 0.1 0.26 -0.25 0.09
autonomy -0.12 -0.18 0.3* 1 0.39 0.36 0.35 0.47 0.64 -0.22 0.14 0.08 0.18 0.16 0.14 0.03 -0.34 0.02 -0.22 -0.09
envmattery -0.24 -0.32* 0.43* 0.39* 1 0.53 0.65 0.65 0.8 -0.04 0.3 0.15 0.26 0.11 0.11 -0.01 -0.17 0.25 0.04 0.03
persgrowth 0.03 -0.27 0.52* 0.36* 0.53* 1 0.57 0.58 0.76 0.23 0.18 0.12 0.4 0.42 0.39 0.25 -0.19 0.26 -0.04 0.16
purposeinlife -0.13 -0.16 0.4* 0.35* 0.65* 0.57* 1 0.6 0.79 0.12 0.06 0.1 0.28 0.16 0.28 0.07 -0.21 0.38 0.03 0.07
selfacceptance -0.31* -0.3* 0.46* 0.47* 0.65* 0.58* 0.6* 1 0.84 0.09 0.34 0.16 0.42 0.4 0.3 0.12 -0.25 0.16 0.03 0.16
pwb -0.18 -0.3* 0.7* 0.64* 0.8* 0.76* 0.79* 0.84* 1 0.05 0.29 0.12 0.45 0.34 0.29 0.12 -0.23 0.29 -0.09 0.09
ei 0.18 0.01 0.09 -0.22 -0.04 0.23 0.12 0.09 0.05 1 0.28 0.27 0.24 0.43 0.26 0.09 0.17 0.18 0.08 0.46
conformity 0.03 -0.23 0.25 0.14 0.3* 0.18 0.06 0.34* 0.29 0.28 1 0.59 0.58 0.47 0.39 0.11 0 0.42 0.15 0.57
tradition 0.06 -0.07 -0.06 0.08 0.15 0.12 0.1 0.16 0.12 0.27 0.59* 1 0.37 0.44 0.41 0.23 -0.11 0.3 0.32 0.33
benevolence -0.04 -0.31* 0.45* 0.18 0.26 0.4* 0.28 0.42* 0.45* 0.24 0.58* 0.37* 1 0.58 0.52 0.16 0.22 0.47 0.07 0.34
universalism 0.15 -0.08 0.31* 0.16 0.11 0.42* 0.16 0.4* 0.34* 0.43* 0.47* 0.44* 0.58* 1 0.57 0.09 -0.02 0.21 -0.02 0.59
selfdirection 0.07 0.09 0.1 0.14 0.11 0.39* 0.28 0.3* 0.29 0.26 0.39* 0.41* 0.52* 0.57* 1 0.32 0.11 0.44 0.16 0.43
stimulation 0.17 0 0.1 0.03 -0.01 0.25 0.07 0.12 0.12 0.09 0.11 0.23 0.16 0.09 0.32* 1 0.23 0.34 0.49 0.12
hedonism 0.12 0.18 0.1 -0.34* -0.17 -0.19 -0.21 -0.25 -0.23 0.17 0 -0.11 0.22 -0.02 0.11 0.23 1 0.24 0.18 0.1
achievement 0.11 -0.01 0.26 0.02 0.25 0.26 0.38* 0.16 0.29* 0.18 0.42* 0.3* 0.47* 0.21 0.44* 0.34* 0.24 1 0.47 0.26
power 0.04 -0.01 -0.25 -0.22 0.04 -0.04 0.03 0.03 -0.09 0.08 0.15 0.32* 0.07 -0.02 0.16 0.49* 0.18 0.47* 1 0.19
security 0.26 -0.02 0.09 -0.09 0.03 0.16 0.07 0.16 0.09 0.46* 0.57* 0.33* 0.34* 0.59* 0.43* 0.12 0.1 0.26 0.19 1

correlation matrix for french

numeric_data <- french[sapply(french, is.numeric)]

correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 -0.08 -0.02 -0.28 0.1 0.09 -0.02 -0.04 -0.05 0.3 0.08 0.05 0.17 -0.02 0.04 0.03 -0.01 0.12 -0.01 0.02
adult -0.08 1 0.06 0.02 -0.17 0.24 0.11 -0.02 0.06 0.19 0.04 0.15 0.15 0.18 0.14 0.13 -0.16 -0.11 -0.06 -0.03
positiverelat -0.02 0.06 1 -0.03 -0.2 -0.15 0.67 0.6 0.73 0.03 0.06 0.07 0.12 0.37 0.35 0.21 -0.02 0.07 -0.08 0.09
autonomy -0.28* 0.02 -0.03 1 0.07 -0.09 -0.18 -0.1 0.07 -0.03 0.23 0.18 -0.08 0.02 -0.06 -0.2 -0.03 0.08 0.21 -0.09
envmattery 0.1 -0.17 -0.2 0.07 1 0.07 -0.04 -0.08 0.15 -0.16 0.03 0.07 -0.14 -0.31 -0.39 -0.5 -0.26 -0.03 -0.06 -0.26
persgrowth 0.09 0.24 -0.15 -0.09 0.07 1 -0.17 -0.4 -0.13 -0.01 0 -0.02 -0.07 -0.1 -0.07 -0.07 -0.28 -0.04 -0.07 -0.18
purposeinlife -0.02 0.11 0.67* -0.18 -0.04 -0.17 1 0.78 0.89 -0.05 0.13 0.15 0.19 0.27 0.44 0.33 0.13 0.31 0.02 0.27
selfacceptance -0.04 -0.02 0.6* -0.1 -0.08 -0.4* 0.78* 1 0.86 0.09 0.19 0.26 0.35 0.32 0.41 0.29 0.28 0.26 0.05 0.32
pwb -0.05 0.06 0.73* 0.07 0.15 -0.13 0.89* 0.86* 1 -0.01 0.22 0.26 0.22 0.28 0.36 0.17 0.08 0.28 0.03 0.18
ei 0.3* 0.19 0.03 -0.03 -0.16 -0.01 -0.05 0.09 -0.01 1 0.27 0.29 0.47 0.49 0.18 0.12 0.07 0.07 -0.17 0.16
conformity 0.08 0.04 0.06 0.23 0.03 0 0.13 0.19 0.22 0.27* 1 0.72 0.56 0.4 0.48 0.29 0.31 0.64 0.59 0.58
tradition 0.05 0.15 0.07 0.18 0.07 -0.02 0.15 0.26* 0.26* 0.29* 0.72* 1 0.52 0.38 0.37 0.24 0.19 0.44 0.4 0.46
benevolence 0.17 0.15 0.12 -0.08 -0.14 -0.07 0.19 0.35* 0.22 0.47* 0.56* 0.52* 1 0.59 0.61 0.46 0.29 0.54 0.16 0.45
universalism -0.02 0.18 0.37* 0.02 -0.31* -0.1 0.27* 0.32* 0.28* 0.49* 0.4* 0.38* 0.59* 1 0.61 0.47 0.34 0.28 0.03 0.31
selfdirection 0.04 0.14 0.35* -0.06 -0.39* -0.07 0.44* 0.41* 0.36* 0.18 0.48* 0.37* 0.61* 0.61* 1 0.55 0.48 0.53 0.25 0.47
stimulation 0.03 0.13 0.21 -0.2 -0.5* -0.07 0.33* 0.29* 0.17 0.12 0.29* 0.24 0.46* 0.47* 0.55* 1 0.52 0.52 0.37 0.53
hedonism -0.01 -0.16 -0.02 -0.03 -0.26* -0.28* 0.13 0.28* 0.08 0.07 0.31* 0.19 0.29* 0.34* 0.48* 0.52* 1 0.54 0.44 0.6
achievement 0.12 -0.11 0.07 0.08 -0.03 -0.04 0.31* 0.26* 0.28* 0.07 0.64* 0.44* 0.54* 0.28* 0.53* 0.52* 0.54* 1 0.62 0.54
power -0.01 -0.06 -0.08 0.21 -0.06 -0.07 0.02 0.05 0.03 -0.17 0.59* 0.4* 0.16 0.03 0.25* 0.37* 0.44* 0.62* 1 0.55
security 0.02 -0.03 0.09 -0.09 -0.26* -0.18 0.27* 0.32* 0.18 0.16 0.58* 0.46* 0.45* 0.31* 0.47* 0.53* 0.6* 0.54* 0.55* 1

correlation matrix for french women and men separately

numeric_data <- french_men[sapply(french_men, is.numeric)]

correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 -0.22 0.15 -0.35 0.3 -0.09 0.21 0.05 0.15 0.39 0.19 0.16 0.17 -0.02 -0.16 -0.11 -0.08 0.32 -0.05 0.21
adult -0.22 1 0.2 -0.01 -0.28 0.25 0.1 -0.06 0.02 0.16 -0.15 0.01 -0.11 -0.03 -0.01 0 -0.39 -0.37 -0.28 -0.03
positiverelat 0.15 0.2 1 -0.09 -0.16 -0.42 0.77 0.63 0.81 0.06 -0.01 0.19 0.28 0.44 0.25 0.34 -0.04 0.09 -0.21 0.21
autonomy -0.35 -0.01 -0.09 1 -0.21 -0.08 -0.36 -0.11 -0.1 -0.16 0.13 -0.09 -0.3 -0.06 -0.31 -0.16 -0.18 -0.17 0.18 -0.21
envmattery 0.3 -0.28 -0.16 -0.21 1 0.53 -0.13 -0.35 0.02 -0.23 0.33 0.19 -0.08 -0.32 -0.28 -0.38 -0.29 0.33 0.13 -0.16
persgrowth -0.09 0.25 -0.42* -0.08 0.53* 1 -0.39 -0.62 -0.32 -0.13 0.01 -0.12 -0.29 -0.53 -0.4 -0.34 -0.29 -0.09 -0.2 -0.21
purposeinlife 0.21 0.1 0.77* -0.36 -0.13 -0.39 1 0.73 0.88 0.19 0.03 0.18 0.26 0.25 0.35 0.48 -0.08 0.31 -0.08 0.36
selfacceptance 0.05 -0.06 0.63* -0.11 -0.35 -0.62* 0.73* 1 0.83 0.3 0.09 0.33 0.44 0.5 0.48 0.36 0.14 0.17 -0.14 0.36
pwb 0.15 0.02 0.81* -0.1 0.02 -0.32 0.88* 0.83* 1 0.15 0.19 0.33 0.3 0.3 0.28 0.31 -0.13 0.3 -0.13 0.3
ei 0.39 0.16 0.06 -0.16 -0.23 -0.13 0.19 0.3 0.15 1 0.31 0.38 0.66 0.5 0.39 0.03 -0.03 0.19 -0.29 0.49
conformity 0.19 -0.15 -0.01 0.13 0.33 0.01 0.03 0.09 0.19 0.31 1 0.68 0.39 0.23 0.15 -0.21 -0.14 0.55 0.46 0.44
tradition 0.16 0.01 0.19 -0.09 0.19 -0.12 0.18 0.33 0.33 0.38 0.68* 1 0.44 0.37 0.33 -0.08 -0.16 0.19 0.11 0.44
benevolence 0.17 -0.11 0.28 -0.3 -0.08 -0.29 0.26 0.44* 0.3 0.66* 0.39 0.44* 1 0.73 0.65 0.25 0.34 0.47 -0.19 0.61
universalism -0.02 -0.03 0.44* -0.06 -0.32 -0.53* 0.25 0.5* 0.3 0.5* 0.23 0.37 0.73* 1 0.61 0.26 0.38 0.05 -0.13 0.39
selfdirection -0.16 -0.01 0.25 -0.31 -0.28 -0.4* 0.35 0.48* 0.28 0.39 0.15 0.33 0.65* 0.61* 1 0.29 0.43 0.31 -0.06 0.5
stimulation -0.11 0 0.34 -0.16 -0.38 -0.34 0.48* 0.36 0.31 0.03 -0.21 -0.08 0.25 0.26 0.29 1 0.33 0.16 0.11 0.4
hedonism -0.08 -0.39 -0.04 -0.18 -0.29 -0.29 -0.08 0.14 -0.13 -0.03 -0.14 -0.16 0.34 0.38 0.43* 0.33 1 0.1 0.08 0.33
achievement 0.32 -0.37 0.09 -0.17 0.33 -0.09 0.31 0.17 0.3 0.19 0.55* 0.19 0.47* 0.05 0.31 0.16 0.1 1 0.39 0.4
power -0.05 -0.28 -0.21 0.18 0.13 -0.2 -0.08 -0.14 -0.13 -0.29 0.46* 0.11 -0.19 -0.13 -0.06 0.11 0.08 0.39 1 0.15
security 0.21 -0.03 0.21 -0.21 -0.16 -0.21 0.36 0.36 0.3 0.49* 0.44* 0.44* 0.61* 0.39 0.5* 0.4* 0.33 0.4* 0.15 1
numeric_data <- french_women[sapply(french_women, is.numeric)]
correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 -0.01 -0.09 -0.23 -0.04 0.19 -0.16 -0.09 -0.15 0.26 0 -0.02 0.16 -0.03 0.14 0.1 0.02 0.02 0.01 -0.09
adult -0.01 1 0.01 0.04 -0.1 0.22 0.13 0.03 0.09 0.18 0.14 0.23 0.31 0.28 0.22 0.21 -0.05 0.02 0.06 -0.03
positiverelat -0.09 0.01 1 0 -0.2 0.04 0.61 0.58 0.7 0.09 0.16 0.07 0.09 0.4 0.38 0.16 -0.03 0.06 -0.02 0.03
autonomy -0.23 0.04 0 1 0.26 -0.09 -0.07 -0.1 0.16 0.07 0.31 0.37 0.07 0.06 0.06 -0.22 0.04 0.21 0.22 -0.02
envmattery -0.04 -0.1 -0.2 0.26 1 -0.3 0.06 0.17 0.26 -0.18 -0.21 -0.05 -0.22 -0.35 -0.46 -0.58 -0.25 -0.24 -0.2 -0.33
persgrowth 0.19 0.22 0.04 -0.09 -0.3 1 0.03 -0.21 0.01 -0.05 -0.1 -0.05 -0.01 0.02 0.12 0.08 -0.26 -0.02 -0.01 -0.15
purposeinlife -0.16 0.13 0.61* -0.07 0.06 0.03 1 0.81 0.9 -0.14 0.26 0.19 0.2 0.35 0.49 0.25 0.21 0.31 0.09 0.21
selfacceptance -0.09 0.03 0.58* -0.1 0.17 -0.21 0.81* 1 0.88 0.03 0.32 0.28 0.34 0.29 0.38 0.24 0.34 0.32 0.19 0.29
pwb -0.15 0.09 0.7* 0.16 0.26 0.01 0.9* 0.88* 1 -0.05 0.28 0.27 0.22 0.31 0.39 0.1 0.14 0.27 0.11 0.12
ei 0.26 0.18 0.09 0.07 -0.18 -0.05 -0.14 0.03 -0.05 1 0.17 0.17 0.29 0.45 0.11 0.19 0.15 0.01 -0.11 0
conformity 0 0.14 0.16 0.31 -0.21 -0.1 0.26 0.32 0.28 0.17 1 0.73 0.65 0.46 0.68 0.57 0.55 0.71 0.68 0.71
tradition -0.02 0.23 0.07 0.37* -0.05 -0.05 0.19 0.28 0.27 0.17 0.73* 1 0.54 0.35 0.43 0.44 0.39 0.59 0.59 0.51
benevolence 0.16 0.31 0.09 0.07 -0.22 -0.01 0.2 0.34* 0.22 0.29 0.65* 0.54* 1 0.5 0.64 0.61 0.3 0.59 0.38 0.38
universalism -0.03 0.28 0.4* 0.06 -0.35* 0.02 0.35* 0.29 0.31 0.45* 0.46* 0.35* 0.5* 1 0.64 0.59 0.36 0.38 0.1 0.31
selfdirection 0.14 0.22 0.38* 0.06 -0.46* 0.12 0.49* 0.38* 0.39* 0.11 0.68* 0.43* 0.64* 0.64* 1 0.65 0.5 0.62 0.4 0.45
stimulation 0.1 0.21 0.16 -0.22 -0.58* 0.08 0.25 0.24 0.1 0.19 0.57* 0.44* 0.61* 0.59* 0.65* 1 0.59 0.66 0.5 0.59
hedonism 0.02 -0.05 -0.03 0.04 -0.25 -0.26 0.21 0.34* 0.14 0.15 0.55* 0.39* 0.3 0.36* 0.5* 0.59* 1 0.7 0.6 0.71
achievement 0.02 0.02 0.06 0.21 -0.24 -0.02 0.31 0.32 0.27 0.01 0.71* 0.59* 0.59* 0.38* 0.62* 0.66* 0.7* 1 0.74 0.61
power 0.01 0.06 -0.02 0.22 -0.2 -0.01 0.09 0.19 0.11 -0.11 0.68* 0.59* 0.38* 0.1 0.4* 0.5* 0.6* 0.74* 1 0.78
security -0.09 -0.03 0.03 -0.02 -0.33* -0.15 0.21 0.29 0.12 0 0.71* 0.51* 0.38* 0.31 0.45* 0.59* 0.71* 0.61* 0.78* 1

correlation matrix for russian women and men separately

numeric_data <- russian_women[sapply(russian_women, is.numeric)]
correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 0.51 0.04 -0.02 -0.21 0.33 -0.14 -0.16 -0.05 0.16 -0.08 -0.1 -0.15 0.22 0.01 0.02 0.07 -0.19 -0.17 0.31
adult 0.51* 1 -0.03 -0.22 -0.24 0.04 -0.18 -0.03 -0.14 0.11 -0.25 -0.04 -0.16 0.24 0.15 -0.14 0.01 -0.32 -0.28 -0.03
positiverelat 0.04 -0.03 1 0.61 0.38 0.63 0.5 0.64 0.82 0.09 0.11 -0.14 0.37 0.21 0.05 0.26 -0.01 0.22 -0.06 0.04
autonomy -0.02 -0.22 0.61* 1 0.56 0.52 0.55 0.57 0.81 -0.14 -0.08 -0.15 0.33 0.05 0.07 0.3 -0.17 0.15 0.03 -0.18
envmattery -0.21 -0.24 0.38 0.56* 1 0.38 0.62 0.6 0.74 -0.15 0.32 0.11 0.35 -0.21 0.17 0.2 -0.01 0.49 0.3 0
persgrowth 0.33 0.04 0.63* 0.52* 0.38 1 0.39 0.49 0.7 0.08 -0.08 -0.17 0.19 0.17 0.25 0.38 -0.11 0.28 -0.06 -0.02
purposeinlife -0.14 -0.18 0.5* 0.55* 0.62* 0.39* 1 0.57 0.78 0 0.02 -0.03 0.34 -0.08 0.12 0.16 -0.19 0.44 -0.08 -0.07
selfacceptance -0.16 -0.03 0.64* 0.57* 0.6* 0.49* 0.57* 1 0.84 0.09 0.23 0.11 0.44 0.3 0.24 0.35 0.08 0.31 0.09 0.1
pwb -0.05 -0.14 0.82* 0.81* 0.74* 0.7* 0.78* 0.84* 1 0 0.12 -0.05 0.44 0.11 0.18 0.35 -0.08 0.39 0.04 -0.02
ei 0.16 0.11 0.09 -0.14 -0.15 0.08 0 0.09 0 1 0.33 0.22 0.24 0.45 0.17 -0.06 0.25 0.05 -0.22 0.46
conformity -0.08 -0.25 0.11 -0.08 0.32 -0.08 0.02 0.23 0.12 0.33 1 0.59 0.5 0.18 0.2 0.16 0.32 0.54 0.41 0.5
tradition -0.1 -0.04 -0.14 -0.15 0.11 -0.17 -0.03 0.11 -0.05 0.22 0.59* 1 0.35 0.28 0.26 0.32 0.15 0.38 0.39 0.18
benevolence -0.15 -0.16 0.37 0.33 0.35 0.19 0.34 0.44* 0.44* 0.24 0.5* 0.35 1 0.45 0.48 0.16 0.46 0.56 0.28 0.24
universalism 0.22 0.24 0.21 0.05 -0.21 0.17 -0.08 0.3 0.11 0.45* 0.18 0.28 0.45* 1 0.46 0.07 0.19 -0.03 -0.18 0.42
selfdirection 0.01 0.15 0.05 0.07 0.17 0.25 0.12 0.24 0.18 0.17 0.2 0.26 0.48* 0.46* 1 0.29 0.22 0.3 -0.04 0.16
stimulation 0.02 -0.14 0.26 0.3 0.2 0.38* 0.16 0.35 0.35 -0.06 0.16 0.32 0.16 0.07 0.29 1 0.08 0.35 0.48 -0.02
hedonism 0.07 0.01 -0.01 -0.17 -0.01 -0.11 -0.19 0.08 -0.08 0.25 0.32 0.15 0.46* 0.19 0.22 0.08 1 0.2 0.38 0.25
achievement -0.19 -0.32 0.22 0.15 0.49* 0.28 0.44* 0.31 0.39* 0.05 0.54* 0.38 0.56* -0.03 0.3 0.35 0.2 1 0.61 0.05
power -0.17 -0.28 -0.06 0.03 0.3 -0.06 -0.08 0.09 0.04 -0.22 0.41* 0.39* 0.28 -0.18 -0.04 0.48* 0.38 0.61* 1 -0.07
security 0.31 -0.03 0.04 -0.18 0 -0.02 -0.07 0.1 -0.02 0.46* 0.5* 0.18 0.24 0.42* 0.16 -0.02 0.25 0.05 -0.07 1
numeric_data <- russian_men[sapply(russian_men, is.numeric)]
correlation_matrix <- cor(numeric_data) 
p_values <- matrix(NA, nrow = ncol(numeric_data), ncol = ncol(numeric_data))
rownames(p_values) <- colnames(p_values) <- colnames(numeric_data)
for (i in 1:(ncol(numeric_data) - 1)) {
  for (j in (i + 1):ncol(numeric_data)) {
    cor_test_result <- cor.test(numeric_data[[i]], numeric_data[[j]])
    p_values[i, j] <- p_values[j, i] <- cor_test_result$p.value
  }
}

significant_level <- 0.05
significant <- ifelse(p_values < significant_level, "*", "")

cor_matrix_with_significance <- round(correlation_matrix, 2)
cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)] <- paste(cor_matrix_with_significance[lower.tri(cor_matrix_with_significance)], significant[lower.tri(significant)], sep = "")

cor_thes <- as.data.frame(cor_matrix_with_significance)
cor_thes <- cor_thes[-1, -1]
kable(cor_thes, format = "html") %>%
  kable_styling()
child adult positiverelat autonomy envmattery persgrowth purposeinlife selfacceptance pwb ei conformity tradition benevolence universalism selfdirection stimulation hedonism achievement power security
child 1 0.55 -0.05 -0.31 -0.25 -0.21 -0.1 -0.46 -0.31 0.25 0.13 0.23 0.08 0.1 0.12 0.38 0.18 0.43 0.24 0.25
adult 0.55* 1 -0.22 -0.29 -0.37 -0.47 -0.12 -0.53 -0.44 -0.1 -0.22 -0.11 -0.44 -0.34 0.03 0.18 0.33 0.28 0.22 0.03
positiverelat -0.05 -0.22 1 0.11 0.47 0.43 0.26 0.25 0.55 0.06 0.44 0.05 0.57 0.42 0.17 -0.18 0.22 0.3 -0.52 0.12
autonomy -0.31 -0.29 0.11 1 0.46 0.35 0.41 0.48 0.61 -0.28 0.4 0.31 0.1 0.47 0.23 -0.45 -0.51 -0.03 -0.34 0.2
envmattery -0.25 -0.37 0.47* 0.46* 1 0.63 0.67 0.7 0.87 0.05 0.3 0.2 0.17 0.36 0.07 -0.29 -0.33 0 -0.22 0.01
persgrowth -0.21 -0.47* 0.43 0.35 0.63* 1 0.71 0.67 0.83 0.43 0.41 0.38 0.6 0.62 0.52 0.12 -0.27 0.24 -0.04 0.33
purposeinlife -0.1 -0.12 0.26 0.41 0.67* 0.71* 1 0.64 0.82 0.26 0.11 0.26 0.21 0.36 0.46 -0.05 -0.27 0.29 0.07 0.16
selfacceptance -0.46* -0.53* 0.25 0.48* 0.7* 0.67* 0.64* 1 0.84 0.09 0.46 0.22 0.39 0.5 0.36 -0.22 -0.57 0 -0.03 0.22
pwb -0.31 -0.44 0.55* 0.61* 0.87* 0.83* 0.82* 0.84* 1 0.13 0.46 0.31 0.45 0.6 0.4 -0.24 -0.39 0.18 -0.23 0.23
ei 0.25 -0.1 0.06 -0.28 0.05 0.43 0.26 0.09 0.13 1 0.22 0.4 0.25 0.4 0.41 0.46 0.05 0.39 0.45 0.45
conformity 0.13 -0.22 0.44 0.4 0.3 0.41 0.11 0.46* 0.46* 0.22 1 0.6 0.69 0.8 0.6 0.04 -0.28 0.32 -0.07 0.68
tradition 0.23 -0.11 0.05 0.31 0.2 0.38 0.26 0.22 0.31 0.4 0.6* 1 0.41 0.64 0.58 0.09 -0.36 0.22 0.3 0.54
benevolence 0.08 -0.44 0.57* 0.1 0.17 0.6* 0.21 0.39 0.45* 0.25 0.69* 0.41 1 0.72 0.57 0.16 -0.02 0.36 -0.14 0.45
universalism 0.1 -0.34 0.42 0.47* 0.36 0.62* 0.36 0.5* 0.6* 0.4 0.8* 0.64* 0.72* 1 0.7 0.16 -0.23 0.44 0.07 0.78
selfdirection 0.12 0.03 0.17 0.23 0.07 0.52* 0.46* 0.36 0.4 0.41 0.6* 0.58* 0.57* 0.7* 1 0.39 0.02 0.58 0.34 0.8
stimulation 0.38 0.18 -0.18 -0.45* -0.29 0.12 -0.05 -0.22 -0.24 0.46* 0.04 0.09 0.16 0.16 0.39 1 0.47 0.34 0.58 0.4
hedonism 0.18 0.33 0.22 -0.51* -0.33 -0.27 -0.27 -0.57* -0.39 0.05 -0.28 -0.36 -0.02 -0.23 0.02 0.47* 1 0.28 0.01 -0.07
achievement 0.43 0.28 0.3 -0.03 0 0.24 0.29 0 0.18 0.39 0.32 0.22 0.36 0.44 0.58* 0.34 0.28 1 0.33 0.5
power 0.24 0.22 -0.52* -0.34 -0.22 -0.04 0.07 -0.03 -0.23 0.45* -0.07 0.3 -0.14 0.07 0.34 0.58* 0.01 0.33 1 0.4
security 0.25 0.03 0.12 0.2 0.01 0.33 0.16 0.22 0.23 0.45* 0.68* 0.54* 0.45* 0.78* 0.8* 0.4 -0.07 0.5* 0.4 1

regression model for env identity without factors interaction

env identity vs well-being

model <- lm(ei ~ positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance + sex + nationality, data = thes)
final_model <- step(model)  
## Start:  AIC=556.66
## ei ~ positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance + sex + nationality
## 
##                  Df Sum of Sq   RSS    AIC
## - nationality     1      0.08 15263 554.66
## - positiverelat   1      2.36 15265 554.68
## - purposeinlife   1    128.77 15392 555.58
## - autonomy        1    169.50 15433 555.86
## <none>                        15263 556.66
## - persgrowth      1    413.18 15676 557.57
## - envmattery      1    461.63 15725 557.91
## - selfacceptance  1    464.78 15728 557.93
## - sex             1    536.20 15799 558.42
## 
## Step:  AIC=554.66
## ei ~ positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance + sex
## 
##                  Df Sum of Sq   RSS    AIC
## - positiverelat   1      2.32 15265 552.68
## - purposeinlife   1    131.58 15395 553.60
## - autonomy        1    173.16 15436 553.89
## <none>                        15263 554.66
## - selfacceptance  1    464.70 15728 555.93
## - envmattery      1    484.89 15748 556.07
## - sex             1    540.56 15804 556.45
## - persgrowth      1    664.64 15928 557.31
## 
## Step:  AIC=552.68
## ei ~ autonomy + envmattery + persgrowth + purposeinlife + selfacceptance + 
##     sex
## 
##                  Df Sum of Sq   RSS    AIC
## - purposeinlife   1    130.04 15396 551.60
## - autonomy        1    170.96 15436 551.89
## <none>                        15265 552.68
## - envmattery      1    491.35 15757 554.13
## - selfacceptance  1    508.76 15774 554.25
## - sex             1    543.00 15808 554.49
## - persgrowth      1    731.21 15997 555.78
## 
## Step:  AIC=551.6
## ei ~ autonomy + envmattery + persgrowth + selfacceptance + sex
## 
##                  Df Sum of Sq   RSS    AIC
## - autonomy        1    135.49 15531 550.56
## <none>                        15396 551.60
## - selfacceptance  1    415.78 15811 552.51
## - envmattery      1    542.79 15938 553.38
## - sex             1    547.81 15943 553.41
## - persgrowth      1    680.13 16076 554.31
## 
## Step:  AIC=550.56
## ei ~ envmattery + persgrowth + selfacceptance + sex
## 
##                  Df Sum of Sq   RSS    AIC
## <none>                        15531 550.56
## - selfacceptance  1    383.11 15914 551.21
## - persgrowth      1    551.09 16082 552.36
## - envmattery      1    758.56 16290 553.76
## - sex             1    826.93 16358 554.21
summary(final_model)
## 
## Call:
## lm(formula = ei ~ envmattery + persgrowth + selfacceptance + 
##     sex, data = thes)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -30.5211  -7.4818   0.2539   7.8906  25.0681 
## 
## Coefficients:
##                Estimate Std. Error t value     Pr(>|t|)    
## (Intercept)     54.5098     9.1213   5.976 0.0000000324 ***
## envmattery      -0.4415     0.1959  -2.254       0.0263 *  
## persgrowth       0.3486     0.1814   1.921       0.0575 .  
## selfacceptance   0.1777     0.1109   1.602       0.1123    
## sexm            -5.6392     2.3965  -2.353       0.0205 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.22 on 104 degrees of freedom
## Multiple R-squared:  0.09838,    Adjusted R-squared:  0.0637 
## F-statistic: 2.837 on 4 and 104 DF,  p-value: 0.02803

env identity vs values

model <- lm(ei ~ conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security + sex + nationality, data = thes)
final_model <- step(model)  
## Start:  AIC=540.7
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security + 
##     sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - stimulation    1      0.04 12251 538.70
## - conformity     1      0.10 12251 538.70
## - achievement    1      0.19 12252 538.70
## - hedonism       1     57.82 12309 539.22
## - benevolence    1    124.51 12376 539.80
## - security       1    148.45 12400 540.02
## - selfdirection  1    181.56 12433 540.31
## - tradition      1    205.32 12457 540.51
## <none>                       12251 540.70
## - sex            1    233.39 12485 540.76
## - power          1    290.45 12542 541.26
## - nationality    1    655.19 12906 544.38
## - universalism   1    887.57 13139 546.33
## 
## Step:  AIC=538.7
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     hedonism + power + achievement + security + sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - conformity     1      0.08 12251 536.70
## - achievement    1      0.21 12252 536.70
## - hedonism       1     58.39 12310 537.22
## - benevolence    1    125.07 12376 537.81
## - security       1    149.91 12401 538.03
## - selfdirection  1    188.74 12440 538.37
## - tradition      1    205.38 12457 538.51
## <none>                       12251 538.70
## - sex            1    235.98 12487 538.78
## - power          1    316.75 12568 539.48
## - nationality    1    700.64 12952 542.76
## - universalism   1    904.09 13155 544.46
## 
## Step:  AIC=536.7
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     hedonism + power + achievement + security + sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - achievement    1      0.30 12252 534.71
## - hedonism       1     60.42 12312 535.24
## - benevolence    1    130.18 12382 535.86
## - security       1    164.34 12416 536.16
## - selfdirection  1    188.77 12440 536.37
## <none>                       12251 536.70
## - sex            1    235.99 12487 536.78
## - tradition      1    249.17 12501 536.90
## - power          1    318.53 12570 537.50
## - nationality    1    701.60 12953 540.77
## - universalism   1    905.23 13157 542.47
## 
## Step:  AIC=534.71
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     hedonism + power + security + sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - hedonism       1     60.72 12312 533.24
## - benevolence    1    144.53 12396 533.98
## - security       1    164.06 12416 534.16
## - selfdirection  1    202.16 12454 534.49
## <none>                       12252 534.71
## - sex            1    236.75 12488 534.79
## - tradition      1    248.88 12501 534.90
## - power          1    406.96 12659 536.27
## - nationality    1    708.46 12960 538.83
## - universalism   1    925.85 13178 540.65
## 
## Step:  AIC=533.24
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     power + security + sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - benevolence    1    170.85 12483 532.75
## - selfdirection  1    173.35 12486 532.77
## - tradition      1    202.90 12515 533.03
## <none>                       12312 533.24
## - sex            1    232.58 12545 533.28
## - security       1    237.05 12550 533.32
## - power          1    355.74 12668 534.35
## - nationality    1    669.20 12982 537.01
## - universalism   1    921.40 13234 539.11
## 
## Step:  AIC=532.75
## ei ~ tradition + universalism + selfdirection + power + security + 
##     sex + nationality
## 
##                 Df Sum of Sq   RSS    AIC
## - selfdirection  1     94.90 12578 531.57
## <none>                       12483 532.75
## - sex            1    276.77 12760 533.14
## - security       1    282.06 12765 533.18
## - tradition      1    311.17 12794 533.43
## - power          1    398.15 12881 534.17
## - nationality    1    606.93 13090 535.92
## - universalism   1   1219.21 13702 540.90
## 
## Step:  AIC=531.57
## ei ~ tradition + universalism + power + security + sex + nationality
## 
##                Df Sum of Sq   RSS    AIC
## <none>                      12578 531.57
## - security      1    233.87 12812 531.58
## - tradition     1    283.39 12862 532.00
## - sex           1    357.75 12936 532.63
## - power         1    461.90 13040 533.50
## - nationality   1    557.19 13135 534.30
## - universalism  1   1199.19 13777 539.50
summary(final_model)
## 
## Call:
## lm(formula = ei ~ tradition + universalism + power + security + 
##     sex + nationality, data = thes)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -30.9955  -6.3726   0.3801   7.2749  19.9856 
## 
## Coefficients:
##                Estimate Std. Error t value    Pr(>|t|)    
## (Intercept)     34.1416     6.0831   5.612 0.000000172 ***
## tradition        1.4427     0.9517   1.516     0.13263    
## universalism     3.6941     1.1846   3.118     0.00236 ** 
## power           -1.5456     0.7986  -1.935     0.05571 .  
## security         1.4435     1.0482   1.377     0.17148    
## sexm            -3.7811     2.2199  -1.703     0.09156 .  
## nationalityrus   5.1673     2.4309   2.126     0.03595 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 11.1 on 102 degrees of freedom
## Multiple R-squared:  0.2698, Adjusted R-squared:  0.2268 
## F-statistic: 6.281 on 6 and 102 DF,  p-value: 0.00001208

env identity vs well-being + values (without nationality)

model <- lm(ei ~ conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security + sex + positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance, data = thes)
final_model <- step(model)  
## Start:  AIC=545.71
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security + 
##     sex + positiverelat + autonomy + envmattery + persgrowth + 
##     purposeinlife + selfacceptance
## 
##                  Df Sum of Sq   RSS    AIC
## - purposeinlife   1      5.95 11709 543.76
## - achievement     1     10.40 11713 543.81
## - conformity      1     12.72 11716 543.83
## - selfacceptance  1     21.46 11724 543.91
## - stimulation     1     48.98 11752 544.16
## - sex             1     55.81 11758 544.23
## - benevolence     1     73.68 11776 544.39
## - security        1     82.30 11785 544.47
## - envmattery      1     99.09 11802 544.63
## - hedonism        1    112.49 11815 544.75
## - tradition       1    132.97 11836 544.94
## - positiverelat   1    137.72 11840 544.98
## - autonomy        1    186.13 11889 545.43
## <none>                        11703 545.71
## - selfdirection   1    235.77 11938 545.88
## - power           1    297.98 12001 546.45
## - universalism    1    876.45 12579 551.58
## - persgrowth      1   1005.95 12709 552.70
## 
## Step:  AIC=543.76
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security + 
##     sex + positiverelat + autonomy + envmattery + persgrowth + 
##     selfacceptance
## 
##                  Df Sum of Sq   RSS    AIC
## - achievement     1      6.06 11715 541.82
## - selfacceptance  1     15.55 11724 541.91
## - conformity      1     20.26 11729 541.95
## - stimulation     1     52.48 11761 542.25
## - sex             1     52.58 11761 542.25
## - security        1     77.10 11786 542.48
## - benevolence     1     87.01 11796 542.57
## - envmattery      1    113.83 11822 542.82
## - tradition       1    128.33 11837 542.95
## - hedonism        1    143.36 11852 543.09
## - positiverelat   1    167.61 11876 543.31
## - autonomy        1    180.61 11889 543.43
## <none>                        11709 543.76
## - selfdirection   1    281.93 11991 544.36
## - power           1    292.21 12001 544.45
## - universalism    1    893.13 12602 549.78
## - persgrowth      1   1059.08 12768 551.20
## 
## Step:  AIC=541.82
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + security + sex + positiverelat + 
##     autonomy + envmattery + persgrowth + selfacceptance
## 
##                  Df Sum of Sq   RSS    AIC
## - selfacceptance  1     14.05 11729 539.95
## - conformity      1     26.45 11741 540.07
## - stimulation     1     48.66 11763 540.27
## - sex             1     50.97 11766 540.29
## - security        1     76.09 11791 540.53
## - benevolence     1    105.45 11820 540.80
## - envmattery      1    108.66 11823 540.83
## - tradition       1    125.70 11840 540.98
## - hedonism        1    159.65 11874 541.30
## - positiverelat   1    162.52 11877 541.32
## - autonomy        1    181.27 11896 541.49
## <none>                        11715 541.82
## - selfdirection   1    277.61 11992 542.37
## - power           1    303.78 12018 542.61
## - universalism    1    890.78 12606 547.81
## - persgrowth      1   1059.10 12774 549.25
## 
## Step:  AIC=539.95
## ei ~ conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + security + sex + positiverelat + 
##     autonomy + envmattery + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - conformity     1     20.62 11749 538.14
## - stimulation    1     43.77 11773 538.36
## - sex            1     46.46 11775 538.38
## - security       1     92.80 11822 538.81
## - envmattery     1     94.64 11823 538.83
## - benevolence    1    115.08 11844 539.02
## - tradition      1    130.24 11859 539.16
## - hedonism       1    147.94 11877 539.32
## - positiverelat  1    150.35 11879 539.34
## - autonomy       1    175.83 11905 539.57
## <none>                       11729 539.95
## - selfdirection  1    263.56 11992 540.37
## - power          1    298.76 12028 540.69
## - universalism   1    918.66 12648 546.17
## - persgrowth     1   1087.48 12816 547.62
## 
## Step:  AIC=538.14
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + security + sex + positiverelat + 
##     autonomy + envmattery + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - sex            1     50.50 11800 536.61
## - stimulation    1     55.61 11805 536.66
## - envmattery     1     87.40 11837 536.95
## - positiverelat  1    140.20 11890 537.44
## - hedonism       1    140.58 11890 537.44
## - security       1    143.21 11893 537.46
## - benevolence    1    162.65 11912 537.64
## - autonomy       1    163.41 11913 537.65
## <none>                       11749 538.14
## - tradition      1    221.26 11971 538.18
## - selfdirection  1    254.50 12004 538.48
## - power          1    278.52 12028 538.70
## - universalism   1    925.03 12674 544.40
## - persgrowth     1   1068.19 12818 545.63
## 
## Step:  AIC=536.61
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + security + positiverelat + 
##     autonomy + envmattery + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - stimulation    1     68.17 11868 535.24
## - envmattery     1     74.82 11875 535.30
## - security       1    128.60 11928 535.79
## - hedonism       1    140.31 11940 535.90
## - positiverelat  1    147.03 11947 535.96
## - benevolence    1    180.04 11980 536.26
## <none>                       11800 536.61
## - tradition      1    225.12 12025 536.67
## - power          1    246.10 12046 536.86
## - autonomy       1    266.51 12066 537.05
## - selfdirection  1    308.44 12108 537.42
## - universalism   1   1125.74 12926 544.54
## - persgrowth     1   1195.03 12995 545.13
## 
## Step:  AIC=535.24
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     hedonism + power + security + positiverelat + autonomy + 
##     envmattery + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - envmattery     1     40.91 11909 533.61
## - hedonism       1    115.28 11983 534.29
## - security       1    125.94 11994 534.39
## - benevolence    1    164.00 12032 534.73
## - positiverelat  1    187.28 12055 534.95
## - tradition      1    217.23 12085 535.22
## <none>                       11868 535.24
## - autonomy       1    250.88 12119 535.52
## - selfdirection  1    367.94 12236 536.57
## - power          1    378.77 12247 536.66
## - universalism   1   1104.71 12973 542.94
## - persgrowth     1   1145.24 13013 543.28
## 
## Step:  AIC=533.61
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     hedonism + power + security + positiverelat + autonomy + 
##     persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - hedonism       1    114.44 12023 532.66
## - security       1    134.69 12044 532.84
## - benevolence    1    150.01 12059 532.98
## - tradition      1    195.78 12105 533.39
## <none>                       11909 533.61
## - positiverelat  1    221.33 12130 533.62
## - autonomy       1    322.28 12231 534.52
## - selfdirection  1    334.52 12244 534.63
## - power          1    393.60 12303 535.16
## - universalism   1   1227.44 13136 542.31
## - persgrowth     1   1298.85 13208 542.90
## 
## Step:  AIC=532.66
## ei ~ tradition + benevolence + universalism + selfdirection + 
##     power + security + positiverelat + autonomy + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## - tradition      1    130.08 12154 531.83
## - benevolence    1    184.39 12208 532.31
## - positiverelat  1    196.49 12220 532.42
## - security       1    212.30 12236 532.56
## <none>                       12023 532.66
## - selfdirection  1    273.72 12297 533.11
## - power          1    310.72 12334 533.44
## - autonomy       1    363.50 12387 533.90
## - persgrowth     1   1197.51 13221 541.00
## - universalism   1   1212.68 13236 541.13
## 
## Step:  AIC=531.83
## ei ~ benevolence + universalism + selfdirection + power + security + 
##     positiverelat + autonomy + persgrowth
## 
##                 Df Sum of Sq   RSS    AIC
## <none>                       12154 531.83
## - power          1    229.77 12383 531.87
## - security       1    245.70 12399 532.01
## - positiverelat  1    258.37 12412 532.12
## - selfdirection  1    270.76 12424 532.23
## - benevolence    1    313.04 12467 532.60
## - autonomy       1    325.21 12479 532.71
## - persgrowth     1   1079.61 13233 539.11
## - universalism   1   1487.59 13641 542.42
summary(final_model)
## 
## Call:
## lm(formula = ei ~ benevolence + universalism + selfdirection + 
##     power + security + positiverelat + autonomy + persgrowth, 
##     data = thes)
## 
## 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 ** 
## benevolence     2.2434     1.3979   1.605  0.11167    
## universalism    4.6865     1.3396   3.499  0.00070 ***
## selfdirection  -2.2173     1.4855  -1.493  0.13869    
## power          -1.0600     0.7710  -1.375  0.17221    
## security        1.5156     1.0659   1.422  0.15818    
## positiverelat  -0.2382     0.1634  -1.458  0.14797    
## autonomy       -0.2588     0.1582  -1.636  0.10502    
## persgrowth      0.4556     0.1529   2.980  0.00361 ** 
## ---
## 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

regression model for env identity with factors interaction

env identity vs well-being:sex

model <- lm(ei ~ (positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance):(sex + nationality), data = thes)
final_model <- step(model)  
## Start:  AIC=567.5
## ei ~ (positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance):(sex + nationality)
## 
##                              Df Sum of Sq   RSS    AIC
## - positiverelat:sex           2    205.12 14238 565.08
## - autonomy:sex                2    240.76 14274 565.36
## - positiverelat:nationality   1      0.01 14033 565.50
## - envmattery:nationality      1     15.03 14048 565.62
## - persgrowth:nationality      1     15.22 14048 565.62
## - persgrowth:sex              2    317.61 14351 565.94
## - selfacceptance:nationality  1    155.19 14188 566.70
## - autonomy:nationality        1    165.00 14198 566.78
## - envmattery:sex              2    438.25 14471 566.85
## - purposeinlife:sex           2    491.64 14525 567.26
## - purposeinlife:nationality   1    238.35 14271 567.34
## <none>                                    14033 567.50
## - selfacceptance:sex          2    623.35 14656 568.24
## 
## Step:  AIC=566.98
## ei ~ positiverelat:nationality + autonomy:sex + autonomy:nationality + 
##     envmattery:sex + envmattery:nationality + persgrowth:sex + 
##     persgrowth:nationality + purposeinlife:sex + purposeinlife:nationality + 
##     selfacceptance:sex + selfacceptance:nationality
## 
##                              Df Sum of Sq   RSS    AIC
## - positiverelat:nationality   2     27.17 14252 563.19
## - sex:persgrowth              2    228.77 14453 564.72
## - nationality:envmattery      1      1.29 14226 564.99
## - nationality:persgrowth      1     21.01 14246 565.14
## - nationality:autonomy        1    101.17 14326 565.75
## - nationality:selfacceptance  1    111.10 14336 565.83
## - autonomy:sex                2    380.54 14605 565.86
## - sex:envmattery              2    454.34 14679 566.41
## - sex:purposeinlife           2    508.96 14734 566.81
## <none>                                    14225 566.98
## - nationality:purposeinlife   1    346.08 14571 567.60
## - sex:selfacceptance          2    697.81 14922 568.20
## 
## Step:  AIC=563.19
## ei ~ autonomy:sex + nationality:autonomy + sex:envmattery + nationality:envmattery + 
##     sex:persgrowth + nationality:persgrowth + sex:purposeinlife + 
##     nationality:purposeinlife + sex:selfacceptance + nationality:selfacceptance
## 
##                              Df Sum of Sq   RSS    AIC
## - sex:persgrowth              2    255.98 14508 561.13
## - nationality:envmattery      1      1.21 14253 561.20
## - nationality:persgrowth      1      6.91 14259 561.24
## - autonomy:sex                2    380.65 14632 562.06
## - nationality:selfacceptance  1    131.84 14384 562.19
## - autonomy:nationality        1    143.66 14395 562.28
## - sex:envmattery              2    484.93 14737 562.84
## - sex:purposeinlife           2    507.57 14759 563.00
## <none>                                    14252 563.19
## - nationality:purposeinlife   1    331.96 14584 563.70
## - sex:selfacceptance          2    801.59 15053 565.15
## 
## Step:  AIC=562.4
## ei ~ autonomy:sex + autonomy:nationality + sex:envmattery + nationality:envmattery + 
##     nationality:persgrowth + sex:purposeinlife + nationality:purposeinlife + 
##     sex:selfacceptance + nationality:selfacceptance
## 
##                              Df Sum of Sq   RSS    AIC
## - autonomy:sex                2    237.46 14648 560.18
## - nationality:envmattery      1      3.36 14414 560.42
## - nationality:persgrowth      2    302.57 14713 560.66
## - nationality:selfacceptance  1    127.85 14539 561.36
## - autonomy:nationality        1    168.28 14579 561.66
## - sex:envmattery              2    474.53 14885 561.93
## <none>                                    14411 562.40
## - sex:purposeinlife           2    542.72 14954 562.43
## - nationality:purposeinlife   1    368.88 14780 563.15
## - sex:selfacceptance          2    817.69 15228 564.41
## 
## Step:  AIC=562.15
## ei ~ autonomy:nationality + sex:envmattery + nationality:envmattery + 
##     nationality:persgrowth + sex:purposeinlife + nationality:purposeinlife + 
##     sex:selfacceptance + nationality:selfacceptance
## 
##                              Df Sum of Sq   RSS    AIC
## - nationality:envmattery      1      3.30 14647 560.17
## - sex:envmattery              2    299.46 14943 560.35
## - nationality:selfacceptance  1     96.83 14741 560.87
## - autonomy:nationality        2    413.31 15057 561.18
## - nationality:persgrowth      2    416.64 15061 561.21
## - sex:purposeinlife           2    452.00 15096 561.46
## <none>                                    14644 562.15
## - nationality:purposeinlife   1    334.08 14978 562.61
## - sex:selfacceptance          2    750.47 15394 563.60
## 
## Step:  AIC=560.17
## ei ~ autonomy:nationality + sex:envmattery + nationality:persgrowth + 
##     sex:purposeinlife + nationality:purposeinlife + sex:selfacceptance + 
##     nationality:selfacceptance
## 
##                              Df Sum of Sq   RSS    AIC
## - nationality:selfacceptance  1    109.06 14756 558.98
## - nationality:persgrowth      2    414.23 15062 559.21
## - autonomy:nationality        2    438.67 15086 559.39
## - sex:purposeinlife           2    454.99 15102 559.51
## <none>                                    14647 560.17
## - sex:envmattery              2    585.24 15232 560.44
## - nationality:purposeinlife   1    335.36 14983 560.64
## - sex:selfacceptance          2    776.24 15424 561.80
## 
## Step:  AIC=558.98
## ei ~ autonomy:nationality + sex:envmattery + nationality:persgrowth + 
##     sex:purposeinlife + nationality:purposeinlife + sex:selfacceptance
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:persgrowth     2    322.60 15079 557.34
## - sex:purposeinlife          2    346.93 15103 557.51
## - nationality:purposeinlife  1    226.66 14983 558.64
## <none>                                   14756 558.98
## - autonomy:nationality       2    575.12 15331 559.15
## - sex:envmattery             2    688.38 15445 559.95
## - sex:selfacceptance         2    757.46 15514 560.44
## 
## Step:  AIC=557.34
## ei ~ autonomy:nationality + sex:envmattery + sex:purposeinlife + 
##     nationality:purposeinlife + sex:selfacceptance
## 
##                             Df Sum of Sq   RSS    AIC
## - sex:purposeinlife          2    468.54 15548 556.67
## <none>                                   15079 557.34
## - sex:envmattery             2    603.18 15682 557.61
## - autonomy:nationality       2    611.33 15690 557.67
## - sex:selfacceptance         2    865.95 15945 559.42
## - nationality:purposeinlife  1    710.50 15789 560.36
## 
## Step:  AIC=555.73
## ei ~ autonomy:nationality + sex:envmattery + nationality:purposeinlife + 
##     sex:selfacceptance
## 
##                             Df Sum of Sq   RSS    AIC
## <none>                                   15133 555.73
## - sex:envmattery             2    608.29 15741 556.02
## - autonomy:nationality       2    652.68 15786 556.33
## - nationality:purposeinlife  2    792.80 15926 557.29
## - sex:selfacceptance         2    823.43 15956 557.50
summary(final_model)
## 
## Call:
## lm(formula = ei ~ autonomy:nationality + sex:envmattery + nationality:purposeinlife + 
##     sex:selfacceptance, data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -32.049  -7.486   0.569   8.177  27.418 
## 
## Coefficients:
##                              Estimate Std. Error t value     Pr(>|t|)    
## (Intercept)                  75.38151   12.96134   5.816 0.0000000726 ***
## autonomy:nationalityfr        0.05325    0.24604   0.216       0.8291    
## autonomy:nationalityrus      -0.44153    0.23755  -1.859       0.0660 .  
## sexf:envmattery              -0.49861    0.25127  -1.984       0.0500 *  
## sexm:envmattery              -0.32935    0.23650  -1.393       0.1668    
## nationalityfr:purposeinlife  -0.33509    0.20244  -1.655       0.1010    
## nationalityrus:purposeinlife  0.22400    0.23876   0.938       0.3504    
## sexf:selfacceptance           0.45371    0.19576   2.318       0.0225 *  
## sexm:selfacceptance           0.24351    0.17843   1.365       0.1754    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.3 on 100 degrees of freedom
## Multiple R-squared:  0.1215, Adjusted R-squared:  0.05121 
## F-statistic: 1.729 on 8 and 100 DF,  p-value: 0.1008
model <- lm(ei ~ (positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance):(nationality*sex), data = thes)
final_model <- step(model)  
## Start:  AIC=574.66
## ei ~ (positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance):(nationality * sex)
## 
##                                  Df Sum of Sq   RSS    AIC
## - persgrowth:nationality:sex      1     1.644 13425 572.68
## - envmattery:nationality:sex      1    20.815 13444 572.83
## - positiverelat:nationality:sex   1    23.341 13447 572.85
## - selfacceptance:nationality:sex  1    37.491 13461 572.97
## - autonomy:nationality:sex        1    77.770 13502 573.29
## - purposeinlife:nationality:sex   1   120.288 13544 573.64
## <none>                                        13424 574.66
## 
## Step:  AIC=572.68
## ei ~ positiverelat:nationality + positiverelat:sex + autonomy:nationality + 
##     autonomy:sex + envmattery:nationality + envmattery:sex + 
##     persgrowth:nationality + persgrowth:sex + purposeinlife:nationality + 
##     purposeinlife:sex + selfacceptance:nationality + selfacceptance:sex + 
##     positiverelat:nationality:sex + autonomy:nationality:sex + 
##     envmattery:nationality:sex + purposeinlife:nationality:sex + 
##     selfacceptance:nationality:sex
## 
##                                  Df Sum of Sq   RSS    AIC
## - nationality:persgrowth          2   102.170 13528 569.50
## - nationality:sex:envmattery      1    19.298 13445 570.83
## - positiverelat:nationality:sex   1    21.770 13447 570.85
## - nationality:sex:selfacceptance  1    36.378 13462 570.97
## - nationality:sex:autonomy        1    78.229 13504 571.31
## - nationality:sex:purposeinlife   1   131.611 13557 571.74
## - sex:persgrowth                  1   147.892 13573 571.87
## <none>                                        13425 572.68
## 
## Step:  AIC=571.03
## ei ~ positiverelat:nationality + positiverelat:sex + nationality:autonomy + 
##     sex:autonomy + nationality:envmattery + sex:envmattery + 
##     sex:persgrowth + nationality:purposeinlife + sex:purposeinlife + 
##     nationality:selfacceptance + sex:selfacceptance + positiverelat:nationality:sex + 
##     nationality:sex:autonomy + nationality:sex:envmattery + nationality:sex:purposeinlife + 
##     nationality:sex:selfacceptance
## 
##                                  Df Sum of Sq   RSS    AIC
## - positiverelat:nationality:sex   1     17.04 13486 569.17
## - nationality:sex:envmattery      1     17.84 13487 569.18
## - nationality:sex:selfacceptance  1     27.87 13497 569.26
## - nationality:sex:autonomy        1     78.27 13547 569.66
## - nationality:sex:purposeinlife   1    127.62 13597 570.06
## - sex:persgrowth                  2    459.35 13928 570.69
## <none>                                        13469 571.03
## 
## Step:  AIC=569.17
## ei ~ positiverelat:nationality + positiverelat:sex + nationality:autonomy + 
##     sex:autonomy + nationality:envmattery + sex:envmattery + 
##     sex:persgrowth + nationality:purposeinlife + sex:purposeinlife + 
##     nationality:selfacceptance + sex:selfacceptance + nationality:sex:autonomy + 
##     nationality:sex:envmattery + nationality:sex:purposeinlife + 
##     nationality:sex:selfacceptance
## 
##                                  Df Sum of Sq   RSS    AIC
## - positiverelat:nationality       2     62.84 13549 565.68
## - nationality:sex:selfacceptance  1     22.97 13509 567.36
## - nationality:sex:envmattery      1     25.99 13512 567.38
## - nationality:sex:purposeinlife   1    111.72 13598 568.07
## - positiverelat:sex               1    120.48 13607 568.14
## - nationality:sex:autonomy        1    174.90 13661 568.57
## - sex:persgrowth                  2    462.52 13949 568.85
## <none>                                        13486 569.17
## 
## Step:  AIC=567.17
## ei ~ positiverelat:sex + nationality:autonomy + sex:autonomy + 
##     nationality:envmattery + sex:envmattery + sex:persgrowth + 
##     nationality:purposeinlife + sex:purposeinlife + nationality:selfacceptance + 
##     sex:selfacceptance + nationality:sex:autonomy + nationality:sex:envmattery + 
##     nationality:sex:purposeinlife + nationality:sex:selfacceptance
## 
##                                  Df Sum of Sq   RSS    AIC
## - positiverelat:sex               2    130.02 13616 564.22
## - sex:nationality:selfacceptance  1     24.59 13511 565.37
## - sex:nationality:envmattery      1     27.57 13514 565.39
## - sex:nationality:purposeinlife   1    111.52 13598 566.07
## - sex:nationality:autonomy        1    174.96 13661 566.58
## - sex:persgrowth                  2    465.21 13952 566.87
## <none>                                        13486 567.17
## 
## Step:  AIC=564.22
## ei ~ nationality:autonomy + sex:autonomy + nationality:envmattery + 
##     sex:envmattery + sex:persgrowth + nationality:purposeinlife + 
##     sex:purposeinlife + nationality:selfacceptance + sex:selfacceptance + 
##     sex:nationality:autonomy + sex:nationality:envmattery + sex:nationality:purposeinlife + 
##     sex:nationality:selfacceptance
## 
##                                  Df Sum of Sq   RSS    AIC
## - nationality:sex:envmattery      1     10.95 13627 562.30
## - nationality:sex:selfacceptance  1     16.38 13633 562.35
## - nationality:sex:purposeinlife   1     83.72 13700 562.89
## - sex:persgrowth                  2    436.61 14053 563.66
## - nationality:autonomy:sex        1    187.24 13804 563.71
## <none>                                        13616 564.22
## 
## Step:  AIC=562.3
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     sex:envmattery + sex:persgrowth + nationality:purposeinlife + 
##     sex:purposeinlife + nationality:selfacceptance + sex:selfacceptance + 
##     nationality:autonomy:sex + nationality:sex:purposeinlife + 
##     nationality:sex:selfacceptance
## 
##                                  Df Sum of Sq   RSS    AIC
## - nationality:sex:selfacceptance  1     12.04 13639 560.40
## - sex:envmattery                  1     15.52 13643 560.43
## - nationality:sex:purposeinlife   1     72.81 13700 560.89
## - sex:persgrowth                  2    426.32 14054 561.66
## <none>                                        13627 562.30
## - nationality:envmattery          2    537.95 14165 562.52
## - nationality:autonomy:sex        1    320.97 13948 562.84
## 
## Step:  AIC=560.4
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     sex:envmattery + sex:persgrowth + nationality:purposeinlife + 
##     sex:purposeinlife + nationality:selfacceptance + sex:selfacceptance + 
##     nationality:autonomy:sex + nationality:sex:purposeinlife
## 
##                                 Df Sum of Sq   RSS    AIC
## - sex:envmattery                 1     13.00 13652 558.50
## - sex:selfacceptance             1    149.79 13789 559.59
## - sex:persgrowth                 2    417.84 14057 559.69
## - nationality:sex:purposeinlife  1    187.62 13827 559.89
## <none>                                       13639 560.40
## - nationality:envmattery         2    548.62 14188 560.70
## - nationality:autonomy:sex       1    317.55 13957 560.91
## - nationality:selfacceptance     2    952.38 14592 563.76
## 
## Step:  AIC=558.5
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     sex:persgrowth + nationality:purposeinlife + sex:purposeinlife + 
##     nationality:selfacceptance + sex:selfacceptance + nationality:autonomy:sex + 
##     nationality:sex:purposeinlife
## 
##                                 Df Sum of Sq   RSS    AIC
## - sex:selfacceptance             1    140.49 13793 557.62
## - nationality:sex:purposeinlife  1    177.31 13830 557.91
## <none>                                       13652 558.50
## - sex:persgrowth                 2    522.14 14174 558.60
## - nationality:autonomy:sex       1    306.75 13959 558.93
## - nationality:envmattery         2    786.37 14439 560.61
## - nationality:selfacceptance     2    948.68 14601 561.83
## 
## Step:  AIC=557.62
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     sex:persgrowth + nationality:purposeinlife + sex:purposeinlife + 
##     nationality:selfacceptance + nationality:autonomy:sex + nationality:sex:purposeinlife
## 
##                                 Df Sum of Sq   RSS    AIC
## - nationality:sex:purposeinlife  1    114.33 13907 556.52
## - sex:persgrowth                 2    485.65 14278 557.39
## - nationality:autonomy:sex       1    228.86 14022 557.41
## <none>                                       13793 557.62
## - nationality:envmattery         2    702.61 14496 559.04
## - nationality:selfacceptance     2    952.48 14745 560.90
## 
## Step:  AIC=556.52
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     sex:persgrowth + nationality:purposeinlife + sex:purposeinlife + 
##     nationality:selfacceptance + nationality:autonomy:sex
## 
##                              Df Sum of Sq   RSS    AIC
## - sex:persgrowth              2    379.56 14287 555.46
## - sex:purposeinlife           1    159.90 14067 555.77
## <none>                                    13907 556.52
## - nationality:envmattery      2    659.93 14567 557.57
## - nationality:purposeinlife   2    696.32 14604 557.85
## - nationality:autonomy:sex    1    454.83 14362 558.03
## - nationality:selfacceptance  2    890.67 14798 559.29
## 
## Step:  AIC=555.46
## ei ~ nationality:autonomy + autonomy:sex + nationality:envmattery + 
##     nationality:purposeinlife + sex:purposeinlife + nationality:selfacceptance + 
##     nationality:autonomy:sex
## 
##                              Df Sum of Sq   RSS    AIC
## - nationality:envmattery      2    511.93 14799 555.29
## <none>                                    14287 555.46
## - sex:purposeinlife           1    354.37 14641 556.13
## - nationality:purposeinlife   2    655.27 14942 556.34
## - nationality:selfacceptance  2    750.64 15037 557.04
## - nationality:autonomy:sex    1    840.88 15128 559.69
## 
## Step:  AIC=555.29
## ei ~ nationality:autonomy + autonomy:sex + nationality:purposeinlife + 
##     sex:purposeinlife + nationality:selfacceptance + nationality:autonomy:sex
## 
##                              Df Sum of Sq   RSS    AIC
## <none>                                    14799 555.29
## - sex:purposeinlife           1    311.17 15110 555.56
## - nationality:selfacceptance  2    699.49 15498 556.33
## - nationality:purposeinlife   2    749.10 15548 556.68
## - nationality:autonomy:sex    1    821.74 15620 559.18
summary(final_model)
## 
## Call:
## lm(formula = ei ~ nationality:autonomy + autonomy:sex + nationality:purposeinlife + 
##     sex:purposeinlife + nationality:selfacceptance + nationality:autonomy:sex, 
##     data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -33.759  -7.385   1.614   7.913  25.847 
## 
## Coefficients:
##                               Estimate Std. Error t value    Pr(>|t|)    
## (Intercept)                   66.83479   11.74503   5.690 0.000000129 ***
## nationalityfr:autonomy         0.08789    0.26242   0.335      0.7384    
## nationalityrus:autonomy       -0.38815    0.30550  -1.271      0.2069    
## autonomy:sexm                 -0.53940    0.27475  -1.963      0.0524 .  
## nationalityfr:purposeinlife   -0.56586    0.26727  -2.117      0.0368 *  
## nationalityrus:purposeinlife   0.06158    0.27252   0.226      0.8217    
## sexm:purposeinlife             0.31392    0.21758   1.443      0.1522    
## nationalityfr:selfacceptance   0.43592    0.22107   1.972      0.0514 .  
## nationalityrus:selfacceptance  0.20240    0.22869   0.885      0.3783    
## nationalityrus:autonomy:sexm   0.25555    0.10899   2.345      0.0210 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.23 on 99 degrees of freedom
## Multiple R-squared:  0.1409, Adjusted R-squared:  0.06279 
## F-statistic: 1.804 on 9 and 99 DF,  p-value: 0.07679

env identity vs values:sex,:nationality,:sex and nationality

model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security) : sex, data = thes)
final_model <- step(model)  
## Start:  AIC=551.09
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security):sex
## 
##                     Df Sum of Sq   RSS    AIC
## - tradition:sex      2     46.69 11683 547.53
## - achievement:sex    2     47.85 11684 547.54
## - stimulation:sex    2     50.44 11687 547.56
## - benevolence:sex    2    187.40 11824 548.83
## - security:sex       2    240.40 11877 549.32
## - selfdirection:sex  2    290.54 11927 549.78
## - conformity:sex     2    324.06 11961 550.09
## - power:sex          2    416.78 12053 550.93
## <none>                           11637 551.09
## - hedonism:sex       2    437.67 12074 551.12
## - universalism:sex   2    708.41 12345 553.53
## 
## Step:  AIC=547.53
## ei ~ conformity:sex + benevolence:sex + universalism:sex + selfdirection:sex + 
##     stimulation:sex + hedonism:sex + power:sex + achievement:sex + 
##     security:sex
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:achievement    2     30.44 11714 543.81
## - sex:stimulation    2     53.85 11737 544.03
## - sex:benevolence    2    199.57 11883 545.37
## - sex:security       2    281.01 11964 546.12
## - sex:selfdirection  2    290.02 11973 546.20
## - conformity:sex     2    321.67 12005 546.49
## <none>                           11683 547.53
## - sex:power          2    438.69 12122 547.55
## - sex:hedonism       2    495.69 12179 548.06
## - sex:universalism   2    734.40 12418 550.17
## 
## Step:  AIC=543.81
## ei ~ conformity:sex + sex:benevolence + sex:universalism + sex:selfdirection + 
##     sex:stimulation + sex:hedonism + sex:power + sex:security
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:stimulation    2     53.22 11767 540.31
## - sex:benevolence    2    258.00 11972 542.19
## - sex:security       2    280.59 11994 542.39
## - conformity:sex     2    303.01 12017 542.60
## - sex:selfdirection  2    314.56 12028 542.70
## <none>                           11714 543.81
## - sex:hedonism       2    488.98 12203 544.27
## - sex:power          2    562.79 12276 544.93
## - sex:universalism   2    730.29 12444 546.40
## 
## Step:  AIC=540.31
## ei ~ conformity:sex + sex:benevolence + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:benevolence    2    226.98 11994 538.39
## - sex:security       2    237.61 12004 538.48
## - conformity:sex     2    263.46 12030 538.72
## - sex:selfdirection  2    347.94 12115 539.48
## <none>                           11767 540.31
## - sex:hedonism       2    561.23 12328 541.38
## - sex:power          2    660.93 12428 542.26
## - sex:universalism   2    732.27 12499 542.89
## 
## Step:  AIC=538.39
## ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:security       2    307.61 12302 537.15
## - sex:selfdirection  2    351.22 12345 537.53
## - conformity:sex     2    363.40 12357 537.64
## <none>                           11994 538.39
## - sex:hedonism       2    526.52 12520 539.07
## - sex:power          2    648.94 12643 540.13
## - sex:universalism   2    988.97 12983 543.02
## 
## Step:  AIC=537.15
## ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power
## 
##                     Df Sum of Sq   RSS    AIC
## - conformity:sex     2    372.79 12674 536.40
## - sex:selfdirection  2    432.70 12734 536.92
## <none>                           12302 537.15
## - sex:hedonism       2    517.78 12819 537.64
## - sex:power          2    674.25 12976 538.96
## - sex:universalism   2   1082.75 13384 542.34
## 
## Step:  AIC=536.4
## ei ~ sex:universalism + sex:selfdirection + sex:hedonism + sex:power
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:selfdirection  2    325.92 13000 535.17
## - sex:power          2    375.33 13050 535.58
## <none>                           12674 536.40
## - sex:hedonism       2    536.95 13211 536.92
## - sex:universalism   2   1669.51 14344 545.89
## 
## Step:  AIC=535.17
## ei ~ sex:universalism + sex:hedonism + sex:power
## 
##                    Df Sum of Sq   RSS    AIC
## - sex:hedonism      2    337.93 13338 533.97
## <none>                          13000 535.17
## - sex:power         2    532.74 13533 535.55
## - sex:universalism  2   2185.31 15186 548.11
## 
## Step:  AIC=533.97
## ei ~ sex:universalism + sex:power
## 
##                    Df Sum of Sq   RSS    AIC
## - sex:power         2     228.0 13566 531.81
## <none>                          13338 533.97
## - sex:universalism  2    3542.4 16881 555.64
## 
## Step:  AIC=531.81
## ei ~ sex:universalism
## 
##                    Df Sum of Sq   RSS    AIC
## <none>                          13566 531.81
## - sex:universalism  2    3659.5 17226 553.85
summary(final_model)
## 
## Call:
## lm(formula = ei ~ sex:universalism, data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -29.847  -7.287   0.851   7.153  21.177 
## 
## Coefficients:
##                   Estimate Std. Error t value        Pr(>|t|)    
## (Intercept)        36.8803     4.9418   7.463 0.0000000000247 ***
## sexf:universalism   4.8987     0.9558   5.125 0.0000013442810 ***
## sexm:universalism   4.2148     1.0547   3.996        0.000119 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 11.31 on 106 degrees of freedom
## Multiple R-squared:  0.2124, Adjusted R-squared:  0.1976 
## F-statistic:  14.3 on 2 and 106 DF,  p-value: 0.000003183
model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security) : nationality, data = thes)
final_model <- step(model)  
## Start:  AIC=544.15
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security):nationality
## 
##                             Df Sum of Sq   RSS    AIC
## - stimulation:nationality    2     10.75 10929 540.25
## - achievement:nationality    2     61.33 10980 540.76
## - conformity:nationality     2    131.20 11050 541.45
## - tradition:nationality      2    140.60 11059 541.54
## - hedonism:nationality       2    203.42 11122 542.16
## - selfdirection:nationality  2    395.06 11313 544.02
## - power:nationality          2    397.83 11316 544.05
## - security:nationality       2    407.36 11326 544.14
## <none>                                   10918 544.15
## - benevolence:nationality    2    501.21 11420 545.04
## - universalism:nationality   2    762.70 11681 547.51
## 
## Step:  AIC=540.25
## ei ~ conformity:nationality + tradition:nationality + benevolence:nationality + 
##     universalism:nationality + selfdirection:nationality + hedonism:nationality + 
##     power:nationality + achievement:nationality + security:nationality
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:achievement    2     70.77 11000 536.96
## - nationality:tradition      2    141.70 11071 537.66
## - conformity:nationality     2    169.32 11098 537.93
## - nationality:hedonism       2    204.84 11134 538.28
## - nationality:security       2    401.80 11331 540.19
## <none>                                   10929 540.25
## - nationality:selfdirection  2    427.41 11356 540.43
## - nationality:power          2    449.16 11378 540.64
## - nationality:benevolence    2    491.79 11421 541.05
## - nationality:universalism   2    767.11 11696 543.65
## 
## Step:  AIC=536.96
## ei ~ conformity:nationality + nationality:tradition + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + nationality:hedonism + 
##     nationality:power + nationality:security
## 
##                             Df Sum of Sq   RSS    AIC
## - conformity:nationality     2    131.44 11131 534.25
## - nationality:tradition      2    133.19 11133 534.27
## - nationality:hedonism       2    206.90 11207 534.99
## - nationality:security       2    390.93 11391 536.76
## <none>                                   11000 536.96
## - nationality:benevolence    2    449.64 11449 537.32
## - nationality:selfdirection  2    462.94 11463 537.45
## - nationality:power          2    600.45 11600 538.75
## - nationality:universalism   2    835.45 11835 540.94
## 
## Step:  AIC=534.25
## ei ~ nationality:tradition + nationality:benevolence + nationality:universalism + 
##     nationality:selfdirection + nationality:hedonism + nationality:power + 
##     nationality:security
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:tradition      2    214.98 11346 532.34
## - nationality:hedonism       2    221.70 11353 532.40
## <none>                                   11131 534.25
## - nationality:selfdirection  2    425.43 11557 534.34
## - nationality:security       2    437.96 11569 534.46
## - nationality:power          2    485.37 11617 534.90
## - nationality:benevolence    2    599.99 11731 535.97
## - nationality:universalism   2    953.12 12084 539.21
## 
## Step:  AIC=532.34
## ei ~ nationality:benevolence + nationality:universalism + nationality:selfdirection + 
##     nationality:hedonism + nationality:power + nationality:security
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:hedonism       2    165.69 11512 529.92
## - nationality:power          2    373.24 11719 531.86
## - nationality:selfdirection  2    422.98 11769 532.33
## <none>                                   11346 532.34
## - nationality:security       2    480.90 11827 532.86
## - nationality:benevolence    2    841.55 12188 536.14
## - nationality:universalism   2   1176.73 12523 539.09
## 
## Step:  AIC=529.92
## ei ~ nationality:benevolence + nationality:universalism + nationality:selfdirection + 
##     nationality:power + nationality:security
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:power          2    413.12 11925 529.76
## <none>                                   11512 529.92
## - nationality:selfdirection  2    482.27 11994 530.39
## - nationality:security       2    586.84 12099 531.34
## - nationality:benevolence    2    803.35 12315 533.27
## - nationality:universalism   2   1093.94 12606 535.81
## 
## Step:  AIC=529.76
## ei ~ nationality:benevolence + nationality:universalism + nationality:selfdirection + 
##     nationality:security
## 
##                             Df Sum of Sq   RSS    AIC
## <none>                                   11925 529.76
## - nationality:security       2    563.16 12488 530.79
## - nationality:selfdirection  2    625.50 12550 531.33
## - nationality:benevolence    2    859.49 12784 533.35
## - nationality:universalism   2   1335.75 13261 537.33
summary(final_model)
## 
## Call:
## lm(formula = ei ~ nationality:benevolence + nationality:universalism + 
##     nationality:selfdirection + nationality:security, data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -26.284  -6.471   1.428   7.015  21.762 
## 
## Coefficients:
##                               Estimate Std. Error t value  Pr(>|t|)    
## (Intercept)                  27.427308   6.676812   4.108 0.0000818 ***
## nationalityfr:benevolence     5.295566   1.972565   2.685   0.00850 ** 
## nationalityrus:benevolence    0.148816   1.910502   0.078   0.93807    
## nationalityfr:universalism    5.295505   1.722797   3.074   0.00272 ** 
## nationalityrus:universalism   2.778431   2.143030   1.296   0.19779    
## nationalityfr:selfdirection  -4.216546   1.845824  -2.284   0.02446 *  
## nationalityrus:selfdirection  0.004102   2.166130   0.002   0.99849    
## nationalityfr:security       -0.059059   1.105326  -0.053   0.95749    
## nationalityrus:security       4.512261   2.077348   2.172   0.03221 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.92 on 100 degrees of freedom
## Multiple R-squared:  0.3077, Adjusted R-squared:  0.2523 
## F-statistic: 5.556 on 8 and 100 DF,  p-value: 0.00000807
model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security) : (sex*nationality), data = thes)
final_model <- step(model)  
## Start:  AIC=556.44
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security):(sex * 
##     nationality)
## 
##                                 Df Sum of Sq    RSS    AIC
## - benevolence:sex:nationality    1      1.84 8469.8 554.47
## - conformity:sex:nationality     1      2.31 8470.2 554.47
## - universalism:sex:nationality   1      2.59 8470.5 554.48
## - tradition:sex:nationality      1     19.65 8487.6 554.70
## - selfdirection:sex:nationality  1     24.34 8492.3 554.76
## - stimulation:sex:nationality    1     64.13 8532.0 555.27
## - security:sex:nationality       1    103.15 8571.1 555.76
## - hedonism:sex:nationality       1    109.19 8577.1 555.84
## - achievement:sex:nationality    1    119.03 8586.9 555.96
## <none>                                       8467.9 556.44
## - power:sex:nationality          1    358.84 8826.8 558.97
## 
## Step:  AIC=554.47
## ei ~ conformity:sex + conformity:nationality + tradition:sex + 
##     tradition:nationality + benevolence:sex + benevolence:nationality + 
##     universalism:sex + universalism:nationality + selfdirection:sex + 
##     selfdirection:nationality + stimulation:sex + stimulation:nationality + 
##     hedonism:sex + hedonism:nationality + power:sex + power:nationality + 
##     achievement:sex + achievement:nationality + security:sex + 
##     security:nationality + conformity:sex:nationality + tradition:sex:nationality + 
##     universalism:sex:nationality + selfdirection:sex:nationality + 
##     stimulation:sex:nationality + hedonism:sex:nationality + 
##     power:sex:nationality + achievement:sex:nationality + security:sex:nationality
## 
##                                 Df Sum of Sq    RSS    AIC
## - conformity:sex:nationality     1      1.65 8471.4 552.49
## - sex:nationality:universalism   1      7.23 8477.0 552.56
## - sex:nationality:tradition      1     20.13 8489.9 552.73
## - sex:nationality:selfdirection  1     22.71 8492.5 552.76
## - sex:nationality:stimulation    1     66.23 8536.0 553.32
## - sex:nationality:security       1    107.72 8577.5 553.84
## - sex:nationality:hedonism       1    116.81 8586.6 553.96
## - sex:nationality:achievement    1    128.11 8597.9 554.10
## <none>                                       8469.8 554.47
## - sex:nationality:power          1    374.01 8843.8 557.18
## - sex:benevolence                2    549.88 9019.6 557.32
## - nationality:benevolence        1    395.41 8865.2 557.44
## 
## Step:  AIC=552.49
## ei ~ conformity:sex + conformity:nationality + sex:tradition + 
##     nationality:tradition + sex:benevolence + nationality:benevolence + 
##     sex:universalism + nationality:universalism + sex:selfdirection + 
##     nationality:selfdirection + sex:stimulation + nationality:stimulation + 
##     sex:hedonism + nationality:hedonism + sex:power + nationality:power + 
##     sex:achievement + nationality:achievement + sex:security + 
##     nationality:security + sex:nationality:tradition + sex:nationality:universalism + 
##     sex:nationality:selfdirection + sex:nationality:stimulation + 
##     sex:nationality:hedonism + sex:nationality:power + sex:nationality:achievement + 
##     sex:nationality:security
## 
##                                 Df Sum of Sq    RSS    AIC
## - sex:nationality:universalism   1      6.18 8477.6 550.57
## - sex:nationality:tradition      1     19.36 8490.8 550.74
## - sex:nationality:selfdirection  1     21.91 8493.3 550.77
## - conformity:nationality         1     25.81 8497.2 550.82
## - sex:nationality:stimulation    1     71.57 8543.0 551.41
## - conformity:sex                 2    279.23 8750.6 552.02
## - sex:nationality:hedonism       1    122.93 8594.3 552.06
## - sex:nationality:security       1    139.12 8610.5 552.26
## <none>                                       8471.4 552.49
## - sex:nationality:achievement    1    160.84 8632.2 552.54
## - sex:nationality:power          1    374.19 8845.6 555.20
## - sex:benevolence                2    548.87 9020.3 555.33
## - nationality:benevolence        1    394.13 8865.5 555.45
## 
## Step:  AIC=550.57
## ei ~ conformity:sex + conformity:nationality + sex:tradition + 
##     nationality:tradition + sex:benevolence + nationality:benevolence + 
##     sex:universalism + nationality:universalism + sex:selfdirection + 
##     nationality:selfdirection + sex:stimulation + nationality:stimulation + 
##     sex:hedonism + nationality:hedonism + sex:power + nationality:power + 
##     sex:achievement + nationality:achievement + sex:security + 
##     nationality:security + sex:nationality:tradition + sex:nationality:selfdirection + 
##     sex:nationality:stimulation + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                                 Df Sum of Sq    RSS    AIC
## - nationality:universalism       1      8.58 8486.2 548.68
## - sex:nationality:selfdirection  1     17.22 8494.8 548.79
## - conformity:nationality         1     22.28 8499.9 548.85
## - sex:nationality:tradition      1     32.48 8510.1 548.98
## - sex:nationality:stimulation    1     84.18 8561.8 549.64
## - conformity:sex                 2    278.51 8756.1 550.09
## - sex:nationality:hedonism       1    126.35 8603.9 550.18
## - sex:nationality:security       1    140.82 8618.4 550.36
## - sex:nationality:achievement    1    156.29 8633.9 550.56
## <none>                                       8477.6 550.57
## - sex:benevolence                2    542.69 9020.3 553.33
## - sex:nationality:power          1    382.58 8860.2 553.38
## - nationality:benevolence        1    388.22 8865.8 553.45
## - sex:universalism               2    658.15 9135.7 554.72
## 
## Step:  AIC=548.68
## ei ~ conformity:sex + conformity:nationality + sex:tradition + 
##     nationality:tradition + sex:benevolence + nationality:benevolence + 
##     sex:universalism + sex:selfdirection + nationality:selfdirection + 
##     sex:stimulation + nationality:stimulation + sex:hedonism + 
##     nationality:hedonism + sex:power + nationality:power + sex:achievement + 
##     nationality:achievement + sex:security + nationality:security + 
##     sex:nationality:tradition + sex:nationality:selfdirection + 
##     sex:nationality:stimulation + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                                 Df Sum of Sq    RSS    AIC
## - sex:nationality:selfdirection  1     14.07 8500.2 546.86
## - conformity:nationality         1     21.28 8507.4 546.95
## - sex:nationality:tradition      1     32.03 8518.2 547.09
## - sex:nationality:stimulation    1     83.66 8569.8 547.75
## - conformity:sex                 2    284.75 8770.9 548.28
## - sex:nationality:hedonism       1    128.56 8614.7 548.32
## - sex:nationality:security       1    146.31 8632.5 548.54
## <none>                                       8486.2 548.68
## - sex:nationality:achievement    1    166.05 8652.2 548.79
## - sex:nationality:power          1    376.14 8862.3 551.41
## - sex:benevolence                2    597.61 9083.8 552.10
## - nationality:benevolence        1    521.39 9007.5 553.18
## - sex:universalism               2    889.99 9376.2 555.55
## 
## Step:  AIC=546.86
## ei ~ conformity:sex + conformity:nationality + sex:tradition + 
##     nationality:tradition + sex:benevolence + nationality:benevolence + 
##     sex:universalism + sex:selfdirection + nationality:selfdirection + 
##     sex:stimulation + nationality:stimulation + sex:hedonism + 
##     nationality:hedonism + sex:power + nationality:power + sex:achievement + 
##     nationality:achievement + sex:security + nationality:security + 
##     sex:nationality:tradition + sex:nationality:stimulation + 
##     sex:nationality:hedonism + sex:nationality:power + sex:nationality:achievement + 
##     sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - nationality:selfdirection    1      0.36 8500.6 544.86
## - conformity:nationality       1     18.01 8518.2 545.09
## - sex:nationality:tradition    1     21.98 8522.2 545.14
## - sex:nationality:stimulation  1     74.04 8574.3 545.80
## - conformity:sex               2    270.89 8771.1 546.28
## - sex:nationality:hedonism     1    114.50 8614.7 546.32
## <none>                                     8500.2 546.86
## - sex:nationality:achievement  1    257.85 8758.1 548.12
## - sex:nationality:security     1    260.39 8760.6 548.15
## - sex:selfdirection            2    556.78 9057.0 549.77
## - sex:benevolence              2    583.54 9083.8 550.10
## - sex:nationality:power        1    452.72 8953.0 550.51
## - nationality:benevolence      1    518.21 9018.4 551.31
## - sex:universalism             2    919.86 9420.1 554.06
## 
## Step:  AIC=544.86
## ei ~ conformity:sex + conformity:nationality + sex:tradition + 
##     nationality:tradition + sex:benevolence + nationality:benevolence + 
##     sex:universalism + sex:selfdirection + sex:stimulation + 
##     nationality:stimulation + sex:hedonism + nationality:hedonism + 
##     sex:power + nationality:power + sex:achievement + nationality:achievement + 
##     sex:security + nationality:security + sex:nationality:tradition + 
##     sex:nationality:stimulation + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - conformity:nationality       1     18.24 8518.8 543.10
## - sex:nationality:tradition    1     22.43 8523.0 543.15
## - sex:nationality:stimulation  1     73.80 8574.4 543.81
## - sex:nationality:hedonism     1    114.43 8615.0 544.32
## - conformity:sex               2    275.81 8776.4 544.34
## <none>                                     8500.6 544.86
## - sex:nationality:achievement  1    262.63 8763.2 546.18
## - sex:nationality:security     1    264.14 8764.7 546.20
## - sex:benevolence              2    621.43 9122.0 548.55
## - sex:nationality:power        1    483.47 8984.1 548.89
## - sex:selfdirection            2    713.67 9214.3 549.65
## - nationality:benevolence      1    589.09 9089.7 550.17
## - sex:universalism             2    921.57 9422.2 552.08
## 
## Step:  AIC=543.1
## ei ~ conformity:sex + sex:tradition + nationality:tradition + 
##     sex:benevolence + nationality:benevolence + sex:universalism + 
##     sex:selfdirection + sex:stimulation + nationality:stimulation + 
##     sex:hedonism + nationality:hedonism + sex:power + nationality:power + 
##     sex:achievement + nationality:achievement + sex:security + 
##     nationality:security + sex:nationality:tradition + sex:nationality:stimulation + 
##     sex:nationality:hedonism + sex:nationality:power + sex:nationality:achievement + 
##     sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - sex:tradition:nationality    1     21.74 8540.6 541.37
## - sex:nationality:stimulation  1     74.54 8593.4 542.05
## - sex:nationality:hedonism     1    121.54 8640.4 542.64
## - conformity:sex               2    283.01 8801.9 542.66
## <none>                                     8518.8 543.10
## - sex:nationality:achievement  1    267.25 8786.1 544.46
## - sex:nationality:security     1    274.23 8793.1 544.55
## - sex:benevolence              2    632.71 9151.6 546.91
## - sex:nationality:power        1    500.11 9019.0 547.32
## - sex:selfdirection            2    695.79 9214.6 547.65
## - nationality:benevolence      1    648.87 9167.7 549.10
## - sex:universalism             2    960.48 9479.3 550.74
## 
## Step:  AIC=541.37
## ei ~ conformity:sex + sex:tradition + tradition:nationality + 
##     sex:benevolence + nationality:benevolence + sex:universalism + 
##     sex:selfdirection + sex:stimulation + nationality:stimulation + 
##     sex:hedonism + nationality:hedonism + sex:power + nationality:power + 
##     sex:achievement + nationality:achievement + sex:security + 
##     nationality:security + sex:nationality:stimulation + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - sex:tradition                2      8.85 8549.4 537.49
## - tradition:nationality        1     52.76 8593.3 540.05
## - sex:nationality:stimulation  1     75.82 8616.4 540.34
## - sex:nationality:hedonism     1    109.35 8649.9 540.76
## - conformity:sex               2    290.53 8831.1 541.02
## <none>                                     8540.6 541.37
## - sex:nationality:achievement  1    250.12 8790.7 542.52
## - sex:nationality:security     1    257.24 8797.8 542.61
## - sex:benevolence              2    646.18 9186.8 545.32
## - sex:nationality:power        1    512.45 9053.0 545.73
## - sex:selfdirection            2    692.99 9233.6 545.88
## - nationality:benevolence      1    637.73 9178.3 547.22
## - sex:universalism             2    952.88 9493.5 548.90
## 
## Step:  AIC=539.49
## ei ~ conformity:sex + tradition:nationality + sex:benevolence + 
##     nationality:benevolence + sex:universalism + sex:selfdirection + 
##     sex:stimulation + nationality:stimulation + sex:hedonism + 
##     nationality:hedonism + sex:power + nationality:power + sex:achievement + 
##     nationality:achievement + sex:security + nationality:security + 
##     sex:nationality:stimulation + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - tradition:nationality        2     82.37 8631.6 536.53
## - sex:nationality:stimulation  1     77.00 8626.2 538.46
## - sex:nationality:hedonism     1    104.53 8653.8 538.81
## <none>                                     8549.2 539.49
## - conformity:sex               2    377.05 8926.3 540.19
## - sex:nationality:security     1    248.66 8797.9 540.61
## - sex:nationality:achievement  1    252.02 8801.3 540.65
## - sex:benevolence              2    651.78 9201.0 543.49
## - sex:nationality:power        1    504.02 9053.3 543.73
## - sex:selfdirection            2    708.01 9257.2 544.16
## - nationality:benevolence      1    638.07 9187.3 545.33
## - sex:universalism             2    978.48 9527.7 547.30
## 
## Step:  AIC=536.53
## ei ~ conformity:sex + sex:benevolence + nationality:benevolence + 
##     sex:universalism + sex:selfdirection + sex:stimulation + 
##     nationality:stimulation + sex:hedonism + nationality:hedonism + 
##     sex:power + nationality:power + sex:achievement + nationality:achievement + 
##     sex:security + nationality:security + sex:nationality:stimulation + 
##     sex:nationality:hedonism + sex:nationality:power + sex:nationality:achievement + 
##     sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - sex:nationality:stimulation  1     62.98 8694.6 535.32
## - sex:nationality:hedonism     1    111.09 8742.7 535.92
## <none>                                     8631.6 536.53
## - sex:nationality:security     1    243.02 8874.6 537.56
## - sex:nationality:achievement  1    258.98 8890.6 537.75
## - conformity:sex               2    493.08 9124.7 538.59
## - sex:benevolence              2    624.17 9255.8 540.14
## - sex:nationality:power        1    529.63 9161.2 541.02
## - sex:selfdirection            2    721.65 9353.3 541.28
## - benevolence:nationality      1    598.69 9230.3 541.84
## - sex:universalism             2   1080.50 9712.1 545.39
## 
## Step:  AIC=535.32
## ei ~ conformity:sex + sex:benevolence + benevolence:nationality + 
##     sex:universalism + sex:selfdirection + sex:stimulation + 
##     nationality:stimulation + sex:hedonism + nationality:hedonism + 
##     sex:power + nationality:power + sex:achievement + nationality:achievement + 
##     sex:security + nationality:security + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - sex:stimulation              2     37.19 8731.8 531.79
## - nationality:stimulation      1     56.38 8751.0 534.03
## <none>                                     8694.6 535.32
## - sex:nationality:security     1    189.75 8884.3 535.68
## - sex:nationality:hedonism     1    193.36 8887.9 535.72
## - sex:nationality:achievement  1    236.76 8931.3 536.25
## - conformity:sex               2    487.45 9182.0 537.27
## - sex:benevolence              2    687.82 9382.4 539.62
## - sex:selfdirection            2    729.54 9424.1 540.11
## - benevolence:nationality      1    604.96 9299.5 540.65
## - sex:nationality:power        1    689.59 9384.2 541.64
## - sex:universalism             2   1085.39 9780.0 544.15
## 
## Step:  AIC=533.72
## ei ~ conformity:sex + sex:benevolence + benevolence:nationality + 
##     sex:universalism + sex:selfdirection + nationality:stimulation + 
##     sex:hedonism + nationality:hedonism + sex:power + nationality:power + 
##     sex:achievement + nationality:achievement + sex:security + 
##     nationality:security + sex:nationality:hedonism + sex:nationality:power + 
##     sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## - nationality:stimulation      2    105.09 8831.2 531.02
## <none>                                     8726.1 533.72
## - sex:nationality:hedonism     1    181.38 8907.5 533.96
## - sex:nationality:security     1    182.01 8908.1 533.97
## - sex:nationality:achievement  1    214.75 8940.9 534.37
## - conformity:sex               2    462.33 9188.4 535.34
## - sex:selfdirection            2    701.28 9427.4 538.14
## - sex:benevolence              2    706.08 9432.2 538.20
## - benevolence:nationality      1    645.25 9371.3 539.49
## - sex:nationality:power        1    660.29 9386.4 539.67
## - sex:universalism             2   1123.64 9849.7 542.92
## 
## Step:  AIC=531.02
## ei ~ conformity:sex + sex:benevolence + benevolence:nationality + 
##     sex:universalism + sex:selfdirection + sex:hedonism + nationality:hedonism + 
##     sex:power + nationality:power + sex:achievement + nationality:achievement + 
##     sex:security + nationality:security + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security
## 
##                               Df Sum of Sq    RSS    AIC
## <none>                                     8831.2 531.02
## - sex:nationality:security     1    184.03 9015.2 531.27
## - sex:nationality:hedonism     1    241.78 9073.0 531.97
## - sex:nationality:achievement  1    250.56 9081.8 532.07
## - conformity:sex               2    421.00 9252.2 532.10
## - sex:selfdirection            2    632.77 9464.0 534.57
## - sex:benevolence              2    662.51 9493.7 534.91
## - benevolence:nationality      1    592.32 9423.5 536.10
## - sex:nationality:power        1    678.22 9509.4 537.09
## - sex:universalism             2   1156.73 9987.9 540.44
summary(final_model)
## 
## Call:
## lm(formula = ei ~ conformity:sex + sex:benevolence + benevolence:nationality + 
##     sex:universalism + sex:selfdirection + sex:hedonism + nationality:hedonism + 
##     sex:power + nationality:power + sex:achievement + nationality:achievement + 
##     sex:security + nationality:security + sex:nationality:hedonism + 
##     sex:nationality:power + sex:nationality:achievement + sex:nationality:security, 
##     data = thes)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -28.7423  -5.5827   0.8924   5.6468  19.2840 
## 
## Coefficients:
##                                 Estimate Std. Error t value  Pr(>|t|)    
## (Intercept)                     30.01546    7.10179   4.226 0.0000606 ***
## conformity:sexf                  4.30588    2.16468   1.989   0.04998 *  
## conformity:sexm                  0.01795    2.59563   0.007   0.99450    
## sexf:benevolence                 4.43775    2.09821   2.115   0.03743 *  
## sexm:benevolence                 7.20730    3.80333   1.895   0.06157 .  
## benevolence:nationalityrus      -6.40013    2.71257  -2.359   0.02065 *  
## sexf:universalism                5.33452    1.65299   3.227   0.00179 ** 
## sexm:universalism                1.88404    3.06189   0.615   0.54003    
## sexf:selfdirection              -4.50260    1.84643  -2.439   0.01688 *  
## sexm:selfdirection              -0.26169    2.62434  -0.100   0.92081    
## sexf:hedonism                    4.59747    1.98498   2.316   0.02302 *  
## sexm:hedonism                   -4.31232    2.37819  -1.813   0.07340 .  
## nationalityrus:hedonism         -0.66022    2.78702  -0.237   0.81333    
## sexf:power                       0.32491    2.16378   0.150   0.88101    
## sexm:power                      -1.25999    1.94852  -0.647   0.51965    
## nationalityrus:power            -5.97135    2.95341  -2.022   0.04641 *  
## sexf:achievement                -4.23441    2.31363  -1.830   0.07081 .  
## sexm:achievement                -1.10047    2.53070  -0.435   0.66480    
## nationalityrus:achievement       8.57061    3.43127   2.498   0.01447 *  
## sexf:security                   -3.78674    2.28997  -1.654   0.10198    
## sexm:security                    2.47945    1.99070   1.246   0.21645    
## nationalityrus:security          5.02950    2.92917   1.717   0.08970 .  
## sexm:nationalityrus:hedonism     5.80710    3.85230   1.507   0.13550    
## sexm:nationalityrus:power        9.41188    3.72787   2.525   0.01348 *  
## sexm:nationalityrus:achievement -6.73244    4.38717  -1.535   0.12869    
## sexm:nationalityrus:security    -6.16767    4.68978  -1.315   0.19209    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.32 on 83 degrees of freedom
## Multiple R-squared:  0.4873, Adjusted R-squared:  0.3329 
## F-statistic: 3.156 on 25 and 83 DF,  p-value: 0.00004648
model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security) : (nationality+sex), data = thes)
final_model <- step(model)  
## Start:  AIC=548.54
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security):(nationality + 
##     sex)
## 
##                             Df Sum of Sq     RSS    AIC
## - tradition:nationality      2     64.82  9526.5 545.28
## - stimulation:nationality    2     92.59  9554.2 545.60
## - security:nationality       2    111.68  9573.3 545.82
## - tradition:sex              1      3.29  9465.0 546.58
## - stimulation:sex            1      4.72  9466.4 546.59
## - achievement:sex            1      9.74  9471.4 546.65
## - conformity:nationality     2    224.06  9685.7 547.09
## - achievement:nationality    2    287.97  9749.6 547.81
## - selfdirection:sex          1    150.02  9611.7 548.25
## - security:sex               1    150.12  9611.8 548.25
## <none>                                    9461.7 548.54
## - benevolence:sex            1    177.01  9638.7 548.56
## - universalism:sex           1    209.68  9671.3 548.93
## - conformity:sex             1    227.56  9689.2 549.13
## - power:sex                  1    298.90  9760.6 549.93
## - selfdirection:nationality  2    572.76 10034.4 550.94
## - benevolence:nationality    2    605.21 10066.9 551.30
## - power:nationality          2    606.28 10067.9 551.31
## - universalism:nationality   2    776.18 10237.8 553.13
## - hedonism:nationality       2    817.63 10279.3 553.57
## - hedonism:sex               1    691.09 10152.7 554.22
## 
## Step:  AIC=546.97
## ei ~ conformity:nationality + conformity:sex + tradition:sex + 
##     benevolence:nationality + benevolence:sex + universalism:nationality + 
##     universalism:sex + selfdirection:nationality + selfdirection:sex + 
##     stimulation:nationality + stimulation:sex + hedonism:nationality + 
##     hedonism:sex + power:nationality + power:sex + achievement:nationality + 
##     achievement:sex + security:nationality + security:sex
## 
##                             Df Sum of Sq     RSS    AIC
## - sex:tradition              2     70.87  9570.0 543.78
## - nationality:security       2     93.06  9592.2 544.03
## - nationality:stimulation    2     93.40  9592.5 544.04
## - sex:stimulation            1      2.70  9501.8 545.00
## - sex:achievement            1      5.40  9504.5 545.03
## - conformity:nationality     2    202.08  9701.2 545.26
## - nationality:achievement    2    263.40  9762.5 545.95
## - sex:selfdirection          1    138.27  9637.4 546.54
## - sex:security               1    141.08  9640.2 546.58
## - sex:benevolence            1    172.86  9672.0 546.93
## <none>                                    9499.1 546.97
## - sex:universalism           1    198.85  9698.0 547.23
## - conformity:sex             1    240.81  9739.9 547.70
## - sex:power                  1    320.57  9819.7 548.59
## - nationality:selfdirection  2    548.58 10047.7 549.09
## - nationality:benevolence    2    582.91 10082.0 549.46
## - nationality:power          2    584.59 10083.7 549.48
## - nationality:universalism   2    783.76 10282.9 551.61
## - nationality:hedonism       2    788.04 10287.2 551.66
## - sex:hedonism               1    666.42 10165.5 552.36
## 
## Step:  AIC=543.78
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:stimulation + 
##     sex:stimulation + nationality:hedonism + sex:hedonism + nationality:power + 
##     sex:power + nationality:achievement + sex:achievement + nationality:security + 
##     sex:security
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:security       2     74.05  9644.1 540.62
## - nationality:stimulation    2     95.64  9665.6 540.86
## - sex:stimulation            1      2.32  9572.3 541.81
## - sex:achievement            1     15.87  9585.9 541.96
## - nationality:achievement    2    278.23  9848.2 542.90
## - sex:security               1    149.31  9719.3 543.47
## <none>                                    9570.0 543.78
## - sex:benevolence            1    183.86  9753.9 543.85
## - conformity:nationality     2    369.70  9939.7 543.91
## - sex:selfdirection          1    193.10  9763.1 543.96
## - sex:universalism           1    207.45  9777.5 544.12
## - conformity:sex             1    298.73  9868.7 545.13
## - sex:power                  1    337.29  9907.3 545.55
## - nationality:power          2    555.50 10125.5 545.93
## - nationality:selfdirection  2    596.94 10166.9 546.37
## - nationality:benevolence    2    667.71 10237.7 547.13
## - nationality:hedonism       2    784.97 10355.0 548.37
## - nationality:universalism   2    906.70 10476.7 549.65
## - sex:hedonism               1    744.57 10314.6 549.95
## 
## Step:  AIC=542.6
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:stimulation + 
##     sex:stimulation + nationality:hedonism + sex:hedonism + nationality:power + 
##     sex:power + nationality:achievement + sex:achievement + sex:security
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:stimulation    2    106.18  9748.6 539.80
## - sex:stimulation            1      1.57  9644.0 540.62
## - sex:security               2    197.75  9840.2 540.81
## - sex:achievement            1     24.16  9666.6 540.87
## - nationality:achievement    2    247.92  9890.3 541.37
## - sex:benevolence            1    170.20  9812.6 542.51
## - sex:universalism           1    173.26  9815.7 542.54
## <none>                                    9642.4 542.60
## - sex:selfdirection          1    199.64  9842.1 542.84
## - conformity:nationality     2    384.37 10026.8 542.86
## - conformity:sex             1    297.14  9939.6 543.91
## - nationality:selfdirection  2    572.93 10215.4 544.89
## - nationality:benevolence    2    666.38 10308.8 545.89
## - sex:power                  1    510.89 10153.3 546.23
## - nationality:power          2    732.37 10374.8 546.58
## - nationality:hedonism       2    800.98 10443.4 547.30
## - sex:hedonism               1    713.89 10356.3 548.39
## - nationality:universalism   2    979.01 10621.4 549.14
## 
## Step:  AIC=541.79
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + sex:stimulation + 
##     nationality:hedonism + sex:hedonism + nationality:power + 
##     sex:power + nationality:achievement + sex:achievement + sex:security
## 
##                             Df Sum of Sq     RSS    AIC
## - sex:stimulation            2      1.03  9748.9 537.80
## - sex:security               2    163.50  9911.4 539.60
## - sex:achievement            1     27.99  9775.9 540.10
## - nationality:achievement    2    267.76 10015.7 540.74
## <none>                                    9747.9 541.79
## - sex:universalism           1    189.77  9937.7 541.89
## - sex:selfdirection          1    192.66  9940.6 541.92
## - sex:benevolence            1    193.73  9941.6 541.93
## - conformity:nationality     2    421.88 10169.8 542.41
## - conformity:sex             1    272.70 10020.6 542.79
## - nationality:selfdirection  2    579.29 10327.2 544.08
## - nationality:benevolence    2    641.38 10389.3 544.73
## - nationality:power          2    670.82 10418.7 545.04
## - sex:power                  1    493.93 10241.8 545.17
## - nationality:hedonism       2    769.25 10517.2 546.07
## - sex:hedonism               1    631.05 10379.0 546.62
## - nationality:universalism   2    952.43 10700.3 547.95
## 
## Step:  AIC=537.8
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:hedonism + 
##     sex:hedonism + nationality:power + sex:power + nationality:achievement + 
##     sex:achievement + sex:security
## 
##                             Df Sum of Sq     RSS    AIC
## - sex:security               2    180.56  9929.5 535.80
## - sex:achievement            1     28.69  9777.6 536.12
## - nationality:achievement    2    276.35 10025.3 536.85
## <none>                                    9748.9 537.80
## - sex:selfdirection          1    193.62  9942.6 537.94
## - sex:benevolence            1    194.53  9943.5 537.95
## - sex:universalism           1    202.20  9951.1 538.04
## - conformity:nationality     2    425.91 10174.9 538.46
## - conformity:sex             1    292.17 10041.1 539.02
## - nationality:selfdirection  2    591.15 10340.1 540.22
## - nationality:benevolence    2    650.21 10399.2 540.84
## - nationality:power          2    717.92 10466.9 541.54
## - sex:power                  1    541.97 10290.9 541.70
## - nationality:hedonism       2    770.51 10519.5 542.09
## - sex:hedonism               1    663.60 10412.5 542.98
## - nationality:universalism   2    998.17 10747.1 544.42
## 
## Step:  AIC=535.8
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:hedonism + 
##     sex:hedonism + nationality:power + sex:power + nationality:achievement + 
##     sex:achievement
## 
##                             Df Sum of Sq     RSS    AIC
## - sex:achievement            1     35.07  9964.6 534.18
## - nationality:achievement    2    301.41 10230.9 535.06
## <none>                                    9929.5 535.80
## - sex:universalism           1    212.00 10141.5 536.10
## - conformity:sex             1    236.75 10166.3 536.37
## - sex:benevolence            1    257.81 10187.3 536.59
## - sex:selfdirection          1    271.73 10201.2 536.74
## - conformity:nationality     2    493.16 10422.7 537.08
## - nationality:selfdirection  2    575.80 10505.3 537.94
## - nationality:power          2    760.28 10689.8 539.84
## - sex:power                  1    623.86 10553.4 540.44
## - nationality:benevolence    2    841.78 10771.3 540.67
## - sex:hedonism               1    672.67 10602.2 540.94
## - nationality:hedonism       2    887.66 10817.2 541.13
## - nationality:universalism   2   1126.14 11055.6 543.51
## 
## Step:  AIC=534.18
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:hedonism + 
##     sex:hedonism + nationality:power + sex:power + nationality:achievement
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:achievement    2    267.46 10232.0 533.07
## - sex:universalism           1    182.78 10147.4 534.16
## <none>                                    9964.6 534.18
## - sex:benevolence            1    222.76 10187.3 534.59
## - sex:selfdirection          1    238.14 10202.7 534.76
## - conformity:sex             1    278.62 10243.2 535.19
## - conformity:nationality     2    519.52 10484.1 535.72
## - nationality:selfdirection  2    543.68 10508.3 535.97
## - nationality:power          2    741.17 10705.8 538.00
## - nationality:benevolence    2    806.70 10771.3 538.67
## - sex:power                  1    631.34 10595.9 538.88
## - nationality:hedonism       2    866.37 10830.9 539.27
## - sex:hedonism               1    674.97 10639.6 539.33
## - nationality:universalism   2   1092.05 11056.6 541.52
## 
## Step:  AIC=533.07
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + sex:universalism + 
##     nationality:selfdirection + sex:selfdirection + nationality:hedonism + 
##     sex:hedonism + nationality:power + sex:power
## 
##                             Df Sum of Sq   RSS    AIC
## - sex:universalism           1     99.93 10332 532.13
## - sex:benevolence            1    123.82 10356 532.38
## <none>                                   10232 533.07
## - sex:selfdirection          1    198.01 10430 533.16
## - conformity:sex             1    296.49 10528 534.18
## - conformity:nationality     2    571.00 10803 534.99
## - nationality:selfdirection  2    571.87 10804 535.00
## - nationality:benevolence    2    601.49 10834 535.30
## - sex:hedonism               1    509.99 10742 536.37
## - nationality:hedonism       2    717.64 10950 536.46
## - sex:power                  1    570.58 10803 536.99
## - nationality:power          2    941.48 11174 538.67
## - nationality:universalism   2    983.59 11216 539.08
## 
## Step:  AIC=532.13
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     sex:benevolence + nationality:universalism + nationality:selfdirection + 
##     sex:selfdirection + nationality:hedonism + sex:hedonism + 
##     nationality:power + sex:power
## 
##                             Df Sum of Sq   RSS    AIC
## - sex:benevolence            1     61.29 10393 530.77
## - sex:selfdirection          1    122.01 10454 531.41
## <none>                                   10332 532.13
## - nationality:selfdirection  2    510.89 10843 533.39
## - nationality:benevolence    2    589.87 10922 534.18
## - conformity:sex             1    394.08 10726 534.21
## - conformity:nationality     2    624.68 10957 534.53
## - sex:hedonism               1    557.36 10889 535.86
## - nationality:hedonism       2    819.87 11152 536.45
## - sex:power                  1    654.57 10986 536.83
## - nationality:universalism   2    932.60 11265 537.55
## - nationality:power          2   1108.05 11440 539.23
## 
## Step:  AIC=530.77
## ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + sex:selfdirection + 
##     nationality:hedonism + sex:hedonism + nationality:power + 
##     sex:power
## 
##                             Df Sum of Sq   RSS    AIC
## <none>                                   10393 530.77
## - conformity:sex             1    339.59 10733 532.28
## - conformity:nationality     2    569.04 10962 532.58
## - sex:selfdirection          1    382.54 10776 532.71
## - nationality:selfdirection  2    623.75 11017 533.13
## - nationality:benevolence    2    662.68 11056 533.51
## - sex:hedonism               1    496.34 10890 533.86
## - nationality:hedonism       2    768.24 11162 534.55
## - sex:power                  1    593.45 10987 534.83
## - nationality:universalism   2    957.38 11351 536.38
## - nationality:power          2   1102.95 11496 537.77
summary(final_model)
## 
## Call:
## lm(formula = ei ~ conformity:nationality + conformity:sex + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + sex:selfdirection + 
##     nationality:hedonism + sex:hedonism + nationality:power + 
##     sex:power, data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -26.841  -7.597   0.702   6.428  23.279 
## 
## Coefficients:
##                              Estimate Std. Error t value  Pr(>|t|)    
## (Intercept)                    29.180      6.900   4.229 0.0000553 ***
## conformity:nationalityfr        4.632      2.454   1.888    0.0622 .  
## conformity:nationalityrus       3.740      2.336   1.601    0.1129    
## conformity:sexm                -4.316      2.489  -1.734    0.0863 .  
## nationalityfr:benevolence       4.891      2.071   2.362    0.0203 *  
## nationalityrus:benevolence     -1.004      2.181  -0.460    0.6463    
## nationalityfr:universalism      3.512      1.812   1.939    0.0556 .  
## nationalityrus:universalism     4.534      1.998   2.269    0.0256 *  
## nationalityfr:selfdirection    -4.912      2.112  -2.326    0.0222 *  
## nationalityrus:selfdirection   -1.806      2.317  -0.780    0.4376    
## sexm:selfdirection              3.982      2.164   1.840    0.0690 .  
## nationalityfr:hedonism          2.015      1.486   1.356    0.1785    
## nationalityrus:hedonism         4.315      1.699   2.539    0.0128 *  
## sexm:hedonism                  -3.908      1.865  -2.096    0.0388 *  
## nationalityfr:power            -4.453      1.466  -3.037    0.0031 ** 
## nationalityrus:power           -2.534      1.451  -1.746    0.0841 .  
## sexm:power                      3.658      1.596   2.292    0.0242 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.63 on 92 degrees of freedom
## Multiple R-squared:  0.3966, Adjusted R-squared:  0.2917 
## F-statistic:  3.78 on 16 and 92 DF,  p-value: 0.00002452

env identity vs (values+well-being):sex:nationality

model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security + positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance) :sex, data = thes)
final_model <- step(model)  
## Start:  AIC=557.35
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security + 
##     positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance):sex
## 
##                      Df Sum of Sq     RSS    AIC
## - selfacceptance:sex  2     36.72  9925.8 553.76
## - stimulation:sex     2     51.39  9940.5 553.92
## - benevolence:sex     2     51.54  9940.6 553.92
## - positiverelat:sex   2     79.88  9969.0 554.23
## - achievement:sex     2    104.95  9994.0 554.51
## - purposeinlife:sex   2    124.98 10014.1 554.72
## - tradition:sex       2    144.52 10033.6 554.94
## - envmattery:sex      2    148.02 10037.1 554.97
## - autonomy:sex        2    187.91 10077.0 555.41
## - security:sex        2    213.78 10102.9 555.69
## - conformity:sex      2    310.48 10199.6 556.72
## <none>                             9889.1 557.35
## - hedonism:sex        2    400.22 10289.3 557.68
## - selfdirection:sex   2    405.04 10294.1 557.73
## - universalism:sex    2    609.91 10499.0 559.88
## - power:sex           2    734.84 10623.9 561.17
## - persgrowth:sex      2   1206.04 11095.1 565.90
## 
## Step:  AIC=553.76
## ei ~ conformity:sex + tradition:sex + benevolence:sex + universalism:sex + 
##     selfdirection:sex + stimulation:sex + hedonism:sex + power:sex + 
##     achievement:sex + security:sex + positiverelat:sex + autonomy:sex + 
##     envmattery:sex + persgrowth:sex + purposeinlife:sex
## 
##                     Df Sum of Sq     RSS    AIC
## - sex:stimulation    2     52.06  9977.9 550.33
## - sex:benevolence    2     62.92  9988.7 550.45
## - sex:achievement    2     74.11  9999.9 550.57
## - sex:positiverelat  2     77.38 10003.2 550.61
## - sex:envmattery     2    128.15 10054.0 551.16
## - sex:tradition      2    129.18 10055.0 551.17
## - sex:purposeinlife  2    159.60 10085.4 551.50
## - sex:autonomy       2    160.06 10085.9 551.50
## - sex:security       2    209.84 10135.7 552.04
## - conformity:sex     2    301.77 10227.6 553.02
## <none>                            9925.8 553.76
## - sex:selfdirection  2    397.47 10323.3 554.04
## - sex:hedonism       2    436.52 10362.3 554.45
## - sex:universalism   2    622.27 10548.1 556.39
## - sex:power          2    721.31 10647.1 557.41
## - sex:persgrowth     2   1181.30 11107.1 562.02
## 
## Step:  AIC=550.33
## ei ~ conformity:sex + sex:tradition + sex:benevolence + sex:universalism + 
##     sex:selfdirection + sex:hedonism + sex:power + sex:achievement + 
##     sex:security + sex:positiverelat + sex:autonomy + sex:envmattery + 
##     sex:persgrowth + sex:purposeinlife
## 
##                     Df Sum of Sq     RSS    AIC
## - sex:benevolence    2     59.62 10037.5 546.98
## - sex:positiverelat  2     77.49 10055.4 547.17
## - sex:achievement    2     82.86 10060.7 547.23
## - sex:envmattery     2     93.50 10071.4 547.35
## - sex:tradition      2    137.39 10115.3 547.82
## - sex:autonomy       2    143.73 10121.6 547.89
## - sex:purposeinlife  2    143.99 10121.9 547.89
## - sex:security       2    174.43 10152.3 548.22
## - conformity:sex     2    301.58 10279.5 549.57
## <none>                            9977.9 550.33
## - sex:selfdirection  2    408.00 10385.9 550.70
## - sex:hedonism       2    434.26 10412.1 550.97
## - sex:universalism   2    601.51 10579.4 552.71
## - sex:power          2    832.22 10810.1 555.06
## - sex:persgrowth     2   1135.73 11113.6 558.08
## 
## Step:  AIC=546.98
## ei ~ conformity:sex + sex:tradition + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:achievement + sex:security + 
##     sex:positiverelat + sex:autonomy + sex:envmattery + sex:persgrowth + 
##     sex:purposeinlife
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:positiverelat  2     83.12 10121 543.88
## - sex:envmattery     2     99.69 10137 544.06
## - sex:achievement    2    105.28 10143 544.12
## - sex:tradition      2    148.98 10186 544.58
## - sex:purposeinlife  2    160.41 10198 544.71
## - sex:security       2    181.92 10219 544.94
## - sex:autonomy       2    193.40 10231 545.06
## - conformity:sex     2    359.11 10397 546.81
## - sex:selfdirection  2    365.31 10403 546.87
## <none>                           10038 546.98
## - sex:hedonism       2    438.15 10476 547.64
## - sex:universalism   2    774.41 10812 551.08
## - sex:power          2    906.77 10944 552.41
## - sex:persgrowth     2   1214.58 11252 555.43
## 
## Step:  AIC=543.88
## ei ~ conformity:sex + sex:tradition + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:achievement + sex:security + 
##     sex:autonomy + sex:envmattery + sex:persgrowth + sex:purposeinlife
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:achievement    2     95.02 10216 540.90
## - sex:envmattery     2    112.84 10234 541.09
## - sex:tradition      2    167.08 10288 541.66
## - sex:purposeinlife  2    184.93 10306 541.85
## - sex:autonomy       2    199.58 10320 542.01
## - sex:security       2    229.81 10350 542.32
## - sex:selfdirection  2    373.06 10494 543.82
## <none>                           10121 543.88
## - conformity:sex     2    457.92 10578 544.70
## - sex:hedonism       2    509.32 10630 545.23
## - sex:universalism   2    786.07 10907 548.03
## - sex:power          2    851.45 10972 548.68
## - sex:persgrowth     2   1246.61 11367 552.54
## 
## Step:  AIC=540.9
## ei ~ conformity:sex + sex:tradition + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security + sex:autonomy + 
##     sex:envmattery + sex:persgrowth + sex:purposeinlife
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:envmattery     2    101.76 10317 537.98
## - sex:tradition      2    131.94 10348 538.29
## - sex:purposeinlife  2    184.40 10400 538.85
## - sex:autonomy       2    237.86 10454 539.40
## - sex:security       2    238.33 10454 539.41
## - sex:selfdirection  2    363.17 10579 540.70
## <none>                           10216 540.90
## - conformity:sex     2    465.85 10682 541.76
## - sex:hedonism       2    538.32 10754 542.49
## - sex:universalism   2    758.78 10974 544.71
## - sex:power          2   1037.94 11254 547.44
## - sex:persgrowth     2   1296.37 11512 549.92
## 
## Step:  AIC=537.98
## ei ~ conformity:sex + sex:tradition + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security + sex:autonomy + 
##     sex:persgrowth + sex:purposeinlife
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:tradition      2    122.06 10439 535.26
## - sex:purposeinlife  2    237.23 10555 536.45
## - sex:autonomy       2    290.00 10607 537.00
## - sex:security       2    332.20 10650 537.43
## - sex:selfdirection  2    383.81 10701 537.96
## <none>                           10317 537.98
## - conformity:sex     2    518.25 10836 539.32
## - sex:hedonism       2    572.37 10890 539.86
## - sex:universalism   2    841.00 11158 542.52
## - sex:power          2   1043.68 11361 544.48
## - sex:persgrowth     2   1423.23 11741 548.06
## 
## Step:  AIC=535.26
## ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security + sex:autonomy + 
##     sex:persgrowth + sex:purposeinlife
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:purposeinlife  2    229.88 10669 533.63
## - sex:autonomy       2    355.54 10795 534.91
## - sex:security       2    361.80 10801 534.97
## <none>                           10439 535.26
## - sex:selfdirection  2    398.11 10838 535.34
## - sex:hedonism       2    574.65 11014 537.10
## - conformity:sex     2    862.50 11302 539.91
## - sex:power          2    962.52 11402 540.87
## - sex:universalism   2   1079.17 11519 541.98
## - sex:persgrowth     2   1323.27 11763 544.27
## 
## Step:  AIC=533.63
## ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:security + sex:autonomy + 
##     sex:persgrowth
## 
##                     Df Sum of Sq   RSS    AIC
## - sex:security       2    356.68 11026 533.22
## - sex:autonomy       2    359.10 11028 533.24
## <none>                           10669 533.63
## - sex:hedonism       2    563.98 11233 535.25
## - sex:selfdirection  2    644.78 11314 536.03
## - conformity:sex     2    851.91 11521 538.01
## - sex:power          2    876.97 11546 538.24
## - sex:universalism   2   1141.33 11811 540.71
## - sex:persgrowth     2   1240.34 11910 541.62
## 
## Step:  AIC=533.22
## ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:autonomy + sex:persgrowth
## 
##                     Df Sum of Sq   RSS    AIC
## <none>                           11026 533.22
## - sex:autonomy       2    418.54 11445 533.28
## - sex:hedonism       2    498.54 11525 534.04
## - sex:selfdirection  2    696.01 11722 535.89
## - conformity:sex     2    838.68 11865 537.21
## - sex:power          2    965.22 11991 538.36
## - sex:persgrowth     2   1136.73 12163 539.91
## - sex:universalism   2   1177.30 12203 540.27
summary(final_model)
## 
## Call:
## lm(formula = ei ~ conformity:sex + sex:universalism + sex:selfdirection + 
##     sex:hedonism + sex:power + sex:autonomy + sex:persgrowth, 
##     data = thes)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -29.8979  -5.5965   0.7302   5.6917  21.2211 
## 
## Coefficients:
##                    Estimate Std. Error t value Pr(>|t|)   
## (Intercept)        24.10787   11.30747   2.132  0.03561 * 
## conformity:sexf     5.18113    1.94406   2.665  0.00906 **
## conformity:sexm     0.48205    1.94277   0.248  0.80458   
## sexf:universalism   4.52020    1.57467   2.871  0.00506 **
## sexm:universalism   3.18055    2.31792   1.372  0.17328   
## sexf:selfdirection -4.33648    1.89333  -2.290  0.02423 * 
## sexm:selfdirection  2.04994    2.48040   0.826  0.41064   
## sexf:hedonism       2.63409    1.31325   2.006  0.04775 * 
## sexm:hedonism      -0.41653    1.35291  -0.308  0.75886   
## sexf:power         -3.64833    1.27201  -2.868  0.00510 **
## sexm:power          0.01077    1.03571   0.010  0.99172   
## sexf:autonomy      -0.14397    0.22811  -0.631  0.52948   
## sexm:autonomy      -0.42056    0.23041  -1.825  0.07114 . 
## sexf:persgrowth     0.44903    0.21060   2.132  0.03560 * 
## sexm:persgrowth     0.54084    0.22619   2.391  0.01879 * 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.83 on 94 degrees of freedom
## Multiple R-squared:  0.3599, Adjusted R-squared:  0.2646 
## F-statistic: 3.775 on 14 and 94 DF,  p-value: 0.00004954
model <- lm(ei ~ (conformity + tradition + benevolence + universalism + selfdirection + stimulation + hedonism + power + achievement + security + positiverelat + autonomy + envmattery + persgrowth + purposeinlife + selfacceptance) :nationality, data = thes)
final_model <- step(model)  
## Start:  AIC=554.74
## ei ~ (conformity + tradition + benevolence + universalism + selfdirection + 
##     stimulation + hedonism + power + achievement + security + 
##     positiverelat + autonomy + envmattery + persgrowth + purposeinlife + 
##     selfacceptance):nationality
## 
##                              Df Sum of Sq     RSS    AIC
## - achievement:nationality     2      1.17  9655.6 550.75
## - selfacceptance:nationality  2     32.88  9687.3 551.11
## - positiverelat:nationality   2     38.50  9692.9 551.17
## - stimulation:nationality     2     50.08  9704.5 551.30
## - tradition:nationality       2    111.03  9765.4 551.98
## - security:nationality        2    119.44  9773.8 552.08
## - conformity:nationality      2    141.00  9795.4 552.32
## - hedonism:nationality        2    227.20  9881.6 553.27
## - purposeinlife:nationality   2    229.16  9883.5 553.29
## - benevolence:nationality     2    299.42  9953.8 554.07
## - persgrowth:nationality      2    305.87  9960.3 554.14
## <none>                                     9654.4 554.74
## - selfdirection:nationality   2    374.25 10028.6 554.88
## - autonomy:nationality        2    376.35 10030.7 554.90
## - envmattery:nationality      2    401.94 10056.3 555.18
## - power:nationality           2    449.62 10104.0 555.70
## - universalism:nationality    2    680.26 10334.6 558.16
## 
## Step:  AIC=550.75
## ei ~ conformity:nationality + tradition:nationality + benevolence:nationality + 
##     universalism:nationality + selfdirection:nationality + stimulation:nationality + 
##     hedonism:nationality + power:nationality + security:nationality + 
##     positiverelat:nationality + autonomy:nationality + envmattery:nationality + 
##     persgrowth:nationality + purposeinlife:nationality + selfacceptance:nationality
## 
##                              Df Sum of Sq     RSS    AIC
## - nationality:selfacceptance  2     32.70  9688.3 547.12
## - nationality:positiverelat   2     38.91  9694.5 547.19
## - nationality:stimulation     2     49.23  9704.8 547.30
## - nationality:tradition       2    112.37  9767.9 548.01
## - nationality:security        2    120.34  9775.9 548.10
## - conformity:nationality      2    169.47  9825.0 548.65
## - nationality:hedonism        2    237.76  9893.3 549.40
## - nationality:purposeinlife   2    290.19  9945.7 549.98
## - nationality:persgrowth      2    307.42  9963.0 550.17
## <none>                                     9655.6 550.75
## - nationality:benevolence     2    371.73 10027.3 550.87
## - nationality:autonomy        2    375.38 10030.9 550.91
## - nationality:selfdirection   2    380.51 10036.1 550.96
## - nationality:envmattery      2    404.24 10059.8 551.22
## - nationality:power           2    505.88 10161.4 552.32
## - nationality:universalism    2    702.94 10358.5 554.41
## 
## Step:  AIC=547.12
## ei ~ conformity:nationality + nationality:tradition + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + nationality:stimulation + 
##     nationality:hedonism + nationality:power + nationality:security + 
##     nationality:positiverelat + nationality:autonomy + nationality:envmattery + 
##     nationality:persgrowth + nationality:purposeinlife
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:positiverelat  2     36.00  9724.3 543.52
## - nationality:stimulation    2     61.65  9749.9 543.81
## - nationality:security       2    108.54  9796.8 544.33
## - nationality:tradition      2    113.42  9801.7 544.39
## - conformity:nationality     2    174.87  9863.1 545.07
## - nationality:hedonism       2    234.12  9922.4 545.72
## - nationality:persgrowth     2    294.21  9982.5 546.38
## - nationality:purposeinlife  2    312.17 10000.4 546.57
## <none>                                    9688.3 547.12
## - nationality:autonomy       2    365.67 10053.9 547.16
## - nationality:envmattery     2    406.58 10094.8 547.60
## - nationality:selfdirection  2    419.37 10107.6 547.74
## - nationality:power          2    499.83 10188.1 548.60
## - nationality:benevolence    2    502.32 10190.6 548.63
## - nationality:universalism   2    733.96 10422.2 551.08
## 
## Step:  AIC=543.52
## ei ~ conformity:nationality + nationality:tradition + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + nationality:stimulation + 
##     nationality:hedonism + nationality:power + nationality:security + 
##     nationality:autonomy + nationality:envmattery + nationality:persgrowth + 
##     nationality:purposeinlife
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:stimulation    2     68.85  9793.1 540.29
## - nationality:security       2    119.46  9843.7 540.85
## - nationality:tradition      2    145.33  9869.6 541.14
## - conformity:nationality     2    154.22  9878.5 541.24
## - nationality:hedonism       2    206.29  9930.5 541.81
## - nationality:persgrowth     2    263.17  9987.4 542.43
## - nationality:purposeinlife  2    294.10 10018.4 542.77
## <none>                                    9724.3 543.52
## - nationality:selfdirection  2    388.90 10113.2 543.80
## - nationality:autonomy       2    393.37 10117.6 543.84
## - nationality:envmattery     2    420.45 10144.7 544.14
## - nationality:power          2    484.11 10208.4 544.82
## - nationality:benevolence    2    520.62 10244.9 545.21
## - nationality:universalism   2    765.15 10489.4 547.78
## 
## Step:  AIC=540.29
## ei ~ conformity:nationality + nationality:tradition + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + nationality:hedonism + 
##     nationality:power + nationality:security + nationality:autonomy + 
##     nationality:envmattery + nationality:persgrowth + nationality:purposeinlife
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:security       2    120.47  9913.6 537.62
## - nationality:tradition      2    139.79  9932.9 537.84
## - conformity:nationality     2    184.19  9977.3 538.32
## - nationality:hedonism       2    186.18  9979.3 538.34
## - nationality:persgrowth     2    245.69 10038.8 538.99
## - nationality:purposeinlife  2    342.04 10135.1 540.03
## - nationality:envmattery     2    359.81 10152.9 540.22
## <none>                                    9793.1 540.29
## - nationality:selfdirection  2    381.74 10174.8 540.46
## - nationality:autonomy       2    428.18 10221.3 540.96
## - nationality:benevolence    2    457.64 10250.7 541.27
## - nationality:power          2    706.20 10499.3 543.88
## - nationality:universalism   2    720.32 10513.4 544.03
## 
## Step:  AIC=537.62
## ei ~ conformity:nationality + nationality:tradition + nationality:benevolence + 
##     nationality:universalism + nationality:selfdirection + nationality:hedonism + 
##     nationality:power + nationality:autonomy + nationality:envmattery + 
##     nationality:persgrowth + nationality:purposeinlife
## 
##                             Df Sum of Sq     RSS    AIC
## - nationality:tradition      2    112.33 10025.9 534.85
## - nationality:persgrowth     2    238.24 10151.8 536.21
## - nationality:hedonism       2    242.51 10156.1 536.26
## <none>                                    9913.6 537.62
## - nationality:selfdirection  2    407.08 10320.7 538.01
## - nationality:purposeinlife  2    421.93 10335.5 538.17
## - conformity:nationality     2    430.78 10344.3 538.26
## - nationality:envmattery     2    489.31 10402.9 538.88
## - nationality:autonomy       2    510.84 10424.4 539.10
## - nationality:benevolence    2    576.23 10489.8 539.78
## - nationality:power          2    659.52 10573.1 540.64
## - nationality:universalism   2    944.75 10858.3 543.55
## 
## Step:  AIC=534.85
## ei ~ conformity:nationality + nationality:benevolence + nationality:universalism + 
##     nationality:selfdirection + nationality:hedonism + nationality:power + 
##     nationality:autonomy + nationality:envmattery + nationality:persgrowth + 
##     nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:hedonism       2    200.79 10227 533.01
## - nationality:persgrowth     2    210.49 10236 533.12
## <none>                                   10026 534.85
## - nationality:selfdirection  2    396.13 10422 535.08
## - nationality:purposeinlife  2    432.12 10458 535.45
## - nationality:envmattery     2    481.42 10507 535.96
## - nationality:autonomy       2    504.71 10531 536.21
## - nationality:power          2    621.43 10647 537.41
## - nationality:benevolence    2    642.12 10668 537.62
## - conformity:nationality     2    777.61 10804 538.99
## - nationality:universalism   2   1071.26 11097 541.92
## 
## Step:  AIC=533.01
## ei ~ conformity:nationality + nationality:benevolence + nationality:universalism + 
##     nationality:selfdirection + nationality:power + nationality:autonomy + 
##     nationality:envmattery + nationality:persgrowth + nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:persgrowth     2    159.99 10387 530.71
## - nationality:selfdirection  2    342.78 10570 532.61
## - nationality:purposeinlife  2    367.38 10594 532.86
## <none>                                   10227 533.01
## - nationality:envmattery     2    447.28 10674 533.68
## - nationality:benevolence    2    534.72 10761 534.57
## - nationality:power          2    644.64 10871 535.68
## - conformity:nationality     2    686.96 10914 536.10
## - nationality:autonomy       2    704.16 10931 536.27
## - nationality:universalism   2   1075.90 11303 539.92
## 
## Step:  AIC=530.71
## ei ~ conformity:nationality + nationality:benevolence + nationality:universalism + 
##     nationality:selfdirection + nationality:power + nationality:autonomy + 
##     nationality:envmattery + nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:selfdirection  2    343.01 10730 530.25
## - nationality:envmattery     2    379.81 10766 530.62
## <none>                                   10387 530.71
## - nationality:purposeinlife  2    487.53 10874 531.71
## - nationality:benevolence    2    492.46 10879 531.75
## - nationality:autonomy       2    636.26 11023 533.19
## - conformity:nationality     2    642.31 11029 533.25
## - nationality:power          2    670.33 11057 533.52
## - nationality:universalism   2   1258.63 11645 539.17
## 
## Step:  AIC=530.25
## ei ~ conformity:nationality + nationality:benevolence + nationality:universalism + 
##     nationality:power + nationality:autonomy + nationality:envmattery + 
##     nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:envmattery     2    232.15 10962 528.58
## - nationality:benevolence    2    289.64 11019 529.15
## <none>                                   10730 530.25
## - conformity:nationality     2    570.56 11300 531.89
## - nationality:autonomy       2    620.00 11350 532.37
## - nationality:power          2    833.88 11564 534.40
## - nationality:purposeinlife  2    844.33 11574 534.50
## - nationality:universalism   2   1110.92 11841 536.99
## 
## Step:  AIC=528.58
## ei ~ conformity:nationality + nationality:benevolence + nationality:universalism + 
##     nationality:power + nationality:autonomy + nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## - nationality:benevolence    2    270.07 11232 527.23
## <none>                                   10962 528.58
## - conformity:nationality     2    424.55 11386 528.72
## - nationality:power          2    754.79 11717 531.84
## - nationality:autonomy       2    775.08 11737 532.03
## - nationality:purposeinlife  2    900.56 11862 533.19
## - nationality:universalism   2   1476.84 12439 538.36
## 
## Step:  AIC=527.23
## ei ~ conformity:nationality + nationality:universalism + nationality:power + 
##     nationality:autonomy + nationality:purposeinlife
## 
##                             Df Sum of Sq   RSS    AIC
## <none>                                   11232 527.23
## - nationality:purposeinlife  2    751.47 11983 530.29
## - conformity:nationality     2    784.31 12016 530.59
## - nationality:power          2    872.94 12105 531.39
## - nationality:autonomy       2    895.85 12128 531.60
## - nationality:universalism   2   2080.83 13313 541.76
summary(final_model)
## 
## Call:
## lm(formula = ei ~ conformity:nationality + nationality:universalism + 
##     nationality:power + nationality:autonomy + nationality:purposeinlife, 
##     data = thes)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -31.177  -5.295  -0.100   7.784  21.449 
## 
## Coefficients:
##                              Estimate Std. Error t value   Pr(>|t|)    
## (Intercept)                   52.4515    10.9583   4.786 0.00000601 ***
## conformity:nationalityfr       4.3752     1.7881   2.447    0.01619 *  
## conformity:nationalityrus      1.7340     1.8800   0.922    0.35861    
## nationalityfr:universalism     4.9084     1.4618   3.358    0.00112 ** 
## nationalityrus:universalism    4.3993     1.6180   2.719    0.00775 ** 
## nationalityfr:power           -3.0948     1.1229  -2.756    0.00698 ** 
## nationalityrus:power          -0.1503     1.0051  -0.150    0.88145    
## nationalityfr:autonomy        -0.1927     0.2281  -0.845    0.40021    
## nationalityrus:autonomy       -0.5296     0.1895  -2.796    0.00623 ** 
## nationalityfr:purposeinlife   -0.2747     0.1336  -2.056    0.04247 *  
## nationalityrus:purposeinlife   0.2263     0.1782   1.269    0.20728    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 10.71 on 98 degrees of freedom
## Multiple R-squared:  0.348,  Adjusted R-squared:  0.2814 
## F-statistic:  5.23 on 10 and 98 DF,  p-value: 0.000003978

dispersion for environmental identity across birth’s places

# Combine the variable from different data frames into a list
data_list <- list(child_place1$ei, child_place2$ei, child_place3$ei, child_place4$ei)

# Create boxplots for the combined data
boxplot(data_list, names = c("child_place1", "child_place2", "child_place3", 'child_place4'), main = "Boxplot Comparison", xlab = "place of birth", ylab = "environmental identity")

#thes$ei <- factor(thes$ei)
#manova_result <- aov(child ~ ei, data = thes)
#summary(manova_result)
#TukeyHSD(manova_result)

#попарно сравниваю ид с природой (значимые различия - 1\4, 2\4)
t.test(child_place1$ei, child_place4$ei)
## 
##  Welch Two Sample t-test
## 
## data:  child_place1$ei and child_place4$ei
## t = -2.7538, df = 21.405, p-value = 0.01177
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -16.773273  -2.349176
## sample estimates:
## mean of x mean of y 
##  57.93878  67.50000
t.test(child_place4$ei[child_place4$nationality == "fr"], 
      child_place4$ei[child_place4$nationality == "rus"])   # между странами в детстве не отличается
## 
##  Welch Two Sample t-test
## 
## data:  child_place4$ei[child_place4$nationality == "fr"] and child_place4$ei[child_place4$nationality == "rus"]
## t = -0.58438, df = 10.383, p-value = 0.5714
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -15.100667   8.800667
## sample estimates:
## mean of x mean of y 
##     66.60     69.75
t.test(adult_place4$ei[adult_place4$nationality == "fr"], 
      adult_place4$ei[adult_place4$nationality == "rus"])  # во взрослости тоже не отличается
## 
##  Welch Two Sample t-test
## 
## data:  adult_place4$ei[adult_place4$nationality == "fr"] and adult_place4$ei[adult_place4$nationality == "rus"]
## t = -1.2139, df = 1.4706, p-value = 0.3844
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -57.92966  38.92966
## sample estimates:
## mean of x mean of y 
##      56.0      65.5

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

save(thes, file = "thes.RData")