## Confronti multipli ## 1) la Anova necessita di 2 ipotesi: ## la normalita' ## la omoschedasticita' ## (per fortuna che nel 2005 hanno "multcomp" "sandwich") ## 2) la vecchia idea di Bonferroni ## la buona idea di Tukey ## l'attualita' attach(iris) relazione = Petal.Length ~ Species modello50 = lm(relazione, data = iris) summary(modello50) boxplot(relazione) attach(warpbreaks) head(warpbreaks) relazione = breaks ~ tension modello51 = lm(relazione, data = warpbreaks) summary(modello51) boxplot(relazione) #### tooth = read.csv( file.choose() , header = TRUE) indirizzo = "http://www.dmi.units.it/borelli/dataset/tooth.csv" tooth = read.csv( indirizzo , header = TRUE) attach(tooth) head(tooth) relazione = areainfl ~ il1b modello52 = lm(relazione, data = tooth) summary(modello52) boxplot(relazione) ## Tukey propone il "Test Onesto delle Differenze Significative" ## TukeyHSD modello52bis = aov(relazione, data = tooth) TukeyHSD(modello52bis) plot( TukeyHSD(modello52bis) ) # la soluzione ottimale di Bretz, Hothorn, Westfall library(multcomp) library(sandwich) posthoc = glht(modello52, linfct = mcp(il1b = "Tukey"), vcov = sandwich) summary(posthoc)