Diagonal-restricted likelihood-ratio test
t2diagcov.lr(x1, x2, include.mean = FALSE)
x1
no descr
x2
include.mean
n.stadler