"diagCOP" <-
function(cop=NULL, para=NULL, secondary=FALSE,
ploton=TRUE, lines=TRUE, delt=0.005, ...) {
if(ploton) {
plot(c(0,1), c(0,1), type="n",
xlab="U=V, NONEXCEEDANCE PROBABILITY",
ylab="H, NONEXCEEDANCE PROBABILITY")
}
TT <- seq(0+delt, 1-delt, delt)
C <- sapply(TT, function(x) {
y <- x
if(secondary) y <- 1 - y
return( cop(x,y, para=para, ...)) } )
if(lines & ! is.null(dev.list())) lines(TT, C, ...)
txt <- "primary"
if(secondary) txt <- "secondary"
return(list(t=TT, diagcop=C, diagtype=txt))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.