output of ld is not symmetric ? #35
Description
trying to improve sparse matrix use in the package, I add the following error in tests:
error: snpClust gives results identical to those of adjclust 0.3.0
Input matrix is not symmetric
1: expect_message(snpClust(ld.ceph, h = 100), "Note: forcing the diagonal of the LD similarity matrix to be 1", all = FALSE) at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/tests/testthat/test_ascendingCompatibility.R:35
2: quasi_capture(enquo(object), capture_messages, label = label)
3: capture(act$val <- eval_bare(get_expr(quo), get_env(quo)))
4: withCallingHandlers(code, message = function(condition) {
out$push(condition)
invokeRestart("muffleMessage")
})
5: eval_bare(get_expr(quo), get_env(quo))
6: snpClust(ld.ceph, h = 100)
7: snpClust.dgCMatrix(ld.ceph, h = 100) at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/R/snpClust.R:78
8: run.snpClust(x, h = h, stats = stats) at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/R/snpClust.R:96
9: adjClust(x, type = "similarity", h = h) at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/R/snpClust.R:140
10: adjClust.dgCMatrix(x, type = "similarity", h = h) at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/R/adjclust.R:81
11: stop("Input matrix is not symmetric") at /home/nathalie/Private/Travail/Recherche/HiC/adjclust/R/adjclust.R:131