Estou há horas esperando meu código terminar de rodar. Imagino que deve haver uma forma mais inteligente de fazer isso, afinal não preciso visitar todos os pontos quando tenho uma função monótona ao redor de um ponto de mínimo.
Meu problema está relacionada a teoria das filas, tenho o tempo médio na fila, o c
e a taxa de saída.
lambda <-seq(0,ceiling(mu*PDV),0.0001)
n.lambda<-seq_along(lambda)
diferenca <-NULL
diferenca1<-300
for (i in n.lambda) {
try(wq<-queueing::QueueingModel(
queueing::NewInput.MMC(lambda=lambda[i], mu=mu, c=PDV))$Wq, silent = TRUE)
if(is.na(wq)) next
diferenca<-(wq-Wq)^2
if(diferenca<diferenca1){
lambdaWin <-lambda[i]
}
diferenca1<-diferenca
}
Alguma sugestão de como melhorá-lo?
optim
doR
? Acredito que se tu editar a pergunta com algumafunção_de_interesse
, alguém aqui vai ser capaz de te ajudar a implementar uma solução utilizandooptim
.queueing::QueueingModel( queueing::NewInput.MMC(lambda=lambda[i], mu=mu, c=PDV))$Lq
Tenho o comprimento da fila é quero saber o lambda que retornaria isso.