[R-br] ajuda com implementação de Programação Linear com o R - constrOptim

Elias Teixeira Krainski eliaskrainski em yahoo.com.br
Quinta Maio 26 04:17:14 BRT 2016


Isso é um problema de Programação Linear (matemática aplicada, não 
estatística). Veja este (free) livro (em R):
http://www.r-bloggers.com/modeling-and-solving-linear-programming-with-r-free-book/

Elias

> 2016-05-14 9:27 GMT-03:00 Cleber N.Borges <klebyn em yahoo.com.br 
> <mailto:klebyn em yahoo.com.br>>:
>
>     Prezados,
>     Bom dia!
>     Como na lista tem pessoas com conhecimento e experiência, peço
>     algumas orientações... :-)
>
>     Imagino que o método de  "Programação Linear (LP)"
>     deve servir para resolver meu problema de achar os parâmetros de
>     uma dada função objetivo.
>
>     Sei que no R tem a opção nativa "constrOptim" (poderia ser outra,
>     posso instalar pacotes)
>
>     Entretanto eu não sei como elaborar minhas condições de contorno
>     em formato matricial
>     que é o modo como o comando constrOptim trabalha (até onde
>     consegui entender)
>
>     Agradeceria muito se alguém puder me ajudar,
>     respondendo se eu realmente estou certo em pensar que a LP é a
>     ferramenta certa para meu caso
>     e /(ou que) me indique um material básico (introdutório).
>
>     Abaixo segue um exemplo simplificado do tipos de função que
>     pretendo obter parâmetros.
>
>     Agradeço antecipadamente pela ajuda e por qualquer dica
>
>     Cleber Borges
>
>     #####################
>     #
>     fobjetivo <- function( x = 0.1 ){
>     k1 <- 0.007585776
>     k2 <- 6.16595e-08
>     k3 <- 2.137962e-13
>     k4 <- 1e-14
>     k <- k1 * k2 * k3 * k4
>     # Determinar as 6 variáveis: a, b, c, d, e, f
>     ###################
>     # 6 Condições de contorno: ( Igualdades )
>     ###################
>     k1 == a * d / c            # 1
>     k2 == a * e / d            # 2
>     k3 == a * f / e            # 3
>     k4 == a * b                # 4
>     4 * a == b + d + e + f     # 5
>     x     == c + d + e + f     # 6
>     return( c( a, b, c, d, e, f ) )
>

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160526/4beb5578/attachment.html>


Mais detalhes sobre a lista de discussão R-br