After we are satisfied with our model, we can jump into the most interesting part: obtaining confidence intervals for the means of the fixed effects that we estimated. This can be achieved using the following function (note that we are indicating the lsmeans function, which we want to use to get the analysis for the Strategy effect):

result_lsmeans = lsmeans(E,pairwise ~ Strategy)
print(result_lsmeans)

The preceding code generate the following output of lsmeans results:

This prints two parts. The first part part has the least-squares means for the mean effects for each level (along with the 95% confidence intervals), while the second part has the estimated differences for all possible contrasts. Contrasts are just differences between all levels. We can see that all of the contrasts are statistically significative.