Density and random generation for the Borel-Tanner distribution.

dbort(x, Qsize = 1, a = 0.5, log = FALSE)
rbort(n, Qsize = 1, a = 0.5)

Arguments

x

vector of quantiles.

n

number of observations. Must be a positive integer of length 1.

Qsize, a

See borel.tanner.

log

Logical. If log = TRUE then the logarithm of the density is returned.

Value

dbort gives the density, rbort generates random deviates.

Author

T. W. Yee

Details

See borel.tanner, the VGAM family function for estimating the parameter, for the formula of the probability density function and other details.

Warning

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.

See also

Examples

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