<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Walmes e demais membros da lista;<br><br>Com base naquelas dados que eu envie a lista ....Ajustei o modelo abaixo e gostaria de testar se o mesmo é compatível......<br><br>Também ajustei um polinômio de ordem 2 e fiz uma anova() entre os dois modelos....Vc acha que esse procedimento é válido?<br><br>Dessa forma, estava pesquisando e encontrei esse teste de sequências da função runs.test(). Dê uma olhada nos comandos abaixo.....<br><br><br>n0 <- nls(UEQ~A*(1-exp(B*UR)), data=da,<br>          start=list(A=-2, B=0.03))<br>summary(n0)<br><br># curva ajustada<br><br>A <- coef(n0)["A"]<br>B <- coef(n0)["B"]<br>plot(UEQ~UR, da)<br>curve(A*(1-exp(B*x)), add=TRUE, col=2,lty=1)<br>title("Tratamento 1 ")<br><br>SQE <- summary(n0)$sigma^2*summary(n0)$df[2]<br>SQE<br><br>SQT <-
 var(da$UR)*(length(da$UR)-1)# soma de quadrado total corrigida<br><br>R2 <- 1 - SQE/SQT<br>R2<br><br>erro1=residuals(n0)<br><br>#pacotes necessários....<br>library(tseries)<br>library(quadprog)<br>library(zoo)<br><br>runs.test(factor(sign(erro1))) # permite testar se os dados diferem da  curva ajustada   <br><br><br><br>########Modelo polinomial<br>model=lm(UEQ~UR+I(UR^2))<br>summary(model)<br><br>anova (model, n0)<br><br><br>Muito obrigado pelas dicas ..vou olhar aqui no banco de teses da UFLA e encontrar a sua dissertação.....<br>Valeu mesmo.....<br>Suas dicas tem me ajudo bastante!                           <br><br>Thiago de Paula Protásio<br>
Acadêmico de Engenharia Florestal<br>
Universidade Federal de Lavras<br>
Laboratório de Energia da Biomassa Florestal<br>
Laboratório de Biomateriais<br>
(035) 9183-2246</td></tr></table>