r - R の dpareto の関数

Question

Rで書かれたdparetoの機能を持っている人がいるのだろうか？

ライブラリをインストールできないため、R の組み込み関数を使用できません。

ありがとう！

Answer 1

たとえば、次のように、R パッケージから関数を取り出すのは比較的簡単です。

http://cran.r-project.org/web/packages/VGAM/VGAM.pdf

dpareto <- function(x, location, shape, log = FALSE) {
    if (!is.logical(log.arg <- log) || length(log) != 1)
        stop("bad input for argument 'log'")
    rm(log)

    L = max(length(x), length(location), length(shape)) 
    x = rep(x, length.out = L);
    location = rep(location, length.out = L);
    shape = rep(shape, length.out = L)

    logdensity = rep(log(0), length.out = L)
    xok = (x > location)
    logdensity[xok] = log(shape[xok]) + shape[xok] * log(location[xok]) -
            (shape[xok]+1) * log(x[xok])
    if (log.arg) logdensity else exp(logdensity)
}

alpha <- 3; k <- exp(1); x <- seq(2.8, 8, len = 300)
## Not run: 
plot(x, dpareto(x, location = alpha, shape = k), type = "l",
        main = "Pareto density split into 10 equal areas")