R/bpbounds_tri_x2y2z3.R

Defines functions bpbounds_tri_x2y2z3

A_tri_x2y2z3 <- matrix(nrow = 117, ncol = 12)
A_tri_x2y2z3[1, ] <- c(0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
A_tri_x2y2z3[2, ] <- c(0, 0, 0, 0, 0, -1, -1, 0, 0, 0, 0, 0)
A_tri_x2y2z3[3, ] <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, -1, -1, 0)
A_tri_x2y2z3[4, ] <- c(0, -1, -1, -1, 0, 0, 0, 0, 0, 1, 0, 0)
A_tri_x2y2z3[5, ] <- c(0, -1, -1, -1, 0, 1, 0, 0, 0, 0, 0, 0)
A_tri_x2y2z3[6, ] <- c(0, 0, 0, 0, 0, -1, -1, -1, 0, 1, 0, 0)
A_tri_x2y2z3[7, ] <- c(0, 0, 0, 0, 0, 1, 0, 0, 0, -1, -1, -1)
A_tri_x2y2z3[8, ] <- c(0, 1, 0, 0, 0, -1, -1, -1, 0, 0, 0, 0)
A_tri_x2y2z3[9, ] <- c(0, 1, 0, 0, 0, 0, 0, 0, 0, -1, -1, -1)
A_tri_x2y2z3[10, ] <- c(0, -1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 1)
A_tri_x2y2z3[11, ] <- c(0, -1, -1, -1, 0, 0, 0, 1, 0, 0, 0, 0)
A_tri_x2y2z3[12, ] <- c(0, -1, -1, 0, 0, -1, 0, -1, 0, 1, 0, 0)
A_tri_x2y2z3[13, ] <- c(0, -1, -1, 0, 0, 1, 0, 0, 0, -1, 0, -1)
A_tri_x2y2z3[14, ] <- c(0, -1, 0, -1, 0, -1, -1, 0, 0, 1, 0, 0)
A_tri_x2y2z3[15, ] <- c(0, -1, 0, -1, 0, 1, 0, 0, 0, -1, -1, 0)
A_tri_x2y2z3[16, ] <- c(0, 0, 0, 0, 0, -1, -1, -1, 0, 0, 0, 1)
A_tri_x2y2z3[17, ] <- c(0, 0, 0, 0, 0, 0, 0, 1, 0, -1, -1, -1)
A_tri_x2y2z3[18, ] <- c(0, 0, 0, 1, 0, -1, -1, -1, 0, 0, 0, 0)
A_tri_x2y2z3[19, ] <- c(0, 0, 0, 1, 0, 0, 0, 0, 0, -1, -1, -1)
A_tri_x2y2z3[20, ] <- c(0, 1, 0, 0, 0, -1, -1, 0, 0, -1, 0, -1)
A_tri_x2y2z3[21, ] <- c(0, 1, 0, 0, 0, -1, 0, -1, 0, -1, -1, 0)
A_tri_x2y2z3[22, ] <- c(0, -1, -1, 0, 0, -1, 0, -1, 0, 0, 0, 1)
A_tri_x2y2z3[23, ] <- c(0, -1, -1, 0, 0, 0, 0, 1, 0, -1, 0, -1)
A_tri_x2y2z3[24, ] <- c(0, -1, -1, 1, 0, -1, 0, -1, 0, 0, 0, 0)
A_tri_x2y2z3[25, ] <- c(0, -1, -1, 1, 0, 0, 0, 0, 0, -1, 0, -1)
A_tri_x2y2z3[26, ] <- c(0, -1, 0, -1, 0, -1, -1, 0, 0, 0, 0, 1)
A_tri_x2y2z3[27, ] <- c(0, -1, 0, -1, 0, -1, -1, 1, 0, 0, 0, 0)
A_tri_x2y2z3[28, ] <- c(0, -1, 0, -1, 0, 0, 0, 0, 0, -1, -1, 1)
A_tri_x2y2z3[29, ] <- c(0, -1, 0, -1, 0, 0, 0, 1, 0, -1, -1, 0)
A_tri_x2y2z3[30, ] <- c(0, 0, 0, 0, 0, -1, -1, 1, 0, -1, 0, -1)
A_tri_x2y2z3[31, ] <- c(0, 0, 0, 0, 0, -1, 0, -1, 0, -1, -1, 1)
A_tri_x2y2z3[32, ] <- c(0, 0, 0, 1, 0, -1, -1, 0, 0, -1, 0, -1)
A_tri_x2y2z3[33, ] <- c(0, 0, 0, 1, 0, -1, 0, -1, 0, -1, -1, 0)
A_tri_x2y2z3[34, ] <- c(0, -1, -1, -1, 0, -1, -1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[35, ] <- c(0, -1, -1, -1, 0, 1, 0, 1, 0, -1, -1, 0)
A_tri_x2y2z3[36, ] <- c(0, -1, -1, 0, 0, -1, -1, -1, 0, 1, 0, 1)
A_tri_x2y2z3[37, ] <- c(0, -1, -1, 0, 0, 1, 0, 1, 0, -1, -1, -1)
A_tri_x2y2z3[38, ] <- c(0, 1, 0, 1, 0, -1, -1, -1, 0, -1, -1, 0)
A_tri_x2y2z3[39, ] <- c(0, 1, 0, 1, 0, -1, -1, 0, 0, -1, -1, -1)
A_tri_x2y2z3[40, ] <- c(0, -2, -2, -1, 0, 0, 0, 0, 0, 1, 0, 1)
A_tri_x2y2z3[41, ] <- c(0, -2, -2, -1, 0, 1, 0, 1, 0, 0, 0, 0)
A_tri_x2y2z3[42, ] <- c(0, 0, 0, 0, 0, -2, -2, -1, 0, 1, 0, 1)
A_tri_x2y2z3[43, ] <- c(0, 0, 0, 0, 0, 1, 0, 1, 0, -2, -2, -1)
A_tri_x2y2z3[44, ] <- c(0, 1, 0, 1, 0, -2, -2, -1, 0, 0, 0, 0)
A_tri_x2y2z3[45, ] <- c(0, 1, 0, 1, 0, 0, 0, 0, 0, -2, -2, -1)
A_tri_x2y2z3[46, ] <- c(0, -2, -2, 0, 0, -1, 0, -1, 0, 1, 0, 1)
A_tri_x2y2z3[47, ] <- c(0, -2, -2, 0, 0, 1, 0, 1, 0, -1, 0, -1)
A_tri_x2y2z3[48, ] <- c(0, -1, 0, -1, 0, -2, -2, 0, 0, 1, 0, 1)
A_tri_x2y2z3[49, ] <- c(0, -1, 0, -1, 0, 1, 0, 1, 0, -2, -2, 0)
A_tri_x2y2z3[50, ] <- c(0, 1, 0, 1, 0, -2, -2, 0, 0, -1, 0, -1)
A_tri_x2y2z3[51, ] <- c(0, 1, 0, 1, 0, -1, 0, -1, 0, -2, -2, 0)
A_tri_x2y2z3[52, ] <- c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0)
A_tri_x2y2z3[53, ] <- c(0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1)
A_tri_x2y2z3[54, ] <- c(0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0)
A_tri_x2y2z3[55, ] <- c(0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0)
A_tri_x2y2z3[56, ] <- c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1)
A_tri_x2y2z3[57, ] <- c(0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0)
A_tri_x2y2z3[58, ] <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0)
A_tri_x2y2z3[59, ] <- c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1)
A_tri_x2y2z3[60, ] <- c(0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0)
A_tri_x2y2z3[61, ] <- c(0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0)
A_tri_x2y2z3[62, ] <- c(0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0)
A_tri_x2y2z3[63, ] <- c(0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0)
A_tri_x2y2z3[64, ] <- c(0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0)
A_tri_x2y2z3[65, ] <- c(0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0)
A_tri_x2y2z3[66, ] <- c(0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0)
A_tri_x2y2z3[67, ] <- c(0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0)
A_tri_x2y2z3[68, ] <- c(0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
A_tri_x2y2z3[69, ] <- c(0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0)
A_tri_x2y2z3[70, ] <- c(0, -1, 0, -1, 0, 0, 0, 0, 0, 2, 1, 0)
A_tri_x2y2z3[71, ] <- c(0, -1, 0, -1, 0, 2, 1, 0, 0, 0, 0, 0)
A_tri_x2y2z3[72, ] <- c(0, 0, 0, 0, 0, -1, 0, -1, 0, 2, 1, 0)
A_tri_x2y2z3[73, ] <- c(0, 0, 0, 0, 0, 2, 1, 0, 0, -1, 0, -1)
A_tri_x2y2z3[74, ] <- c(0, 2, 1, 0, 0, -1, 0, -1, 0, 0, 0, 0)
A_tri_x2y2z3[75, ] <- c(0, 2, 1, 0, 0, 0, 0, 0, 0, -1, 0, -1)
A_tri_x2y2z3[76, ] <- c(0, -1, -1, 0, 0, 0, 1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[77, ] <- c(0, -1, -1, 0, 0, 1, 0, 1, 0, 0, 1, 0)
A_tri_x2y2z3[78, ] <- c(0, -1, 0, -1, 0, 0, 0, 1, 0, 1, 1, 0)
A_tri_x2y2z3[79, ] <- c(0, -1, 0, -1, 0, 1, 1, 0, 0, 0, 0, 1)
A_tri_x2y2z3[80, ] <- c(0, 0, 0, 1, 0, -1, 0, -1, 0, 1, 1, 0)
A_tri_x2y2z3[81, ] <- c(0, 0, 0, 1, 0, 1, 1, 0, 0, -1, 0, -1)
A_tri_x2y2z3[82, ] <- c(0, 0, 1, 0, 0, -1, -1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[83, ] <- c(0, 0, 1, 0, 0, 1, 0, 1, 0, -1, -1, 0)
A_tri_x2y2z3[84, ] <- c(0, 1, 0, 1, 0, -1, -1, 0, 0, 0, 1, 0)
A_tri_x2y2z3[85, ] <- c(0, 1, 0, 1, 0, 0, 1, 0, 0, -1, -1, 0)
A_tri_x2y2z3[86, ] <- c(0, 1, 1, 0, 0, -1, 0, -1, 0, 0, 0, 1)
A_tri_x2y2z3[87, ] <- c(0, 1, 1, 0, 0, 0, 0, 1, 0, -1, 0, -1)
A_tri_x2y2z3[88, ] <- c(0, -1, 0, -1, 0, 1, 0, 0, 0, 1, 1, 0)
A_tri_x2y2z3[89, ] <- c(0, -1, 0, -1, 0, 1, 1, 0, 0, 1, 0, 0)
A_tri_x2y2z3[90, ] <- c(0, 1, 0, 0, 0, -1, 0, -1, 0, 1, 1, 0)
A_tri_x2y2z3[91, ] <- c(0, 1, 0, 0, 0, 1, 1, 0, 0, -1, 0, -1)
A_tri_x2y2z3[92, ] <- c(0, 1, 1, 0, 0, -1, 0, -1, 0, 1, 0, 0)
A_tri_x2y2z3[93, ] <- c(0, 1, 1, 0, 0, 1, 0, 0, 0, -1, 0, -1)
A_tri_x2y2z3[94, ] <- c(0, -1, -1, -1, 0, 1, 0, 1, 0, 1, 1, 0)
A_tri_x2y2z3[95, ] <- c(0, -1, -1, -1, 0, 1, 1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[96, ] <- c(0, 1, 0, 1, 0, -1, -1, -1, 0, 1, 1, 0)
A_tri_x2y2z3[97, ] <- c(0, 1, 0, 1, 0, 1, 1, 0, 0, -1, -1, -1)
A_tri_x2y2z3[98, ] <- c(0, 1, 1, 0, 0, -1, -1, -1, 0, 1, 0, 1)
A_tri_x2y2z3[99, ] <- c(0, 1, 1, 0, 0, 1, 0, 1, 0, -1, -1, -1)
A_tri_x2y2z3[100, ] <- c(0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 2, 0)
A_tri_x2y2z3[101, ] <- c(0, 0, 0, 0, 0, 1, 2, 0, 0, 1, 0, 1)
A_tri_x2y2z3[102, ] <- c(0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 2, 0)
A_tri_x2y2z3[103, ] <- c(0, 1, 0, 1, 0, 1, 2, 0, 0, 0, 0, 0)
A_tri_x2y2z3[104, ] <- c(0, 1, 2, 0, 0, 0, 0, 0, 0, 1, 0, 1)
A_tri_x2y2z3[105, ] <- c(0, 1, 2, 0, 0, 1, 0, 1, 0, 0, 0, 0)
A_tri_x2y2z3[106, ] <- c(0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0)
A_tri_x2y2z3[107, ] <- c(0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[108, ] <- c(0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0)
A_tri_x2y2z3[109, ] <- c(0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0)
A_tri_x2y2z3[110, ] <- c(0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1)
A_tri_x2y2z3[111, ] <- c(0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0)
A_tri_x2y2z3[112, ] <- c(0, -1, 0, -1, 0, 1, 0, 1, 0, 2, 2, 0)
A_tri_x2y2z3[113, ] <- c(0, -1, 0, -1, 0, 2, 2, 0, 0, 1, 0, 1)
A_tri_x2y2z3[114, ] <- c(0, 1, 0, 1, 0, -1, 0, -1, 0, 2, 2, 0)
A_tri_x2y2z3[115, ] <- c(0, 1, 0, 1, 0, 2, 2, 0, 0, -1, 0, -1)
A_tri_x2y2z3[116, ] <- c(0, 2, 2, 0, 0, -1, 0, -1, 0, 1, 0, 1)
A_tri_x2y2z3[117, ] <- c(0, 2, 2, 0, 0, 1, 0, 1, 0, -1, 0, -1)

ice_tri_x2y2z3 = c(
  -1,
  -1,
  -1,
  0,
  0,
  0,
  0,
  0,
  0,
  -1,
  -1,
  0,
  0,
  0,
  0,
  -1,
  -1,
  -1,
  -1,
  0,
  0,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  0,
  0,
  0,
  0,
  0,
  0,
  1,
  0,
  1,
  1,
  1,
  0,
  1,
  1,
  1,
  1,
  1,
  0,
  1,
  1,
  1,
  1,
  1,
  1,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  1,
  1,
  1,
  1,
  1,
  1,
  0,
  0,
  0,
  0,
  0,
  0,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1,
  1
)

cons_tri_x2y2z3 = c(
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  0,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -1,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2,
  -2
)

bpbounds_tri_x2y2z3 <- function(p,
                                A = A_tri_x2y2z3,
                                ice = ice_tri_x2y2z3,
                                cons = cons_tri_x2y2z3) {
  prod = A %*% p
  prod = prod + cons
  ivinequality <- prod[ice == 0]
  low <- prod[ice == -1]
  upp <- -1 * prod[ice == 1]
  inequality <- max(ivinequality) <= 0
  bplow <- max(low)
  bpupp <- min(upp)

  return(
    list(
      "inequality" = inequality,
      "bplow" = bplow,
      "bpupp" = bpupp,
      "bplower" = low,
      "bpupper" = upp
    )
  )
}

Try the bpbounds package in your browser

Any scripts or data that you put into this service are public.

bpbounds documentation built on Jan. 21, 2020, 5:07 p.m.