P <-matrix(0, nrow(A), nrow(A))rownames(P) <- wcolnames(P) <- w G <-matrix(0, ncol(A), ncol(A))rownames(G) <- ecolnames(G) <- e At <-t(A) s1 <-cbind(P, A) s2 <-cbind(At, G) B <-rbind(s1, s2) B2 <- B%*%B
sim.res <-gen.sim.corr.abs(B) mds.res <-get.mds(sim.res$col.sims) mds.res[, 2] <- mds.res[, 2] *10000 clus.res <-get.clus(mds.res, k =6) p <-get.plot(clus.res) p