Решение нелинейного уравнения с помощью итерации Ньютона-Рафсона в R

Я биолог и пытаюсь решить следующее нелинейное уравнение:

R^2*exp(c*R) =Ccf*Acf*E*(I/Ke+I)

где c, Ccf, Acf,E и Ke — фиксированные значения, а I — вектор со значениями освещенности.

Я предполагаю, что итерация Ньютона-Рафсона будет подходящим способом, но я не могу заставить ее работать ни с мультирутом {rootSolve}, ни с newtonraphson {spuRs}.

Любая помощь будет принята с благодарностью! Заранее спасибо!


r iteration exp solver nonlinear-functions
person Larusson    schedule 19.08.2015    source источник
comment
Можете показать, что вы уже пробовали? (т.е. воспроизводимый пример)   -  person Heroka    schedule 19.08.2015
comment
Надеюсь, что это дает и идею...library(spuRs) getR <- function (R) { Ccf = c(0.3) Acf = c(3*10^-6) E = c(0.00006) Ke = c(1000) I = c(150) c = c(0.1) R^2*exp(c*R) =Ccf*Acf*E*(I/Ke+I) } newtonraphson(getR, 2, 1e-6)   -  person Larusson    schedule 19.08.2015
comment
Здесь - это то, что wolframalpha.com предлагает для уравнения этого типа.   -  person RHertel    schedule 19.08.2015
comment
Я знаком с wolframalpha.com, но еще не понял, как использовать его в R.   -  person Larusson    schedule 19.08.2015
comment
Для полноты ссылка на решение, которое я, наконец, внедрил в свою модель function/32167555#32167555" title="использовать значения из разных фреймов данных в качестве входных данных для вложенной функции"> stackoverflow.com/questions/32167058/   -  person Larusson    schedule 23.08.2015