set.seed(45644529)
t_collected <- c()
t_notcollected <- c()
WRT <- 8
CT<- 0## tiempo acumulado de trabajo
WH <- 8
Trec <- 0
Tret_land <- 0
TrecACT <- 0
TTB<-array(
rnorm(20,5,.5),
dim = c(20, 20, 1),
dimnames = list(paste0("In", 1:20), paste0("Out", 1:20), paste0("Truck"))
)
#(5,max(20,rnorm(5,5,10)))
TC <- 10
ACT <-replicate(20,10)
t_landfill <-0
t_returning <-0
t_recollection <- 0
# repeat{if (CT[i]>= WH){
# break
# }
for (i in 1:ncol(TTB)){ #primer set the condiciones
if(is.na(TTB[i])<= is.na(ACT[i])){
ACT[i]<-ACT[i]-TTB[i]
t_recollection[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
Trec[i]=t_recollection[i]+CT[i]
if(is.na(Trec[i]) <= WH){
CT[i]<- sum(t_recollection,t_landfill,t_returning)
Trec[i]=t_recollection[i]+CT[i]
t_collected[i]<-TTB[i]
t_notcollected[i] <-ACT[i]
TrecACT[i] <- (ACT[i]*t_recollection[i])/t_collected[i]
}else if (is.na(TrecACT[i])<= is.na(WH)){
t_recollection[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
t_collected[i]<-ACT[i]
TTB[i] <- TTB[i]-ACT[i]
t_notcollected[i] <-TTB[i]
ACT[i]<-0
CT[i]<- sum(t_recollection,t_landfill,t_returning)
Tret_land[i] <-t_returning[i]+t_landfill[i]+CT[i]
}else if(Tret_land[i]<=WH){
t_recollection[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
t_landfill[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
t_returning[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
t_collected[i]<-ACT[i]
CT[i]<- sum(t_recollection,t_landfill,t_returning)
ACT[i]=TC
}else{
t_recollection[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
t_landfill[i]<-min(0.5,max(1.5,rnorm(1,1,0.5)))
ACT[i]=ACT[i]-(ACT[i]*(WRT[i]/t_recollection[i]))
TTB[i]=TTB[i]-(ACT[i]*(WRT[i]/t_recollection[i]))
CT[i]=8
}
}
t_recollection[i]=t_recollection[i]
t_landfill[i]=t_landfill[i]
t_returning[i]=t_returning[i]
CT[i] <- CT[i]
WRT[i] <- WH-CT[i]
Trec[i] <- Trec[i]
TrecACT[i] <- TrecACT[i]
Tret_land[i] <- Tret_land[i]
t_notcollected[i] <-t_notcollected[i]
TTB[i] <- TTB[i]
ACT[i] <- ACT[i]
ACT[i+1]<-ACT[i]
WRT[i+1] <- WRT[i]
t_recollection[i+1]=t_recollection[i]
t_landfill[i+1]=t_landfill[i]
t_returning[i+1]=t_returning[i]
CT[i+1]=CT[i]
Trec[i+1] <- t_recollection[i]+CT[i]
t_notcollected[i+1] <-t_notcollected[i]
TrecACT[i+1] <- (ACT[i]*t_recollection[i])/t_collected[i]
Tret_land[i+1] <-t_returning[i]+t_landfill[i]+CT[i]
TTB[i+1] <- TTB[i]
if(CT[i]==WH){
break
}
}
Result <-cbind(TTB=TTB[1:20],t_collected=t_collected[1:20],#AcumTcollec,
CT=CT[1:10],t_recollection=t_recollection[1:20],
t_landfill=t_landfill[1:20],
t_returning=t_returning[1:20],
ACT=ACT[1:20],Trec=Trec[1:20],WRT=WRT[1:20],
t_notcollected=t_notcollected[1:20],TrecACT=TrecACT[1:20],
Tret_land= Tret_land[1:20] )
This topic was automatically closed 21 days after the last reply. New replies are no longer allowed.
If you have a query related to it or one of the replies, start a new topic and refer back with a link.