Я не могу понять, как применить функцию HWf из пакета HardyWeinberg к «двойному» объекту или матрице данных. Функция занимает только одну строку с тремя столбцами, но я бы хотел применить ее к каждой строке по очереди с вектором с тремя столбцами. Я хотел бы, чтобы выходные данные были вектором или столбцом результатов HWf (обычно каждый является десятичным, например, 0,25). Вот как далеко я дошел, но, похоже, это не работает. Любая помощь будет высоко ценится!
library(HardyWeinberg)
a <- data.frame(AA=c(1,1,2,1), AB=c(2,3,5,4), BB=c(5,6,8,7))
b <- data.matrix(a)
c <- lapply(b, function(x) HWf(x))