z = bn.moments(learning.test, algorithm = "gs", R = 100)
bn.var(z, method = "tvar")
# statistic normalized
# 1.29060 0.34416
bn.var.test(z, method = "nvar")
#
# Squared Frobenius Norm
#
# data: covariance matrix
# nvar = 0.5471, B = 5000, R = 100, p-value < 2.2e-16
# alternative hypothesis: true value is greater than 0Run the code above in your browser using DataLab