<div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif;font-size:large">Prezados  boa tarde,</div><div class="gmail_default"><span style="font-family:georgia,serif;font-size:large">estou estudando a função </span><font face="georgia, serif" size="4">BiCopKPlot utilizada para a modelagem de cópulas. </font></div>
<div class="gmail_default"><font face="georgia, serif" size="4">Nesta função tem-se:</font></div><div class="gmail_default"><font face="georgia, serif" size="4"><br></font></div><div class="gmail_default"><font face="georgia, serif" size="4"><div class="gmail_default">
function (u1, u2, PLOT = TRUE, ...) </div><div class="gmail_default">{</div><div class="gmail_default">    if (is.null(u1) == TRUE || is.null(u2) == TRUE) </div><div class="gmail_default">        stop("u1 and/or u2 are not set or have length zero.")</div>
<div class="gmail_default">    if (length(u1) != length(u2)) </div><div class="gmail_default">        stop("Lengths of 'u1' and 'u2' do not match.")</div><div class="gmail_default">    if (length(u1) < 2) </div>
<div class="gmail_default">        stop("Number of observations has to be at least 2.")</div><div class="gmail_default">    if (PLOT != TRUE && PLOT != FALSE) </div><div class="gmail_default">        stop("The parameter 'PLOT' has to be set to 'TRUE' or 'FALSE'.")</div>
<div class="gmail_default">    Wi <- <span style="background-color:rgb(255,255,0)">W</span>(u1, u2)</div><div class="gmail_default">    Hi <-<span style="background-color:rgb(255,255,0)"> H</span>(u1, u2)</div><div class="gmail_default">
    Hi.sort <- sort(Hi)</div><div class="gmail_default">    n <- length(u1)</div><div class="gmail_default">    W.in <- rep(NA, n)</div><div class="gmail_default">    for (i in 1:n) {</div><div class="gmail_default">
        f = function(w) {</div><div class="gmail_default">            w * (-log(w)) * (w - w * log(w))^(i - 1) * (1 - w + </div><div class="gmail_default">                w * log(w))^(n - i)</div><div class="gmail_default">
        }</div><div class="gmail_default">        W.in[i] <- n * choose(n - 1, i - 1) * (integrate(f, lower = 0, </div><div class="gmail_default">            upper = 1)$value)</div><div class="gmail_default">    }</div>
<div class="gmail_default">    g <- function(w) {</div><div class="gmail_default">        w - w * log(w)</div><div class="gmail_default">    }</div><div class="gmail_default">    if (PLOT) {</div><div class="gmail_default">
        plot(g, xlim = c(0, 1), ylim = c(0, 1), pch = "x", xlab = expression(W[1:n]), </div><div class="gmail_default">            ylab = "H", ...)</div><div class="gmail_default">        points(W.in, Hi.sort, pch = "x", cex = 0.4, ...)</div>
<div class="gmail_default">        abline(a = 0, b = 1)</div><div class="gmail_default">    }</div><div class="gmail_default">    else {</div><div class="gmail_default">        kendall.plot.output <- list(W.in, Hi.sort)</div>
<div class="gmail_default">        names(kendall.plot.output) <- c("W.in", "Hi.sort")</div><div class="gmail_default">        return(kendall.plot.output)</div><div class="gmail_default">    }</div><div class="gmail_default">
}</div><div class="gmail_default"><br></div><div class="gmail_default">Porem, não encontrei no pdf da library(CDVine) como é definida a função H e W (em destaque).</div><div class="gmail_default">Por favor se alguem souber e puder me ajudar.</div>
<div class="gmail_default"><br></div><div class="gmail_default">Muito obrigada!</div></font></div><div class="gmail_default" style="font-family:georgia,serif;font-size:large"><br></div>-- <br><div dir="ltr"><br>Natália Martins<br>
Contato: (35) 91812482 <br></div>
</div>