[R-br] Coeficientes dos modelos polinomiais ajustados

[Wenceslau]

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)


-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110504/5715e806/attachment.html>
-------------- Próxima Parte ----------
Um texto embutido e sem conjunto de caracteres especificado foi limpo...
Nome: dosagem.txt
Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110504/5715e806/attachment.txt>