bortUC.RdDensity and random generation for the Borel-Tanner distribution.
dbort(x, Qsize = 1, a = 0.5, log = FALSE)
rbort(n, Qsize = 1, a = 0.5)number of observations. Must be a positive integer of length 1.
See borel.tanner.
Logical.
If log = TRUE then the logarithm of the density is
returned.
dbort gives the density,
rbort generates random deviates.
See borel.tanner, the VGAM family function
for estimating the parameter,
for the formula of the probability density function and other
details.
Looping is used for rbort, therefore
values of a close to 1 will result in long (or infinite!)
computational times.
The default value of a is subjective.
if (FALSE) qsize <- 1; a <- 0.5; x <- qsize:(qsize+10)
#> Error: object 'qsize' not found
plot(x, dbort(x, qsize, a), type = "h", las = 1, col = "blue",
ylab = paste("fbort(qsize=", qsize, ", a=", a, ")"),
log = "y", main = "Borel-Tanner density function") # \dontrun{}
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'plot': object 'x' not found