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