dbort(x, Qsize=1, a=0.5)rbort(n, Qsize=1, a=0.5)
borel.tanner
.dbort
gives the density, rbort
generates random deviates.rbort
, therefore
values of a
close to 1 will result in long (or infinite!)
computational times.
The default value of a
is subjective.borel.tanner
, the borel.tanner
.qsize = 1; a = 0.5
x = qsize:(qsize+10)
plot(x, dbort(x, qsize, a), type="h", las=1, col="blue",
ylab=paste("fbort(qsize=", qsize, ", a=", a, ")"),
main="Borel-Tanner density function")
Run the code above in your browser using DataLab