我已经写了这个函数
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
你能告诉我如何规避吗?我希望将此过程推广到数千行,我无法想象手动完成。
谢谢你。
最佳答案
我想你可以改变你的功能...
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]))
我认为这应该都有效。抱歉,当前正在运行一些代码,无法验证/测试它。
关于r - 将具有两个参数的函数应用于 n × 2 矩阵,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/29542130/