# 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)