# NOT RUN {
x <- ExtremalDep:::simplex(2)
data <- evd::rbvevd(50, dep = 0.4, model = "log", mar1 = c(1,1,1))
Amd <- beed(data, x, 2, "md", "emp", 20, plot=TRUE)
Acfg <- beed(data, x, 2, "cfg", "emp", 20)
Aht <- beed(data, x, 2, "ht", "emp", 20)
lines(x[,1], Aht$A, lty = 1, col = 3)
lines(x[,1], Acfg$A, lty = 1, col = 2)
##################################
# Trivariate case
##################################
# }
# NOT RUN {
x <- ExtremalDep:::simplex(3)
data <- evd::rmvevd(50, dep = 0.8, model = "log", d = 3, mar = c(1,1,1))
op <- par(mfrow=c(1,3))
Amd <- beed(data, x, 3, "md", "emp", 18, plot=TRUE)
Acfg <- beed(data, x, 3, "cfg", "emp", 18, plot=TRUE)
Aht <- beed(data, x, 3, "ht", "emp", 18, plot=TRUE)
par(op)
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace