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!
_______________________________________________
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.
