[R-br] RES.: Função runs.test()!!!

[Thiago De paula protásio]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110409/09498246/attachment.html>