R - мыши - tukeyHSD

Question

У меня есть база данных с качественными и количественными переменными. Я вменяю недостающие данные с помощью мышей.

Я хочу сделать анализ с TukeyHSD, но не могу найти способ сделать это с мышами. TukeyHSD использует только объекты aov или lm, поэтому при попытке использовать его с mi.anova выдается следующая ошибка.

    bdd=data.frame(
    gender=c("M","F","M", "M", "M", "F"),
    choice=c(1,2,NA,1,1,1),
    gardes=c(0,0,0,5,7,NA),
    EE1=c(3,4,1,NA,3,0),
    EE2=c(2,5,1,3,3,0),
    EE3=c(3,NA,1,5,3,0),
    EE4=c(3,6,1,2,3,0),
    EE5=c(1,4,1,2,3,5),
    EE6=c(3,1,1,3,3,4),
    EE7=c(5,0,1,5,3,5),
    EE8=c(2,6,1,1,3,3),
    EE9=c(3,4,1,6,3,4)
    )

    imput = mice(bdd, m=2, printFlag=FALSE,  seed=2525)
    TukeyHSD(mi.anova(mi.res=imput,formula="EE1~choice"))


    Univariate ANOVA for Multiply Imputed Data (Type 2)  

    lm Formula:  AC~psychotherapie
    R^2=0 
    ..........................................................................
    ANOVA Table 
                          SSQ df1      df2 F value Pr(>F)  eta2 partial.eta2
    psychotherapie 7.6538e-01   1 2952.643 -0.0374      1 4e-05        4e-05
    Residual       1.9894e+04  NA       NA      NA     NA    NA           NA
    Error in UseMethod("TukeyHSD") : 
      no applicable method for 'TukeyHSD' applied to an object of class "list"

Кто-нибудь нашел обходной путь? TIA!!