[R-br] função BiCopKPlot

FHRB Toledo fernandohtoledo em gmail.com
Segunda Março 24 16:38:42 BRT 2014


Consulte o arquivo fonte:

CDVine/R/BiCopChiPlot.R

Na linha 154 é definida a função H() e na linha 170 a função W()!

att,
FH


On 24 March 2014 16:36, Marcus Nunes <marcus.nunes em gmail.com> wrote:

> 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/
>
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140324/1eff80e7/attachment.html>


Mais detalhes sobre a lista de discussão R-br