y<-read.csv("y.csv",header=T)
y<-as.matrix(y)
b<-read.csv("b.csv",header=T)
b<-as.matrix(b)
x1<-read.csv("x1.csv",header=T)
x1<-as.matrix(x1)
x2<-read.csv("x2.csv",header=T)
x2<-as.matrix(x2)
x3<-read.csv("x3.csv",header=T)
x3<-as.matrix(x3)
library(lpSolve)
re<-matrix(rep(0,14*27),nrow=14,ncol=27)
for(t in 1:14){
for(k in 1:27){
f.obj<-c(rep(0,27),1)
f.con<-matrix(rbind(c(y[t,],-y[t,k]),c(b[t,],b[t,k]),c(x1[t,],x1[t,k]),c(x2[t,],x2[t,k]),c(x3[t,],x3[t,k])))
f.dir<-c(">=","=","<=","<=","<=")
f.rhs<-c(y[t,k],b[t,k],x1[t,k],x2[t,k],x3[t,k])
re[t,k]<-lp("max",f.obj,f.con,f.dir,f.rhs)$solution[28]
}
}
re
请问循环结构有啥问题,R语言一运行就崩溃了,是有死循环吗?


雷达卡




京公网安备 11010802022788号







