options(scipen = 999)

setwd("C:/Users/irina/Documents/DND/EOHI/eohi1")

data <- read.csv("ehi1.csv")

# Check the levels of the demo_edu variable
print("Current levels of demo_edu variable:")
print(levels(factor(data$demo_edu)))

# Also show the unique values and their frequencies
print("\nUnique values and frequencies:")
print(table(data$demo_edu, useNA = "ifany"))