# NOT RUN {
Lambda <- matrix(c(.82, .10, 0, 0,
.77, .35, 0, 0,
.79, .32, 0, 0,
.66, .39, 0, 0,
.51, 0, .71, 0,
.56, 0, .43, 0,
.68, 0, .13, 0,
.60, 0, .50, 0,
.83, 0, 0, .47,
.60, 0, 0, .27,
.78, 0, 0, .28,
.55, 0, 0, .75),
ncol = 4, byrow = TRUE)
colnames(Lambda) <- c("General", "SF1", "SF2", "SF3")
Thresh = list(c(-1, 0, 1), c(-0.5, 0, 0.5),
c(0, 1, 2), c(0, 0.5, 1),
c(-2, -1, 0), c(-1, -0.5, 0),
c(-1, 0, 2), c(-0.5, 0, 1),
c(-2, 0, 1), c(-1, 0, 0.5),
c(-1, 0, 1), c(-0.5, 0, 0.5))
cat_Omega_H(Lambda, Thresh)
# }
Run the code above in your browser using DataLab