Я построил следующую функцию для преобразования матрицы данных.
Кажется, мои ошибки зависят от ненулевого минимума, почти в каждом столбце есть ноль в моем data.frame / matrix.
Моя функция следующая:
Tfun = function(x) {log10(x + 10^as.integer(log10(min(x[x > 0])))) - as.integer(log10(min(x[x > 0])))}
Сейчас я не могу передать матрицу с NA, потому что base::min, где na.rm = TRUE не будет работать.
Есть ли альтернатива min(x[x>0]), чтобы я мог затем использовать следующее для создания преобразованного набора данных:
xmNA <- sapply(tapply(m, col(m), na.action = na.exclude), Tfun)