Hi, I am trying to use mice() with an aov() model.I can fun mice() fine but when I try to pool the results it doesn't work. For example, I can run the following: fit=with(data=imp,exp=aov(Y~B*P*T + Error(S/(B*P*T)))) pool(fit) # doesn't work Where B, P, and T are within subject (S) factors. However, I can't go the next step to pool these results as I am guessing this hasn't been implemented in mice yet. Any ideas if it is possible to pool these results and how I might go about doing this? Or how might I present some sort of pooled results instead of 10 separate imputations? Thanks, Chris [[alternative HTML version deleted]]