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_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
