pareto_density: Pareto density

In ACSWR: A Companion Package for the Book "A Course in Statistics with R"

Description

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.

Usage

pareto_density(x, scale, shape)

Arguments

x	the x value
scale	the scale parameter of Pareto RV
shape	the shape parameter of Pareto RV

Author(s)

Prabhanjan N. Tattar

See Also

VGAM::dpareto

Examples

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

Example output

[1] 9273.165

ACSWR documentation built on May 2, 2019, 6:53 a.m.