[R-br] função BiCopKPlot

[Marcus Nunes]

Natália,

tenta rodar

library(CDVine)
CDVine:::H
CDVine:::W

Att,


On Mon, Mar 24, 2014 at 4:30 PM, Natalia Martins <nsmbarreto em 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 em 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/
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140324/78866732/attachment.html>