Erro nls: matriz gradiente singular com estimativas de parâmetros iniciais
Prezados, Estou tentando ajustar dados a um modelo de cinética química do tipo y = K*x(^a)*w(^b)/(z^c) x, w e z são as concentrações dos meus reagentes e produtos (dados) K, a, b e c são os parâmetros que eu preciso determinar no modelo. Utilizei o seguinte comando dados <- read.table("Teste2.txt",header=TRUE) modelo <- nls(y~K*(x^a)*(w^b)*(z^c),data=dados,start=list(K=1,a=0.33,b=0.66,c=-1.33)) Porém, ao executar a última linha, recebo a seguinte mensagem: *matriz gradiente singular com estimativas de parâmetros iniciais* Já procurei alguns tópicos e algumas pessoas tiveram o mesmo problema em ajustes diferentes, e encontrei diversas "soluções" diferentes, como usar o argumento alg="plinear", trocar os chutes iniciais etc, porém nenhuma delas resolveu meu problema. Alguém sabe como me ajudar? Devo trocar de comando ou fazer alguma manipulação no modelo? Obrigado!
O problema está nos chutes, que avaliando a matriz do modelo gerada com os valores iniciais resulta em uma matriz singular, não pode inverter, dentro do passo de otimização. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Prezado Walmes, Eu também já tentei trocar os chutes iniciais, mas o problema persiste. Os chutes que apresentei no tópico são os valores exatos do modelo (os valores de y no conjunto de dados são exatamente aqueles quando eu defino K=1, a=0.33, b=0.66, c=-1.33 ) visto que criei dados exatos apenas para testar a utilidade do comando e da regressão não linear para o meu modelo de cinética química. Estou tentando pensar numa manipulação algébrica e procurando outros comandos/pacotes que possam resolver o meu problema. Se alguém souber de algo, por favor me avise. Em 8 de agosto de 2012 10:22, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
O problema está nos chutes, que avaliando a matriz do modelo gerada com os valores iniciais resulta em uma matriz singular, não pode inverter, dentro do passo de otimização.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
participantes (2)
-
Gustavo Dias Azevedo -
Walmes Zeviani