```
# NOT RUN {
set.seed(1)
pneumo <- transform(pneumo, let = log(exposure.time),
x3 = rnorm(nrow(pneumo)))
(pfit <- vglm(cbind(normal, mild, severe) ~ let + x3,
propodds, data = pneumo))
score.stat(pfit) # No HDE here; should be similar to the next line:
coef(summary(pfit))[, "z value"] # Wald statistics computed at the MLE
summary(pfit, score0 = TRUE)
# }
```

Run the code above in your browser using DataLab