I am confused with the following:
I found the next examples in a website https://stats.oarc.ucla.edu/r/dae/zero-truncated-negative-binomial/
library(foreign)
dat <- read.dta("https://stats.idre.ucla.edu/stat/data/ztp.dta")
dat <- within(dat, {
hmo <- factor(hmo)
died <- factor(died)
})
summary(dat)
library(VGAM)
m1 <- vglm(stay ~ age + hmo + died, family = posnegbinomial(),
data = dat)
summary(m1)
m2 <- glm(stay ~ age + hmo + died, family = gaussian(),
data = dat)
summary(m2)
Why in model 1 m1 the summary comes with a z value and in the glm with familiy gaussian the summary comes with t value??
I have seen in other examples that the first one always gives the z value in the third column and in the model 2 m2 the summary comes with t value.
m2, but with apoissonfamily. And check?summary.glmand?VGAM::summaryvglm– Alex J Jun 06 '23 at 22:36glmwith Gaussian, the dispersion (residual variance) is estimated. So according to its documentation, it'st. Forglmwith Poisson, the dispersion is assumed to be 1, so according to its documentation, it'st. For anything from VGAM, it "is labelled z value regardless of whether the dispersion is estimated or known" – Alex J Jun 07 '23 at 23:16