Natália,

tenta rodar

library(CDVine)

CDVine:::H

CDVine:::W

Att,

On Mon, Mar 24, 2014 at 4:30 PM, Natalia Martins <nsmbarreto@gmail.com> wrote:

Prezados  boa tarde,

estou estudando a função BiCopKPlot utilizada para a modelagem de cópulas. 

Nesta função tem-se:

function (u1, u2, PLOT = TRUE, ...) 

{

    if (is.null(u1) == TRUE || is.null(u2) == TRUE) 

        stop("u1 and/or u2 are not set or have length zero.")

    if (length(u1) != length(u2)) 

        stop("Lengths of 'u1' and 'u2' do not match.")

    if (length(u1) < 2) 

        stop("Number of observations has to be at least 2.")

    if (PLOT != TRUE && PLOT != FALSE) 

        stop("The parameter 'PLOT' has to be set to 'TRUE' or 'FALSE'.")

    Wi <- W(u1, u2)

    Hi <- H(u1, u2)

    Hi.sort <- sort(Hi)

    n <- length(u1)

    W.in <- rep(NA, n)

    for (i in 1:n) {

        f = function(w) {

            w * (-log(w)) * (w - w * log(w))^(i - 1) * (1 - w + 

                w * log(w))^(n - i)

        }

        W.in[i] <- n * choose(n - 1, i - 1) * (integrate(f, lower = 0, 

            upper = 1)$value)

    }

    g <- function(w) {

        w - w * log(w)

    }

    if (PLOT) {

        plot(g, xlim = c(0, 1), ylim = c(0, 1), pch = "x", xlab = expression(W[1:n]), 

            ylab = "H", ...)

        points(W.in, Hi.sort, pch = "x", cex = 0.4, ...)

        abline(a = 0, b = 1)

    }

    else {

        kendall.plot.output <- list(W.in, Hi.sort)

        names(kendall.plot.output) <- c("W.in", "Hi.sort")

        return(kendall.plot.output)

    }

}

Porem, não encontrei no pdf da library(CDVine) como é definida a função H e W (em destaque).

Por favor se alguem souber e puder me ajudar.

Muito obrigada!

--

Natália Martins
Contato: (35) 91812482 

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

--

Marcus Nunes
http://marcusnunes.me/

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.