library(pglm)
data(Unions)
head(Unions,10)
anb <- pglm(union ~ wage + exper + rural, Unions, family = binomial('probit'),
            model = "pooling",  method = "bfgs", print.level = 3, R = 5)
summary(anb)
predict(anb,data=Unions)
#######然后报错
Error in UseMethod("predict") :
  "predict"没有适用于"c('maxLik', 'maxim')"目标对象的方法

请问这是什么情况？另外，如果不能用predict的话，那该怎么做预测，谢谢。