Natália, tenta rodar
library(CDVine) CDVine:::H CDVine:::W Att, On Mon, Mar 24, 2014 at 4:30 PM, Natalia Martins <nsmbarr...@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@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@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.