"bicCOP" <-
function(u, v=NULL, cop=NULL, para=NULL, m=NA, ...) {
if(is.null(v)) {
v <- u[,2]
u <- u[,1]
}
if(is.null(cop)) {
warning("must have copula argument specified, returning NULL")
return(NULL)
}
theo.cop <- COP(u, v, cop=cop, para=para, ...)
emp.cop <- EMPIRcop(u,v, para=data.frame(u,v), ...)
n <- length(theo.cop); k <- length(emp.cop)
if(n != k) return(Inf)
mse <- mean((theo.cop - emp.cop)^2)
if(is.na(m)) m <- length(para)
bic <- m*log(n) + n*log(mse)
return(bic)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.