Caros colegas,

estou tentando fazer uma analise comparando modelos lineares da evolução do pH vs tempo, num experimento com diferentes doses de carbonato de cálcio (7 doses).

Tenho este Tenho este script que faz a analise e também os gráficos. Mas não consegui obter os coeficientes de cada modelo individualmente (dose 0 vs tempo; dose 4 vs tempo, etc).

Agradeço se puderem me ajudar

## Experimento pH vs tempo

dados <- read.table("dosagem.txt", h=T)
names(dados)
attach(dados)

B <- factor(Rep)
D <- factor(Dose)
T <- factor(Tempo)


modelo <- lm(Leitura ~ D*T+B:D)
anova(modelo)

# Pode-se considerar o modelo fatorial diretamente

fatorial <- lm(Leitura ~ D*(I(Tempo)+I(Tempo^2)))
anova(fatorial)
summary(fatorial)
#plot(fatorial)

#SELECIONANDO O MODELO QUADRATICO E GERANDO INTERVALOS DE CONFIANÇA DE 0.95

preditos <- predict(fatorial,interval="confidence")

media <- preditos[,1]
LI <- preditos[,2]
LS <- preditos[,3]

interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias", ylim=c(4,8), pch=0:7, legend=T)

#Plotando os pontos e as curvas ajustadas - Grafico 1
par(new=TRUE)
interaction.plot(Tempo,D,Leitura,ylab="pH",xlab="Dias", type="p",pch=0:7,col=3:9, ylim=c(4,7.5),legend=T)
par(new=TRUE)
interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias", ylim=c(4,7.5), pch=0:7, legend=F)

#Plotando a curva ajustada e os IC

par(new=TRUE)
interaction.plot(Tempo,D,ylab="pH", media,ylim=c(4,7.5),legend=F)
par(new=TRUE)
interaction.plot(Tempo,D,ylab="pH",xlab="Dias",LI,col=2,ylim=c(4,7.5),legend=F)
par(new=TRUE)
interaction.plot(Tempo,D,ylab="pH",xlab="Dias", LS,col=3,ylim=c(4,7.5),legend=F)
legend("topleft", c("IC - Superior", "IC - Inferior"), col=c("green", "red"),
title="Legenda", lty=1, lwd=2)