library(copula)
library(missMethods)
n = 100
cp = claytonCopula(param = c(1), dim = 5)
P = mvdc(copula = cp, margins = c("exp", "exp", "exp", "exp", "exp"),
paramMargins = list(list(1), list(1), list(1), list(1), list(1)))
X = rMvdc(n, P)
X = delete_MCAR(X, 0.1, c(1,4,5))
get_SigmaS(X)
get_SigmaS(X, min_diff = 20)
Run the code above in your browser using DataLab