I am analysing for differences in seagrass cover (cover) over time (period), specifically if the seagrass cover differed significantly from the baseline (periodbaseline1 and periodbaseline2).
The data was found to be non-normal (shapiro-wilkins) and variance was not homogenous (Levene's Test). The data was log-transformed (cover_log) and although the data remained non-normal, the variance was homogenous. A glm was carried out and the output showed that qr6 was significantly different from periodbaseline1 (P=0.04322). I carried out a pair-wise comparison with the the multcomp package as I also needed to compare qr6 against periodbaseline2. Unexpectedly, the result for qr6-periodbaseline1 changed to being not significantly different (p=0.4665). In a case such as this, which result should I be taking, the one from the glm output or the multcomp output?