<font face="trebuchet ms,sans-serif">Thiago,<br><br>Segue abaixo um CMR no estilo do que deveria ser (pelo menos para mim) as trocas de ajuda dentro da lista. Esclarecimentos estão comentados dentro do código.<br><br><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># exemplo repruzível (exemplo para os demais inscritos da lista)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># dados hopesdados pelo site <a href="http://www.datafilehost.com/">http://www.datafilehost.com/</a></span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># link para a página dos dados <a href="http://www.datafilehost.com/download-a7e94384.html">http://www.datafilehost.com/download-a7e94384.html</a></span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># link do download <a href="http://www.datafilehost.com/get.php?file=a7e94384">http://www.datafilehost.com/get.php?file=a7e94384</a></span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># o melhor é fazer com que os dados sejam baixados na leitura</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># esse procedimento não funcionou para os seus dados. porque?</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">da <- read.table("<a href="http://www.datafilehost.com/get.php?file=688e5891">http://www.datafilehost.com/get.php?file=688e5891</a>", header=TRUE) # falha</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># hospedei os dados no site do LEG, assim eu posso baixar direto da web</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># temos que arrumar algo assim para que todos possam usar</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">da <- read.table("<a href="http://www.leg.ufpr.br/~walmes/docs/umi.txt">http://www.leg.ufpr.br/~walmes/docs/umi.txt</a>", header=TRUE) # falha</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">str(da)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># diagnose gráfica</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">plot(UR~UEQ, da)</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># pelo jeito, você fixou a UR para obter os valores de UEQ. então não seria UEQ em função</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># de UR? vou continuar a análise para mostrar os procedimentos, mas pense sobre isso.</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># modelo logístico não vai se ajustar aos seus dados porque eles não tem informação</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"># para ajustar uma curva sigmoide. os dados tem padrão de crescimento assintótico</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># por isso ajustar um modelo que alçance essa comportamento</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">A <- 100   # mudar até obter bons chutes</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">B <- 0.11  # mudar até obter bons chutes</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">plot(UR~UEQ, da)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">curve(A*(1-exp(-B*x)), add=TRUE)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># ajuste do modelo com os chutes</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">n0 <- nls(UR~A*(1-exp(-B*UEQ)), data=da,</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">          start=list(A=100, B=0.11))</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">summary(n0)</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># curva ajustada</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">A <- coef(n0)["A"]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">B <- coef(n0)["B"]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">plot(UR~UEQ, da)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">curve(A*(1-exp(-B*x)), add=TRUE, col=2)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"># Código mínimo reproduzível: tudo se resolve em questão de segundos!</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">#------------------------------------------------------------------------------------------</span><br style="font-family: courier new,monospace;">
<br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Walmes Marques Zeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">fone: (+55) 41 3361 3573</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">VoIP: (3361 3600) 1053 1173</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">twitter: @walmeszeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">linux user number: 531218</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br>

<br>