Prezados,<br><br>Estou tentando ajustar dados a um modelo de cinética química do tipo<br><br>y = K*x(^a)*w(^b)/(z^c)<br><br>x, w e z são as concentrações dos meus reagentes e produtos (dados)<br>K, a, b e c são os parâmetros que eu preciso determinar no modelo.<br>
<br>Utilizei o seguinte comando<br><br>dados <- read.table("Teste2.txt",header=TRUE)<br>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))<br><br>Porém, ao executar a última linha, recebo a seguinte mensagem:<br>
<br><b>matriz gradiente singular com estimativas de parâmetros iniciais</b><br><br>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.<br>
<br>Alguém sabe como me ajudar? Devo trocar de comando ou fazer alguma manipulação no modelo?<br><br>Obrigado!<br>