次のコマンドを実行しています。
>simplex(a = a, A2 = A2, b2 = b2, A3 = A3, b3 = b3)
エラーが発生する:
Error in pivot(tableau, prow, pcol) :
NAs are not allowed in subscripted assignments
ここに私の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)
変数 "i" が 5 より大きい場合、上記のエラーが発生します。それ以外の場合は、次のエラーが発生します。
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
これらのエラーが LP 問題を解決できないことを意味するのか、それとも関数の呼び出し方法に誤りがあることを意味するのか、私には理解できません。
前もって感謝します。