Итак, мне нужно что-то вроде microbenchmark
, но вместо того, чтобы рассчитывать свой код, я хочу знать, как часто он терпит неудачу (если вообще). По сути, я хочу протестировать что-то вроде
f1 <- function(x){
x^2 - 1
}
lims <- rnorm(2)
try(uniroot(f1, lims))
большое количество раз (скажем, 10 ^ 6 тестов), а затем записать частоту, с которой это не удается. Конечно, в этом простом примере все сводится к одному из 1 или -1 в lims
. В общем, я этого не знаю.
Мне нужно решить многочлен как часть блестящего приложения. Я хочу, чтобы приложение было как можно более чистым, поэтому я хочу, чтобы пользователь не указывал ограничения uniroot
. Поэтому я написал адаптивный метод, несколько симуляций показывают, что он работает, но я хочу протестировать его много раз.