从R中的完全矢量化函数子集？

我写了Vectorize d “主”功能下面list的三个对象：BF10,p.value和d。从R中的完全矢量化函数子集？

我不知道节约“主”函数调用像b对象（见下文）时，那么，为什么我不能子集BF10，从bp.value或d？

也就是说，运行主要功能后，当我当我运行：

b = BF.d.pvalue(t = c(2.46, 3.21), n1 = c(20, 30), n2 = c(20, NA)) 
b$BF10 ; b$p.value ; b$d ## None of these return anything !!!

这里是 “主要” 功能：

BF.d.pvalue = Vectorize(function(t, n1, n2 = NA, scale = sqrt(2)/2){ 

     options(warn = -1) 
     N = ifelse(is.na(n2), n1, (n1*n2)/(n1+n2)) 
    df = ifelse(is.na(n2), n1 - 1, (n1 + n2) - 2) 
     d = t/sqrt(N) 

    H1 = integrate(function(delta)dcauchy(delta, 0, scale)*dt(t, df, delta*sqrt(N)), -Inf, Inf)[[1]] 
    H0 = dt(t, df) 
    BF10 = H1/H0 
p.value = 2*(1-pt(abs(t), df)) 

cbind(BF10 = BF10, p.value = p.value, d = d) 

}, vectorize.args = c("t", "n1", "n2", "scale"))

来源

2017-06-25 rnorouzian

mapply默认完成了简化。这里的结果已经简化为一个矩阵。你可以做

b["BF10", ] 
b["p.value", ] 
b["d", ]

虽然你可以设置SIMPLIFY = FALSE，但结果是一个嵌套列表。

如果你想要的东西就像一个列表/数据帧在那里你可以使用$，怎么样

b <- data.frame(t(b))

来源

2017-06-25 03:45:05

从R中的完全矢量化函数子集？

回答

相关问题