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 <[email protected] <mailto:[email protected]>>:

    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 ) )


_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Responder a