Correlation of the red probes, first four eigengenes with the batch effect:
kruskal.test(mredsvd$v[,2],as.factor(batch)) Kruskal-Wallis rank sum test data: mredsvd$v[, 2] and as.factor(batch) Kruskal-Wallis chi-squared = 10.8903, df = 12, p-value = 0.5383 kruskal.test(mredsvd$v[,3],as.factor(batch)) Kruskal-Wallis rank sum test data: mredsvd$v[, 3] and as.factor(batch) Kruskal-Wallis chi-squared = 21.7447, df = 12, p-value = 0.04048 kruskal.test(mredsvd$v[,4],as.factor(batch)) Kruskal-Wallis rank sum test data: mredsvd$v[, 4] and as.factor(batch) Kruskal-Wallis chi-squared = 35.8388, df = 12, p-value = 0.0003439
After removing the batch effect p-value for the association of the first principal component with the batch is 0.7748.