Displaying 1 result from an estimated 1 matches for "paretofit".
2003 Aug 28
2
ks.test()
...mbda, b = b)
class(result) <- "expt"
result}
dexpt<-function(x,lambda,b) dexp(x-b,rate=lambda)
pexpt<-function(x,lambda,b) pexp(x-b,rate=lambda)
qexpt<-function(p,lambda,b) qexp(p,rate=lambda)+b
rexpt<-function(n,lambda,b) rexp(n,rate=lambda)+b
#PARETO
paretofit<-function(x,b, plot.it = F, lty = 1){
x <- sort(x)
alpha <- length(x)/sum(log(x/b))
if(plot.it) { lines(x, dpareto(x, alpha = alpha, b = b), lty =
lty,col="red")}
result <- list(alpha = alpha, b = b)
class(result) <- "...