データセットの各行で関数を呼び出す関数を作成しました。ベクトルとして出力したいと思います。以下に示すように、関数は結果を画面に出力しますが、関数の外部で使用できるベクトルに出力をリダイレクトする方法がわかりません。
n_markers <- nrow(data)
p_values <-rep(0, n_markers)
test_markers <- function()
{
for (i in 1:n_markers)
{
hets <- data[i, 2]
hom_1 <- data[i, 3]
hom_2 <- data[i, 4]
p_values[i] <- SNPHWE(hets, hom_1, hom_2)
}
return(p_values)
}
test_markers()