n <- 100
d <- 4
p <- .8
Gamma <- cbind(c(0, 1.5, 1.5, 2),
c(1.5, 0, 2, 1.5),
c(1.5, 2, 0, 1.5),
c(2, 1.5, 1.5, 0))
set.seed(123)
my_data = rmstable(n, "HR", d = d, par = Gamma)
emp_chi_mat(my_data, p)
Run the code above in your browser using DataCamp Workspace