## Ler os dados dados<-read.table("teste1.txt",h=T) ## define a variável X como fator para fazer a análise de variância x1<-as.factor(dados$x) ## define a variável X como numérica para o modelo de regressão x<-dados$x ## variável resposta y<-dados$y ## Modelo de regressão mod<-lm(y~x) anova(mod) summary(mod) ## Modelo de análise de variância mod2<-lm(y~x1) anova(mod2) summary(mod2) ## Faz o desdobramento de tratamento em regressão e desvio de regressão mod3<-lm(y~x+x1) anova(mod3) #summary(mod3) ## teste para falta de ajuste do desvio de regressão ## modelo de regressão e modelo de análise de variância anova(mod,mod2) ## Cria um vetor de médias para a variável resposta para cada nível de X med<-tapply(y,x1,mean);med ## Cria uma sequencia dos níveis de X xp<-c(0,30,40,50,60,70,80) # Plota os valores da variável resposta plot(dados) # Adiciona a reta de regressão ajustada no modelo de regressão abline(mod,lwd=2) # Adiciona a sequencia de valores X e médias da variável resposta lines(xp,med,type="p",col="red",pch=19)