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))