When groupBy input is a factor class, the table() output in Line 821 will include groups that has 0 frequencies, resulting in different vector lengths in estbgd (only non-zero groups) and obsbgd, and other objects derived from these.
|
estbgd <- sort(floor(minTotal * bgdProbx)) |
|
obsbgd <- rep(0, length(estbgd)) |
|
names(obsbgd) <- names(estbgd) |
|
tabGroups <- table(groups[idY]) |
|
obsbgd[names(tabGroups)] <- tabGroups |
|
estbgdP <- round(100 * estbgd / sum(estbgd),3) |
|
obsbgdP <- round(100 * obsbgd / sum(obsbgd),3) |
During the "Matching Known Biases" step, we can see these warning messages when ArchR is creating the out object in the following lines.
Warning message in rbind(estbgd, obsbgd):
“number of columns of result is not a multiple of vector length (arg 1)”
Warning message in rbind(estbgdP, obsbgdP):
“number of columns of result is not a multiple of vector length (arg 1)”
Due to the different lengths of estbgdP and obsbgdP, the correlation calculation in Line 847 fail to execute.
|
.logThis(out, paste0("MatchSummary : Pre", useGroups[x]), logFile = logFile) |
|
|
|
out$corbgdGroups <- suppressWarnings(cor(estbgdP, obsbgdP)) |
Error in cor(estbgdP, obsbgdP) : incompatible dimensions
This has been reported in these issues in ArchR repo:
When
groupByinput is a factor class, thetable()output in Line 821 will include groups that has 0 frequencies, resulting in different vector lengths inestbgd(only non-zero groups) andobsbgd, and other objects derived from these.ArchR/R/MarkerFeatures.R
Lines 818 to 824 in 8b7c4bc
During the "Matching Known Biases" step, we can see these warning messages when ArchR is creating the
outobject in the following lines.Due to the different lengths of
estbgdPandobsbgdP, the correlation calculation in Line 847 fail to execute.ArchR/R/MarkerFeatures.R
Lines 845 to 847 in 8b7c4bc
This has been reported in these issues in ArchR repo: