One-truncated Zeta Distribution Family Function
Estimates the parameter of the 1-truncated zeta distribution.
otzeta(lshape = "loglink", ishape = NULL, gshape = exp((-4:3)/4), zero = NULL)
lshape, ishape, gshape, zero |
Same as |
The 1-truncated zeta distribution is the ordinary zeta
distribution but with the probability of one being 0. Thus the
other probabilities are scaled up (i.e., divided by 1-P[Y=1]).
The mean is returned by default as the fitted values.
More details can be found at zetaff.
Multiple responses are handled.
An object of class "vglmff" (see vglmff-class).
The object is used by modelling functions such as vglm,
and vgam.
T. W. Yee
## Not run: odata <- data.frame(x2 = runif(nn <- 1000)) # Artificial data odata <- transform(odata, shape = loglink(-0.25 + x2, inverse = TRUE)) odata <- transform(odata, y1 = rotzeta(nn, shape)) with(odata, table(y1)) ofit <- vglm(y1 ~ x2, otzeta, data = odata, trace = TRUE, crit = "coef") coef(ofit, matrix = TRUE) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.