## TASK 1
termometers<-read.delim2("clipboard")
summary(termometers)
var.test(temp~manufacturer, data=termometers)
boxplot(temp~manufacturer, data=termometers)

## task 2
fert<-read.delim2("clipboard")
summary(fert)

aov.f<-aov(Biomass~Fertilizer, data=fert)
summary(aov.f)
# Df Sum Sq Mean Sq F value  Pr(>F)    
# Fertilizer   4  26740    6685   45.96 3.6e-11 ***
#   Residuals   25   3636     145                    
# ---
#   Signif. codes:  
#   0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

plot(aov.f)

boxplot(Biomass~Fertilizer, data=fert,
        las=2, ylim=c(70, 200))
text(1:5, tapply(fert$Biomass, fert$Fertilizer, 
     max)+8, c("a", "c", "c", "b", "c"))

TukeyHSD(aov.f)
