data(galo, package = "homals")
galo <- galo[, -(ncol(galo))]
ngalo <- nrow(galo)
mgalo <- ncol(galo)
kgalo <- sapply(galo, function(x) length(unique(x)))
ggalo <- NULL
for (j in 1:4) {
  ggalo <- cbind(ggalo, ifelse (outer (galo[[j]], unique(galo[[j]]), "=="), 1, 0))
}
cgalo <- crossprod(ggalo)
dgalo <- diag(cgalo)
egalo <- cgalo / sqrt(outer(dgalo, dgalo))
sgalo <- eigen(egalo)
vgalo <- sgalo$vectors
lgalo <- sgalo$values / mgalo
ygalo <- (1 / sqrt(dgalo)) * vgalo
hgalo <- DMCA(cgalo, kgalo, itmax = 100, verbose = TRUE)
wgalo <- diag(hgalo$lpkekpl) / mgalo
zgalo <- (1 / sqrt(dgalo)) * hgalo$kpl
rgalo <- crossprod (vgalo, hgalo$kpl)
pgalo1 <- hgalo$pkekp[1:4, 1:4]
pgalo2 <- hgalo$pkekp[5:8, 5:8]
pgalo3 <- hgalo$pkekp[9:11, 9:11]
pgalo4 <- hgalo$pkekp[12:14, 12:14]
pgalo5 <- hgalo$pkekp[15:17, 15:17]
pgalo6 <- hgalo$pkekp[18:20, 18:20]
pgalo7 <- hgalo$pkekp[21:22, 21:22]
pgalo8 <- hgalo$pkekp[23:24, 23:24]