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