H.6 <- matrix(rbind(diag(3), c(0, 0, 0), c(0, 0, 0)),
              nrow=5, ncol=3)
H6 <- summary(bh6lrtest(z = H1, H = H.6,
                        r = 2, r1 = 1))