title: “Statistics project” author: “Sofia Miraglia” date: “2025-12-28” output: html_document summary(data[, c(“population”, “net_migration”, “gdp_pc”, “homicide_rate”, “total_employment”)]) hist(data\(population, main = "Histogram of Population", xlab = "Population", col = "lightblue", border = "white") hist(data\)net_migration, main = “Histogram of Net Migration”, xlab = “Net Migration”, col = “lightgreen”, border = “white”) hist(data\(gdp_pc, main = "Histogram of GDP per capita", xlab = "GDP per capita", col = "lightgray", border = "white") hist(data\)homicide_rate, main = “Histogram of Homicide Rate”, xlab = “Homicide Rate”, col = “pink”, border = “white”) hist(data\(total_employment, main = "Histogram of Total Employment", xlab = "Total Employment", col = "orange", border = "white") shapiro.test(data\)population) shapiro.test(data\(net_migration) shapiro.test(data\)gdp_pc) shapiro.test(data\(homicide_rate) shapiro.test(data\)total_employment) plot(data\(net_migration, data\)population, xlab = “Net Migration”, ylab = “Population”, main = “Population vs Net Migration”) plot(data\(gdp_pc, data\)population, xlab = “GDP per capita”, ylab = “Population”, main = “Population vs GDP per capita”) plot(data\(homicide_rate, data\)population, xlab = “Homicide Rate”, ylab = “Population”, main = “Population vs Homicide Rate”) plot(data\(total_employment, data\)population, xlab = “Total Employment”, ylab = “Population”, main = “Population vs Total Employment”) cor.test(data\(population, data\)net_migration, method = “pearson”) cor.test(data\(population, data\)gdp_pc, method = “pearson”) cor.test(data\(population, data\)homicide_rate, method = “pearson”) cor.test(data\(population, data\)total_employment, method = “pearson”) cor.test(data\(population, data\)net_migration, method = “spearman”) cor.test(data\(population, data\)gdp_pc, method = “spearman”) cor.test(data\(population, data\)total_employment, method = “spearman”) model <- lm(population ~ total_employment, data = data) summary(model) plot(data\(total_employment, data\)population, xlab = “Total Employment”, ylab = “Population”, main = “Regression: Population vs Total Employment”) abline(model, col = “red”, lwd = 2) aggregate(cbind(population, net_migration, gdp_pc, homicide_rate, total_employment) ~ region, data = data, FUN = mean)