Pareto density
A simple function is given here which returns the density function values for a Pareto RV. A more efficient implementation is obtainable in the function "dpareto" from the "VGAM" package.
pareto_density(x, scale, shape)
x |
the x value |
scale |
the scale parameter of Pareto RV |
shape |
the shape parameter of Pareto RV |
Prabhanjan N. Tattar
VGAM::dpareto
m <- 9184 n <- 103 b <- 10000 K <- 10 theta <- seq(1000,20000,500) plot(theta,as.numeric(sapply(theta,pareto_density,scale=b,shape=K)),"l", xlab=expression(theta),ylab="The Posterior Density") (n+1)*m/n
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.