<div dir="ltr">Natália,<div><br></div><div>tenta rodar</div><div><br></div><div>library(CDVine)</div><div>CDVine:::H<br></div><div>CDVine:::W<br></div><div><br></div><div>Att,</div><div class="gmail_extra"><br><br><div class="gmail_quote">

On Mon, Mar 24, 2014 at 4:30 PM, Natalia Martins <span dir="ltr"><<a href="mailto:nsmbarreto@gmail.com" target="_blank">nsmbarreto@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div dir="ltr"><div style="font-family:georgia,serif;font-size:large">Prezados  boa tarde,</div><div><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><font face="georgia, serif" size="4">Nesta função tem-se:</font></div><div><font face="georgia, serif" size="4"><br></font></div><div><font face="georgia, serif" size="4"><div>
function (u1, u2, PLOT = TRUE, ...) </div><div>{</div><div>    if (is.null(u1) == TRUE || is.null(u2) == TRUE) </div><div>        stop("u1 and/or u2 are not set or have length zero.")</div>
<div>    if (length(u1) != length(u2)) </div><div>        stop("Lengths of 'u1' and 'u2' do not match.")</div><div>    if (length(u1) < 2) </div>
<div>        stop("Number of observations has to be at least 2.")</div><div>    if (PLOT != TRUE && PLOT != FALSE) </div><div>        stop("The parameter 'PLOT' has to be set to 'TRUE' or 'FALSE'.")</div>


<div>    Wi <- <span style="background-color:rgb(255,255,0)">W</span>(u1, u2)</div><div>    Hi <-<span style="background-color:rgb(255,255,0)"> H</span>(u1, u2)</div><div>
    Hi.sort <- sort(Hi)</div><div>    n <- length(u1)</div><div>    W.in <- rep(NA, n)</div><div>    for (i in 1:n) {</div><div>
        f = function(w) {</div><div>            w * (-log(w)) * (w - w * log(w))^(i - 1) * (1 - w + </div><div>                w * log(w))^(n - i)</div><div>
        }</div><div>        W.in[i] <- n * choose(n - 1, i - 1) * (integrate(f, lower = 0, </div><div>            upper = 1)$value)</div><div>    }</div>
<div>    g <- function(w) {</div><div>        w - w * log(w)</div><div>    }</div><div>    if (PLOT) {</div><div>
        plot(g, xlim = c(0, 1), ylim = c(0, 1), pch = "x", xlab = expression(W[1:n]), </div><div>            ylab = "H", ...)</div><div>        points(W.in, Hi.sort, pch = "x", cex = 0.4, ...)</div>


<div>        abline(a = 0, b = 1)</div><div>    }</div><div>    else {</div><div>        kendall.plot.output <- list(W.in, Hi.sort)</div>
<div>        names(kendall.plot.output) <- c("W.in", "Hi.sort")</div><div>        return(kendall.plot.output)</div><div>    }</div><div>
}</div><div><br></div><div>Porem, não encontrei no pdf da library(CDVine) como é definida a função H e W (em destaque).</div><div>Por favor se alguem souber e puder me ajudar.</div>
<div><br></div><div>Muito obrigada!</div></font></div><span class="HOEnZb"><font color="#888888"><div style="font-family:georgia,serif;font-size:large"><br></div>-- <br><div dir="ltr"><br>Natália Martins<br>
Contato: (35) 91812482 <br></div>
</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>

<div dir="ltr">Marcus Nunes<br><a href="http://marcusnunes.me/" target="_blank">http://marcusnunes.me/</a><br></div>
</div></div>