I was solving a LGP problem which is very basic.
Find x0 = [x1; x2], n0 = [n1; n2; n3] and p0 = [p1; p2; p3] that minimize a
= [(2p1); (n2); (n3)]
The objectives are as follows
10x1 + 15x2 + n1 - p1 = 40
100x1 + 100x2 + n2 - p2 = 1000
x2 + n3 - p3 = 7
x; n; p >= 0
The solution is x' = [4; 0] and a = [0; 600; 7]
> local({pkg <- select.list(sort(.packages(all.available =
TRUE)),graphics=TRUE)
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
> local({pkg <- select.list(sort(.packages(all.available =
TRUE)),graphics=TRUE)
+ if(nchar(pkg)) library(pkg, character.only=TRUE)})
> coeff<-matrix (c(10,15,100,100,0,1), nrow=3, ncol=2, byrow=TRUE)
> target<-c(40,1000,7)
> p1<-c(2,0,0,0,0,0)
> p2 <- c(0,0,0,0,1,0)
> p3<- c(0,0,0,0,0,1)
> achievement <- data.frame(p1,p2,p3)
> achievement
p1 p2 p3
1 2 0 0
2 0 0 0
3 0 0 0
4 0 0 0
5 0 1 0
6 0 0 1
> llgp(coeff,target,achievement)
Do you have any idea why I am seeing below error ?
Error in matrix(0, nrow = levels, ncol = nonbasics) :
invalid 'nrow' value (too large or NA)
In addition: Warning messages:
1: In max(achievements$priority) :
no non-missing arguments to max; returning -Inf
2: In matrix(0, nrow = levels, ncol = nonbasics) :
NAs introduced by coercion to integer range
Regards,
Debasis Ghosh, Ph.D
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
