Escrevi uma função para estimar o valor de p de correlações por permutação, usando pontos para passar argumentos adicionais para cor
. Aqui uma versão simples para teste:
set.seed(86)
permCorTest <- function(x, y, n = 1e3, ...) {
obs <- cor(x, y, ...)
per <- replicate(n, cor(x, sample(y), ...))
c(cor = obs, p = mean(abs(per) >= abs(obs)))
}
with(mtcars, permCorTest(mpg, qsec, method = "k"))
#> Error in cor(x, sample(y), ...) : invalid 'use' argument
O erro é gerado pelo uso dos pontos em cor
dentro da replicate
. Posso separar as duas para a função funcionar:
permCorTest2 <- function(x, y, n = 1e3, ...) {
obs <- cor(x, y, ...)
yr <- replicate(n, sample(y))
per <- apply(yr, 2, function(s) cor(x, s, ...))
c(cor = obs, p = mean(abs(per) >= abs(obs)))
}
with(mtcars, permCorTest2(mpg, qsec, method = "k"))
#> cor p
#> 0.3153652 0.0180000
Ou usar sapply
ao invés de replicate
. Mas gostaria de entender o que leva a esse erro nessa situação. Tem alguma forma de passar os argumentos adicionais para cor
dentro da replicate
?