dane = read.csv('~/Pobrane/dane1.csv', sep = ';')
library(nlme)
library(multcomp)

model = lme(imp.1.20 ~ szczep, random = ~1|mysz, data = dane)
print(model)
anova(model)

summary(glht(model, linfct=mcp(szczep="Tukey")))