Re: [R-br] Erro nls: matriz gradiente singular com estimativas de parâmetros iniciais
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 deve ser este!! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}
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 <ivanalaman@yahoo.com.br>escreveu:
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 deve ser este!!
(S,f,P) Allaman * * \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}
_______________________________________________ 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.
Alguém tem ideia de como resolver este problema? Será que o problema está no modelo? Em 8 de agosto de 2012 14:25, Gustavo Dias Azevedo <gustavoazevedo@id.uff.br
escreveu:
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 < ivanalaman@yahoo.com.br> escreveu:
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 deve ser este!!
(S,f,P) Allaman * * \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}
_______________________________________________ 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.
CMR....
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 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 ==========================================================================
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 software não consegue fazer a minha regressão, mostrando a mensagem de erro **matriz gradiente singular com estimativas de parâmetros iniciais.* Em 9 de agosto de 2012 10:37, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
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 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.
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 <gustavoazevedo@id.uff.br>:
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 software não consegue fazer a minha regressão, mostrando a mensagem de erro matriz gradiente singular com estimativas de parâmetros iniciais.
Em 9 de agosto de 2012 10:37, Walmes Zeviani <walmeszeviani@gmail.com> escreveu:
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 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.
_______________________________________________ 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.
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 seu CMR. Acredito que os seus dados não sejam muitos, menos de 50 linhas talvez. Sendo assim você usar dput() para gerar um CMR, assim incorporando os dados no corpo da mensagem. Veja este exemplo http://r-br.2285057.n4.nabble.com/R-br-homocedasticidade-modelo-nao-linear-t... À 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 ==========================================================================
participantes (4)
-
Benilton Carvalho -
Gustavo Dias Azevedo -
Ivan Bezerra Allaman -
Walmes Zeviani