<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
<font face="Comic Sans MS">Caros colegas,<br>
<br>
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).<br>
<br>
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).<br>
<br>
Agradeço se puderem me ajudar<br>
<br>
## Experimento pH vs tempo<br>
<br>
dados <- read.table("dosagem.txt", h=T)<br>
names(dados)<br>
attach(dados)<br>
<br>
B <- factor(Rep)<br>
D <- factor(Dose)<br>
T <- factor(Tempo)<br>
<br>
<br>
modelo <- lm(Leitura ~ D*T+B:D)<br>
anova(modelo)<br>
<br>
# Pode-se considerar o modelo fatorial diretamente<br>
<br>
fatorial <- lm(Leitura ~ D*(I(Tempo)+I(Tempo^2)))<br>
anova(fatorial)<br>
summary(fatorial)<br>
#plot(fatorial)<br>
<br>
#SELECIONANDO O MODELO QUADRATICO E GERANDO INTERVALOS DE
CONFIANÇA DE 0.95<br>
<br>
preditos <- predict(fatorial,interval="confidence")<br>
<br>
media <- preditos[,1]<br>
LI <- preditos[,2]<br>
LS <- preditos[,3]<br>
<br>
interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias", ylim=c(4,8),
pch=0:7, legend=T)<br>
<br>
#Plotando os pontos e as curvas ajustadas - Grafico 1<br>
par(new=TRUE)<br>
interaction.plot(Tempo,D,Leitura,ylab="pH",xlab="Dias",
type="p",pch=0:7,col=3:9, ylim=c(4,7.5),legend=T)<br>
par(new=TRUE)<br>
interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias",
ylim=c(4,7.5), pch=0:7, legend=F)<br>
<br>
#Plotando a curva ajustada e os IC<br>
<br>
par(new=TRUE)<br>
interaction.plot(Tempo,D,ylab="pH", media,ylim=c(4,7.5),legend=F)<br>
par(new=TRUE)<br>
interaction.plot(Tempo,D,ylab="pH",xlab="Dias",LI,col=2,ylim=c(4,7.5),legend=F)<br>
par(new=TRUE)<br>
interaction.plot(Tempo,D,ylab="pH",xlab="Dias",
LS,col=3,ylim=c(4,7.5),legend=F)<br>
legend("topleft", c("IC - Superior", "IC - Inferior"),
col=c("green", "red"),<br>
title="Legenda", lty=1, lwd=2)<br>
<br>
<br>
</font>
</body>
</html>