"qua.regressCOP" <-
function(f=0.5, u=seq(0.01,0.99, by=0.01), cop=NULL, para=NULL, ...) {
if(is.null(cop)) {
warning("must have copula argument specified, returning NULL")
return(NULL)
}
V <- sapply(1:length(u), function(i) {
v <- derCOPinv(cop=cop, u[i], f, para=para, ...)
if(is.na(v)) {
warning("could not uniroot in derCOPinv, skipping i=",i," with U=", u[i])
return(NA)
}
return(v) } )
z <- data.frame(U=u,V=V)
return(z)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.