CMR
___
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.
Não temos o dádiva de adivinhar as coisas a partir do vácuo. Por isso envie
um CMR.
À disposição.
Walmes.
==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de
Segue o comando utilizado. Os dados encontram-se hospedados no DataFileHost
e o link encontra-se imbutido no comando.
dados - read.table(http://www.datafilehost.com/download-c367edc8.html
,header=TRUE)
modelo - nls(y~K*(x^a)*(w^b)*(z^c),data=dados,start=list(K=4,a=3,b=2,c=1))
*# a partir daqui o
seu exemplo nao e' reproduzivel como dado abaixo... mas apos baixar o
arquivo, observe que vc tem 3 variaveis preditoras: x, w e z
dai', note que
w = 11-x
z = 2x
entao, com esse nivel de colinearidade, nao vejo mesmo como ter alguma
convergencia.
b
2012/8/9 Gustavo Dias Azevedo
Gustavo,
O seu exemplo não é reproduzível!
dados - read.table(http://www.datafilehost.com/download-c367edc8.html;,
header=TRUE)
Erro em scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings,
:
linha 1 não tinha 6 elementos
Antes de submeter inicie uma nova sessão do R e teste o
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 -
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
Se você indicou exatamente os valores não deveria dar erro!! Perceba que
primeiro vc nos informa o seguinte modelo:
y = K*x(^a)*w(^b)/(z^c)
e depois na função você informa outro modelo:
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))
Provavelmente o erro
Foi apenas um erro meu de digitação. O modelo correto encontra-se como
escrito na função: (y~K*(x^a)*(w^b)*(z^c).
Em 8 de agosto de 2012 12:18, Ivan Bezerra Allaman
ivanala...@yahoo.com.brescreveu:
Se você indicou exatamente os valores não deveria dar erro!! Perceba que
primeiro vc nos informa