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)
Rep Dose Tempo Leitura
A 0 10 5.2
B 0 10 4.5
C 0 10 4.5
A 0 21 4.5
B 0 21 4.2
C 0 21 4.2
A 0 30 4.4
B 0 30 4
C 0 30 3.9
A 0 37 4.11
B 0 37 3.93
C 0 37 3.91
A 0 50 4.14
B 0 50 3.9
C 0 50 3.88
A 0 62 4.07
B 0 62 3.84
C 0 62 3.82
A 0 70 4.44
B 0 70 4.17
C 0 70 4.05
A 0 82 4.23
B 0 82 3.99
C 0 82 3.97
A 0 98 4.06
B 0 98 3.92
C 0 98 4.03
A 0 108 4.05
B 0 108 3.94
C 0 108 3.9
A 0 118 4.04
B 0 118 4.3
C 0 118 4.19
A 0 129 4.42
B 0 129 4
C 0 129 3.9
A 0 144 4.19
B 0 144 4.08
C 0 144 4.08
A 4 10 5.8
B 4 10 5.8
C 4 10 6
A 4 21 5.6
B 4 21 5.6
C 4 21 5.7
A 4 30 5.7
B 4 30 5.8
C 4 30 5.8
A 4 37 5.92
B 4 37 5.85
C 4 37 5.96
A 4 50 5.97
B 4 50 5.97
C 4 50 6
A 4 62 5.78
B 4 62 5.98
C 4 62 5.96
A 4 70 5.91
B 4 70 6.01
C 4 70 6.15
A 4 82 5.53
B 4 82 5.66
C 4 82 5.83
A 4 98 5.8
B 4 98 5.78
C 4 98 6
A 4 108 5.84
B 4 108 5.81
C 4 108 5.9
A 4 118 6.39
B 4 118 6.17
C 4 118 6.1
A 4 129 5.7
B 4 129 5.85
C 4 129 5.95
A 4 144 5.95
B 4 144 5.98
C 4 144 6.02
A 8 10 5.9
B 8 10 6.2
C 8 10 6.3
A 8 21 5.9
B 8 21 5.9
C 8 21 6.1
A 8 30 6.2
B 8 30 6.3
C 8 30 6.2
A 8 37 6.41
B 8 37 6.39
C 8 37 6.44
A 8 50 6.53
B 8 50 6.57
C 8 50 6.5
A 8 62 6.53
B 8 62 6.63
C 8 62 6.54
A 8 70 6.29
B 8 70 6.59
C 8 70 6.38
A 8 82 6.15
B 8 82 6.35
C 8 82 6.4
A 8 98 6.32
B 8 98 6.46
C 8 98 6.49
A 8 108 6.6
B 8 108 6.67
C 8 108 6.69
A 8 118 6.41
B 8 118 6.36
C 8 118 6.5
A 8 129 6.79
B 8 129 6.8
C 8 129 6.68
A 8 144 6.42
B 8 144 6.67
C 8 144 6.69
A 12 10 6.1
B 12 10 6.2
C 12 10 6.1
A 12 21 6.1
B 12 21 6.2
C 12 21 6.3
A 12 30 6.5
B 12 30 6.6
C 12 30 6.6
A 12 37 6.57
B 12 37 6.58
C 12 37 6.63
A 12 50 6.71
B 12 50 6.7
C 12 50 6.71
A 12 62 6.81
B 12 62 6.84
C 12 62 6.9
A 12 70 6.87
B 12 70 6.76
C 12 70 6.79
A 12 82 6.57
B 12 82 6.55
C 12 82 6.63
A 12 98 6.84
B 12 98 6.89
C 12 98 6.92
A 12 108 6.79
B 12 108 6.9
C 12 108 6.97
A 12 118 6.48
B 12 118 6.48
C 12 118 6.54
A 12 129 6.8
B 12 129 6.78
C 12 129 6.81
A 12 144 6.91
B 12 144 6.71
C 12 144 6.9
A 16 10 6.5
B 16 10 6.6
C 16 10 6.6
A 16 21 6.5
B 16 21 6.6
C 16 21 6.2
A 16 30 6.6
B 16 30 6.7
C 16 30 6.7
A 16 37 6.78
B 16 37 6.8
C 16 37 6.79
A 16 50 6.8
B 16 50 6.79
C 16 50 6.84
A 16 62 6.87
B 16 62 6.92
C 16 62 6.88
A 16 70 6.79
B 16 70 6.89
C 16 70 6.95
A 16 82 6.8
B 16 82 6.61
C 16 82 6.77
A 16 98 7.08
B 16 98 6.98
C 16 98 6.95
A 16 108 7.12
B 16 108 7.04
C 16 108 7.08
A 16 118 6.82
B 16 118 6.82
C 16 118 6.76
A 16 129 7.02
B 16 129 7.19
C 16 129 7.11
A 16 144 6.92
B 16 144 7.11
C 16 144 7.16
A 20 10 6.3
B 20 10 6.2
C 20 10 6.2
A 20 21 6.4
B 20 21 6.2
C 20 21 6.2
A 20 30 6.7
B 20 30 6.7
C 20 30 6.7
A 20 37 6.79
B 20 37 6.79
C 20 37 6.85
A 20 50 6.82
B 20 50 6.85
C 20 50 6.89
A 20 62 6.87
B 20 62 6.92
C 20 62 7.01
A 20 70 7.03
B 20 70 6.97
C 20 70 7.02
A 20 82 6.71
B 20 82 6.78
C 20 82 6.83
A 20 98 7.09
B 20 98 7
C 20 98 7.01
A 20 108 7.06
B 20 108 7.12
C 20 108 7.05
A 20 118 6.83
B 20 118 6.97
C 20 118 7.09
A 20 129 6.96
B 20 129 7.04
C 20 129 7.06
A 20 144 6.91
B 20 144 7.12
C 20 144 7.16
A 24 10 6.3
B 24 10 6.2
C 24 10 6.5
A 24 21 6.4
B 24 21 6.2
C 24 21 6.4
A 24 30 6.6
B 24 30 6.7
C 24 30 6.7
A 24 37 6.86
B 24 37 6.85
C 24 37 6.89
A 24 50 6.93
B 24 50 6.88
C 24 50 6.93
A 24 62 6.93
B 24 62 6.96
C 24 62 7.02
A 24 70 7.03
B 24 70 6.94
C 24 70 7
A 24 82 6.73
B 24 82 6.78
C 24 82 6.82
A 24 98 7.04
B 24 98 7.05
C 24 98 7.04
A 24 108 7.15
B 24 108 7.21
C 24 108 7.23
A 24 118 6.82
B 24 118 6.99
C 24 118 7.04
A 24 129 7.09
B 24 129 7.15
C 24 129 7.16
A 24 144 7.2
B 24 144 7.17
C 24 144 7.19
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
