Как решить два уравнения c с двумя символами с нормальным CDF (pnorm) с использованием R - PullRequest
Новая
       60

Как решить два уравнения c с двумя символами с нормальным CDF (pnorm) с использованием R

2 голосов
/ 05 мая 2020

Я хочу решить эти два уравнения численно, используя R

Уравнение 1: pnorm ((c2 + 1) / 5) -pnorm ((c1 + 1) / 5) = 0,025

Уравнение 2: pnorm ((c2-1) / 5) -pnorm ((c1-1) / 5) = 0,025

Мне нужно найти значения для c1 и c2

I пытался использовать rSymPy и Ryacas, но не мог понять, как это сделать.

Код R Я пробовал

решить ([Eq (pnorm ((c2 + 1) / 5) - pnorm ((c1 + 1) / 5), 0,025), Eq (pnorm ((c2-1) / 5) -pnorm ((c1-1) / 5), 0,02)]

1 Ответ

3 голосов
/ 05 мая 2020

Вы смотрели пакет nleqslv в R? Используя его, вы легко можете получить решение. Решение не сработает, потому что у вас есть нелинейная система уравнений, которую нельзя представить матрицами. См. Ниже описание вашей проблемы: 

library(nleqslv)


eqs_fun <- function(x) {
  y = numeric(2)
  y[1] = pnorm(q = (x[2]+1)/5) - pnorm(q = (x[1]+1)/5) - .025
  y[2] = pnorm(q = (x[2]-1)/5) - pnorm(q = (x[1]-1)/5) - .025
  y
}

xstart = c(0,1)
nleqslv(xstart, eqs_fun)

Запуск этого должен дать вам следующий результат: 

list(x = c(-0.159856055122849, 0.159854416799285), fvec = c(1.006390991376e-09, 
-6.31374431903087e-10), termcd = 1L, message = "Function criterion near zero", 
    scalex = c(1, 1), nfcnt = 8L, njcnt = 1L, iter = 6L)

, и вы можете сами убедиться, что c1 = -.16, c2 = .16 будет дать вам решение!

...