shape <- 0.75; pobs1 <- 0.10; x <- (-1):7
doalog(x, shape = shape, pobs1 = pobs1)
table(roalog(100, shape = shape, pobs1 = pobs1))
if (FALSE) x <- 0:10
barplot(rbind(doalog(x, shape = shape, pobs1 = pobs1),
dlog(x, shape = shape)),
beside = TRUE, col = c("blue", "orange"), cex.main = 0.7, las = 1,
ylab = "Probability", names.arg = as.character(x),
main = paste("OAL(shape = ", shape, ", pobs1 = ", pobs1,
") [blue] vs", " Logarithmic(shape = ", shape,
") [orange] densities", sep = ""))
Run the code above in your browser using DataLab