1
Je me demande comment je pourrais faire ma fonction Bpp
à accepter un vecteur pour son premier argumentt
?Comment faire `integra()` pour accepter un vecteur dans une fonction R?
Bpp = function(t, n1, n2 = NULL){
N = ifelse(is.null(n2), n1, n1*n2/(n1+n2))
df = ifelse(is.null(n2), n1 - 1, n1 + n2 - 2)
H1 = integrate(function(delta)dcauchy(delta, 0, sqrt(2)/2)*dt(t, df, delta*sqrt(N)), -Inf, Inf)[[1]]
H0 = dt(t, df)
BF10 = H1/H0
p.value = 2*(1-pt(abs(t), df))
list(BF10 = BF10, p.value = p.value)
}
Bpp(t = -6:6, 20, 20) ## This will give error because `t` is now a vector?
'Vectorize' est généralement utile dans ces situations. Bien que réécrire pour accepter directement un vecteur est préférable. – Dason