Fits a 1-inflated Zipf distribution.
1 2 3
For initial values.
The 1-inflated Zipf distribution is a mixture
distribution of the
distribution with some probability of obtaining a (structural) 1.
Thus there are two sources for obtaining the value 1.
This distribution is written here
in a way that retains a similar notation to the
zero-inflated Poisson, i.e., the
probability P[Y=1] involves another parameter phi.
This family function can handle multiple responses.
An object of class
The object is used by modelling functions such as
Under- or over-flow may occur if the data is ill-conditioned. Lots of data is needed to estimate the parameters accurately. Usually, probably the shape parameter is best modelled as intercept-only.
Thomas W. Yee
1 2 3 4 5 6 7 8 9 10
## Not run: odata <- data.frame(x2 = runif(nn <- 1000)) # Artificial data odata <- transform(odata, pstr1 = logitlink(-1 + x2, inverse = TRUE), myN = 10, shape = exp(-0.5)) odata <- transform(odata, y1 = roizipf(nn, N = myN, s = shape, pstr1 = pstr1)) with(odata, table(y1)) fit1 <- vglm(y1 ~ x2, oizipf(zero = "shape"), data = odata, trace = TRUE) coef(fit1, matrix = TRUE) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.