# Student's t-test
t.test(netustm, mu = 30.6)
t.test(netustm, mu = 30.6, alternative = "greater")
t.test(netustm, mu = 32, alternative = "less")

# test to compare variances
var.test(netustm[gndr == "Male"], netustm[gndr == "Female"])

# two sample Student's/Welsh's t-test to compare means
t.test(netustm[gndr == "Male"], netustm[gndr == "Female"], var.equal = TRUE)
t.test(netustm[gndr == "Male"], netustm[gndr == "Female"])

# paired.test (not run!)
# t.test(x, y, paired = TRUE)

### one-way ANOVA
model1 <- lm(netustm ~ nuts3)
anova(model1)
TukeyHSD(aov(model1), which = "nuts3")

bartlett.test(netustm ~ nuts3)
library(car)
leveneTest(netustm ~ nuts3, center=mean)