r - 将具有两个参数的函数应用于 n × 2 矩阵

Question

我已经写了这个函数

tippett <- function(p1, p2) {
  t <- 1 - (1 - min(p1, p2))^2
  return(list(p.value=t))
}

我想将它应用在水平生成数字的矩阵中。例如在这个矩阵中: z <- matrix(c(rbeta(5, 1, 5), rbeta(5, 1, 10)), ncol=2) .

不幸的是 apply似乎不起作用并以这种方式使用它 s<-apply(z,1,tippett)
产生这个错误

Error in FUN(newX[, i], ...) : argument "p2" is missing, with no default

你能告诉我如何规避吗？我希望将此过程推广到数千行，我无法想象手动完成。

谢谢你。

Answer 1

我想你可以改变你的功能...

tippett<-function(p){
t<-1-(1-min(p))^2
return(list(p.value=t))
}

那么你的函数调用应该可以工作。 (虽然我没有测试过这个......)

s<-apply(z,1,tippett)

或者，如果您无法更改此功能，您可以尝试类似...

s <- sapply(seq_along(nrow(z)), function(i) tippett(z[i,1], z[i,2]))

我认为这应该都有效。抱歉，当前正在运行一些代码，无法验证/测试它。