Impostazione dei vincoli in constrOptim

Question

Esiste un metodo semplice per impostare theta, ui, ci per i seguenti vincoli nella funzione constrOptim?

c1<x1<=c2 
x1+1<x2<=c2+1 
x2+1<x3<=c2+2 
x3+1<x4<=c2+3 

Ho considerato l'utilizzo di simplex ma richiede solo 3 vincoli.

Grazie

Answer 1

Basta riscrivere i vincoli nella forma desiderata, ui %*% theta >= ci.

# Initial formulation of the constraints 
c1 <= x1 
     x1 <= c2 
x1+1 <= x2 
     x2 <= c2+1 
x2+1 <= x3 
     x3 <= c2+2 
x3+1 <= x4 
     x4 <= c2+3 

# Rewrite them 
    x1    >= c1 
- x1    >= -c2 
- x1 + x2   >= 1 
    - x2   >= -c2 - 1 
    - x2 + x3  >= 1 
      - x3  >= -c2 - 2 
      - x3 + x4 >= 1 
       - x4 >= -c2 - 3 

# In matrix form 
ui <- matrix(c(
    1, 0, 0, 0, 
    -1, 0, 0, 0, 
    -1, 1, 0, 0, 
    0, -1, 0, 0, 
    0, -1, 1, 0, 
    0, 0, -1, 0, 
    0, 0, -1, 1, 
    0, 0, 0, -1 
), 
    ncol = 4, 
    byrow = TRUE 
) 
ci <- c(c1, -c2, 1, -c2-1, 1, -c2-2, 1, -c2-3)