[R-br] ajuda com implementação de Programação Linear com o R - constrOptim
Cleber N.Borges
klebyn em yahoo.com.br
Terça Maio 24 23:18:25 BRT 2016
Prezado Cesar
Boa noite,
no caso a variável x também é uma igualdade a uma constante
(na função foi passado como valor padrão - default )
seria possível usar otimização não linear?
por acaso já viu algo parecido com o problema?
Agradeço por sua atenção
abraço
Cleber
Em 24/05/2016 23:05, Cesar Rabak escreveu:
> Cleber,
>
> Olhando de "soslaio" as suas condições de contorno vê-se que não se
> trata de um problema de LP. . .
>
> Prá complicar suas condições de contorno têm um grau de liberdade (a
> variável x é livre[?])...
>
>
> 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 ) )
>
---
Este email foi escaneado pelo Avast antivírus.
https://www.avast.com/antivirus
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160524/cea0e5d9/attachment.html>
Mais detalhes sobre a lista de discussão R-br