获得；订阅的工作分配中不允许使用NAs“；在R中使用单纯形时

我正在运行以下命令：

>simplex(a = a, A2 = A2, b2 = b2, A3 = A3, b3 = b3)

我从中得到错误信息：

Error in pivot(tableau, prow, pcol) : 
  NAs are not allowed in subscripted assignments

Error in simplex1(c(a, rep(0, m1 + 2 * m2 + m3)), cbind(rbind(A1, A2,  : 
  number of items to replace is not a multiple of replacement length
In addition: Warning message:
In simplex1(c(a, rep(0, m1 + 2 * m2 + m3)), cbind(rbind(A1, A2,  :
  number of items to replace is not a multiple of replacement length

这是我的R脚本：

library("boot")     # For simplex

set.seed(100)

i <- 6
X <- rnorm(i, mean = 4, sd = 1)
e <- rnorm(i, mean = 0, sd = 1)
Y <- 5*X + e

sum <- 0
Y_a <- matrix(0,1,floor(i/2))                 
for(miss in 1:floor(i/2)){
    sum <- sum + Y[miss*2]
    Y_a[miss] <- Y[miss*2]
}

Y_p <- 5*rnorm(floor(i/2), mean = 4, sd = 1) + rnorm(floor(i/2), mean = 0, sd = 1)          
a <- matrix(0,1,2*floor(i/2))
for(miss in 1:floor(i/2)){
    a[miss] <- 1
}
A3 <- t(apply(a,1, rev))
A2 <- matrix(0,2*floor(i/2),2*floor(i/2))

for(miss in 1:floor(i/2)){
    index <- 2*miss
    index_ <- index - 1
    A2[index_,miss] <- 1
    A2[index_,miss+floor(i/2)] <- 1
    A2[index,miss] <- 1
    A2[index,miss+floor(i/2)] <- -1
}
b2 <- matrix(0,2*floor(i/2),1)
b2[(1:(2*floor(i/2)))%%2==1] <- Y_p
b2[(1:(2*floor(i/2)))%%2==0] <- -Y_p
simplex(a = a, A2 = A2, b2 = b2, A3 = A3, b3 = sum)

我无法理解这些错误是否意味着LP问题无法解决，或者意味着函数调用方式存在错误

---

提前感谢。

关于参数

b2

的单纯形

文档说明

长度为m2的向量，表示>=
限制。如果给定并忽略A2，则需要此参数
否则b2中的所有值必须为非负值。请注意
约束x>=0将自动包括在内，因此不应
这里重复

尝试确保
b2
始终为非负值。
单工
b2
参数的文档说：
b2中的所有值都必须为非负值。
但您有3个负值。我不确定这是否是罪魁祸首，但这是一个开始。谢谢罗曼，似乎是罪魁祸首：）。在我确定b2是正的，并将负的b2转化成方程得到b1之后，我得到了答案。罗曼，你能把它作为答案贴出来吗。