#boxes_exp$sample.type:

#boxes_exp$outcome
[1] 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 1 0 1 0 0 0 0 0 0 1
1 1 1 1 0 1 1 1 1 0 1 1 1 1 0 1 0 0 0 0 0..... (levels: 0, 1)

#plot(effect("sample.type", m1), x.var="sample.type",
 main="outcome in relation to sample type",
 multiline=TRUE,confint=TRUE,ci.style="none")

#Using package visreg:
visreg(m1, "sample.type", type = "conditional",main="Outcome in relation
to Sample type", xlab = "Sample type", ylab = "outcome")

#mosaicplot(table(boxes_exp$sample.type, boxes_exp$outcome))