A1 <- matrix(c(1,0,0,0,0, 0,0,1,0,0, 0,0,0,1,0, 0,0,0,0,1), c(5,4))
A2 <- matrix(c(1,0,0,0,0, 0,1,0,0,0, 0,0,1,0,0, 0,0,0,1,0), c(5,4))
H41 <- summary(alrtest(z=H1, A=A1, r=2))
H42 <- summary(alrtest(z=H1, A=A2, r=2))