Функция для dpareto в R

Я задаюсь вопросом, есть ли у кого-либо функция для dpareto, написанного в R?

Я не в состоянии использовать в построенной функции в R, поскольку я не могу установить библиотеки.

Спасибо!

1
"анализ первопричины" применяется здесь: удар и крик до безотносительно гения блокируют ваш веб-доступ, и/или местные права решает проблему. Можно знать это, но если вы загружаете двоичный файл "почтового индекса", можно установить в местном масштабе, не нуждаясь в активной ссылке с хранилищем.
добавлено автор Carl Witthoft, источник

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")
2
добавлено
Правда, но it' s еще быстрее и легче захватить $FUNCTION.R из zip-файла (или источник tarball, если код был предварительно собран), чем попытаться разорвать его из PDF.
добавлено автор Carl Witthoft, источник