Estimating the (single) parameter of the 1-truncated logarithmic distribution.
otlog(lshape = "logitlink", gshape = ppoints(8), zero = NULL)
The 1-truncated logarithmic distribution is a logarithmic
distribution but with the probability of a one being zero. The
other probabilities are scaled to add to unity.
Some more details can be found at
Multiple responses are permitted.
An object of class
The object is used by modelling functions
T. W. Yee
## Not run: odata <- data.frame(y1 = rotlog(1000, shape = logitlink(1/3, inv = TRUE))) ofit <- vglm(y1 ~ 1, otlog, data = odata, trace = TRUE, crit = "c") coef(ofit, matrix = TRUE) Coef(ofit) with(odata, hist(y1, shape = TRUE, breaks = seq(0.5, max(y1) + 0.5, by = 1), border = "blue")) x <- seq(1, with(odata, max(y1)), by = 1) with(odata, lines(x, dotlog(x, Coef(ofit)), col = "orange", type = "h", lwd = 2)) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.