glht for lme object with significant interaction term
I'm working on some data from an experiment on the breeding behavior of birds. In short, I have been measuring how the time spent on performing a certain task (variable 'mean_on_active') differs over time (variable 'day', 2 levels) across three experimental categories (variable 'treat'; levels 'C', 'R', 'E'). The model shows a significant interaction between treatment and day. To have a closer look at this, I would like to do multiple comparisons for treatment levels within day (i.e. across treatments for each day in turn). I have browsed the forum for a way of doing this using the glht function, but haven't found a good solution to my problem.
Any hints on how to proceed, or on other methods that might be (more) appropriate? Final model output below, data attached.sub.data.txt
Linear mixed-effects model fit by REML
AIC BIC logLik
1051.779 1075.757 -517.8896
Formula: ~1 | id
StdDev: 6.276309 5.473167