if (FALSE) {
data("bioChemists")
zip <- zeroinfl(art ~ . | ., data = bioChemists, EM = TRUE)
phat <- predprob(zip)
newdata <- expand.grid(list(fem="Men",mar="Married",
kid5=1,phd=3.103,
ment=0:77))
phat <- predprob(zip, newdata = newdata)
}
Run the code above in your browser using DataLab