<div dir="ltr">Consulte o arquivo fonte:<div><br></div><div>CDVine/R/BiCopChiPlot.R</div><div><br></div><div>Na linha 154 é definida a função H() e na linha 170 a função W()!</div><div><br></div><div>att,</div><div>FH</div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On 24 March 2014 16:36, Marcus Nunes <span dir="ltr"><<a href="mailto:marcus.nunes@gmail.com" target="_blank">marcus.nunes@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">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"><div><div class="h5">
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></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div class="h5">
<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><font color="#888888"><div style="font-family:georgia,serif;font-size:large"><br></div>-- <br><div dir="ltr"><br>Natália Martins<br>
Contato: <a href="tel:%2835%29%2091812482" value="+553591812482" target="_blank">(35) 91812482</a> <br></div>
</font></span></div>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.<span class="HOEnZb"><font color="#888888"><br></font></span></blockquote>
</div><span class="HOEnZb"><font color="#888888"><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>
</font></span></div></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></div>