ourMELONS/R/computeDiffInCounts.R
2020-11-09 14:58:16 +01:00

17 lines
561 B
R
Raw Blame History

computeDiffInCounts <- function(rows, max_noalle, nloci, data) {
# % Muodostaa max_noalle*nloci taulukon, jossa on niiden alleelien
# % lukum<75><6D>r<EFBFBD>t (vastaavasti kuin COUNTS:issa), jotka ovat data:n
# % riveill<6C> rows. rows pit<69><74> olla vaakavektori.
diffInCounts <- zeros(max_noalle, nloci)
for (i in rows) {
row <- data[i, ]
notEmpty <- find(row>=0)
if (length(notEmpty) > 0) {
diffInCounts(row(notEmpty) + (notEmpty - 1) * max_noalle) <-
diffInCounts(row(notEmpty) + (notEmpty - 1) * max_noalle) + 1
}
}
return(diffInCounts)
}