x <- rnsbeta(1e5, 5, 13, -4, 8)
xx <- seq(-20, 20, by = 0.1)
hist(x, 100, freq = FALSE)
lines(xx, dnsbeta(xx, 5, 13, -4, 8), col = "red")
hist(pnsbeta(x, 5, 13, -4, 8))
plot(ecdf(x))
lines(xx, pnsbeta(xx, 5, 13, -4, 8), col = "red", lwd = 2)
Run the code above in your browser using DataLab