<font face="trebuchet ms,sans-serif">Thiago,<br><br>Você tem que debruçar-se mais sobre os seus dados/modelos. Pega o comando e vai testando modelo, inventa formas funcionais, não fique preso a modelos de páginas de livros. Eu ajustei outro modelo e observei um melhor ajuste. Ainda assim há falta de ajuste com relação ao modelo de médias, veja com relação ao modelo de médias, que diga-se, usa todos os graus de liberdade possíveis e não permite extrapolação. O modelo de médias é sem dúvida o melhor modelo em termos de ajuste, vai dá melhor R2 que qualquer outro, mas não dá para extrapolar. Já o não linear permite isso.<br>
<br>Causas da falta de ajuste:<br>* o modelo não ser mesmo o correto para representar, ser uma versão aproximada de um modelo com mais parâmetros que traduzam aquela forma dos dados;<br>* haver uma causa de variação associada a cada UR que tá deslocando todas as repetições de uma UR, coisa que acontece com frequência, por exemplo se você medisse todas as repetições da UR_i em um dia, da UR_j no dia seguinte, estaria aí o efeito do dia, ou da pessoa que faz as determinações no laboratório, ou do local de onde vieram, ou do aparelho usado, das condições do momento.<br>
<br>O que fazer:<br>* se o foco são as estimativas dos parâmetros, usar o modelo como está;<br>* se o foco é apenas fazer predição, usar uma função suave da UR, lowess, spline, etc, ou uma linear segmentada;<br>* tentar uma função com mais parâmetros para alcançar mais o comportamento dos dados;<br>
<br>Em todos os casos, você vai ter que sentar e espremer possibilidades sobre isso, fazer o seu trabalho.<br><br><span style="font-family: courier new,monospace;">#-----------------------------------------------------------------------------</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)</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;">plot(UEQ~UR, da)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">with(da, tapply(UEQ, UR, length))</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">A <- 5</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">B <- 0.01   # mudar até obter bons chutes</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">C <- 0.085  # mudar até obter bons chutes</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">plot(UEQ~UR, da)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">curve(A+B*exp(C*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;">n0 <- nls(UEQ~A+B*exp(C*UR), data=da,</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">          start=list(A=A, B=B, C=C))</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;">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;">C <- coef(n0)["C"]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">plot(UEQ~UR, da)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">curve(A+B*exp(C*x), add=TRUE)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">abline(v=c(40,50,60), col=3) # causas da falta de ajuste</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">m0 <- lm(UEQ~factor(UR), da)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">anova(n0, m0)</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">qqnorm(residuals(n0))</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">plot(residuals(n0)~fitted(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;">
<br>À disposição.<br>Walmes.<br><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>