Нелинейная функция Optimize / Optimize () для Ordinal Logisti c Функция регрессии в R

Question

Я хочу максимизировать эту нелинейную функцию, но я не могу понять, что я делаю неправильно. Вот код:

max.fun<-function(y0,y1,y2,pi0,pi1,pi2){y0*log(pi0)+y1*log(pi1)+y2*log(pi2)} 

th = matrix(c(-1,-3,1),nrow=1,ncol = 3)
y0 = matrix(c(26,104),nrow=1,ncol = 2)
y1 = matrix(c(33,72),nrow=1,ncol = 2)
y2 = matrix(c(22,31),nrow=1,ncol = 2)
race = matrix(c(1,0),nrow=1,ncol = 2)

alpha1<-th[1]
alpha2<-th[2]
beta<-th[3]

ex1<-exp(alpha1+beta*race)
ex2<-exp(alpha2+beta*race)

ee1<-ex1/(1+ex1)
ee2<-ex2/(1+ex2)

pi0<-1/(1+ex1)
pi1<-ee1-ee2
pi2<-ee2


optimize(max.fun,
         c(0,1),
         y0=y0,
         y1=y1,
         y2=y2,
         pi0=pi0,
         pi1=pi1,
         pi2=pi2,
         tol = 0.0001)

Когда я запускаю функцию оптимизации, я получаю это сообщение об ошибке:

Error in f(arg, ...) : unused argument (arg)

Любая помощь очень ценится.

Answer 1

Что вы пытаетесь оптимизировать? Если x, то он должен быть указан в LHS, и должны быть предоставлены другие значения, например, если вы хотите максимизировать y0:

max.fun<-function(y0,y1,y2,pi0,pi1,pi2){y0*log(pi0)+y1*log(pi1)+y2*log(pi2)} 

th = matrix(c(-1,-3,1),nrow=1,ncol = 3)
y0 = matrix(c(26,104),nrow=1,ncol = 2)
y1 = matrix(c(33,72),nrow=1,ncol = 2)
y2 = matrix(c(22,31),nrow=1,ncol = 2)
race = matrix(c(1,0),nrow=1,ncol = 2)

alpha1<-th[1]
alpha2<-th[2]
beta<-th[3]

ex1<-exp(alpha1+beta*race)
ex2<-exp(alpha2+beta*race)

ee1<-ex1/(1+ex1)
ee2<-ex2/(1+ex2)

pi0<-1/(1+ex1)
pi1<-ee1-ee2
pi2<-ee2

# Using purely illustrative values
y0m <- optimize(max.fun,
         c(0,1),
         y1=2,
         y2=3,
         pi0=4,
         pi1=5,
         pi2=6,
         tol = 0.0001,
         maximum = T)
y0m

> y0m
$maximum
[1] 0.9999339

$objective
[1] 9.980357

Если этот / любой другой предоставленный вам ответ полезен для вас, примите, пожалуйста, ответить / проголосовать +!