I have just run a repeated measures factorial ANOVA using the aov() function in r. The one problem I am having is that the output does not include partial eta squared, so I was wondering what would be the best way to get that information.