Как избежать ошибки uniroot, которая останавливает цикл [дубликат]

Question

Я выполняю функцию в цикле, но получаю ошибку, и код останавливается. Код приведен ниже;

      func <-function(f) -a*b/c*0.5*d*e^2 + (d/f-1)*g*sin(h*(pi/180))-i
dat <- data.frame(a = c(0.99,0.99,0.99),
                  b = c(0.1986572,0.1986572,0.1986572),
                  c = c(237.5,237.5,237.5),
                  d = c(1028.372, 1028.711, 1028.372),
                  e = c(2.46261, 2.986461, 2.46261),
                  f = c(-1,-1,-1),
                  g = c(9.8,9.8,9.8),
                  h = c(-54.97964, -51.65978, -54.97964),
                  i = c(0.03699588, -0.0375189, 0.03699588))

for(j in 1:length(dat$a)){
   a <- dat$a[j]
   b <- dat$b[j]
   c <- dat$c[j]
   d <- dat$d[j]
   e <- dat$e[j]
   #f: this should be solved by uniroot
   g <- dat$g[j]
   h <- dat$h[j]
   i <- dat$i[j]
   sol <- uniroot(func,c(0, 2000),extendInt = "yes") 
   dat$f[j] <- sol$root
   print(j)
}

Запустив приведенный выше код, нажмите на ошибку ниже:

      [1] 1
Error in uniroot(func, c(0, 2000), extendInt = "yes") : 
      no sign change found in 1000 iterations

Код остановился на j = 1 и не перешел на j = 2 & 3. Следовательно, dat$f показывает

      > dat$f
[1] 1526.566   -1.000   -1.000

Моя цель - когда uniroot обнаруживает ошибку в заданном j, положил NA в dat$f[j], и продолжаем цикл до конца.

Если это сработает, dat$f[1] а также dat$f[3] должен иметь такое же значение (= 1526,566) с использованием вышеуказанного фрейма данных.

Пожалуйста, посоветуйте мне, как бороться с ошибкой uniroot.

(Это сообщение было изменено, чтобы включить в него код для воспроизведения для всех, как было предложено в моем последнем сообщении )

Answer 1

Попробуйте расширить диапазон интервала. Например:

      sol <- uniroot(func, c(0, 5000), extendInt = "yes") 

Ваш код останавливается, потому что в пределах заданного вами диапазона unirootне нашел решения