inserir intervalo de confiança em grafico: interaction.plot(z$dia, z$trat, z$A)

Boa tarde, tenho uma tabela com dados em série temporal e fiz gráfico com dados de Amônia do controle e tratamento com réplicas. Preciso plota rno gráfico o intervalo de confiança referente as réplicas do tratamento E e do tratamento C para cada dia. Alguém pode me auxilar de como faço isso.
trat dia replica A
E 1 1 0,31
E 1 2 0,31
E 1 3 0,31
E 1 4 0,31
C 1 1 0,31
C 1 2 0,31
C 1 3 0,31
C 1 4 0,31
E 7 1 0,01
E 7 2 0,02
E 7 3 0,04
E 7 4 0,03
C 7 1 0,52
C 7 2 0,63
C 7 3 0,58
C 7 4 0,64
E 14 1 0,02
E 14 2 0,05
E 14 3 0,06
E 14 4 0,04
C 14 1 0,62
C 14 2 0,63
C 14 3 0,68
C 14 4 0,64
E 21 1 2,70
E 21 2 0,40
E 21 3 0,00
E 21 4 0,04
C 21 1 35,00
C 21 2 23,10
C 21 3 33,60
C 21 4 30,80

Olá, Roberta.

Você pode fazer isso com o pacote sciplot. O gráfico dos seus dados fica ligeiramente estranho por não existir variação no dia 1, mas isso é da natureza dos dados, e não um problema em si.

Veja aqui como fica:

Código:

dados<-structure(list(trat = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L), .Label = c("C", "E"), class = "factor"),
    dia = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 7L, 7L, 7L, 7L, 7L,
    7L, 7L, 7L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 21L,
    21L, 21L, 21L, 21L, 21L, 21L, 21L), replica = c(1L, 2L, 3L,
    4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L,
    3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L),
    A = c(0.31, 0.31, 0.31, 0.31, 0.31, 0.31, 0.31, 0.31, 0.01,
    0.02, 0.04, 0.03, 0.52, 0.63, 0.58, 0.64, 0.02, 0.05, 0.06,
    0.04, 0.62, 0.63, 0.68, 0.64, 2.7, 0.4, 0, 0.04, 35, 23.1,
    33.6, 30.8)), class = "data.frame", row.names = c(NA, -32L
))

summary(dados)


library(sciplot)

lineplot.CI(dados$dia, dados$A, dados$trat, xlab="dia", ylab="A")

Basta copiar e colar no R para ver, ok? Qualquer coisa vamos nos falando.

Abraços

Grata. Deu certo. Tem como mudar a posição da legenda C e E? E deixar as letras com outra cor?

Boa noite,
A dúvida que ficou é como posso fazer o script do sciplot baixando os dados de uma tabela em .txt
Fiz novamente e deu certo.

Segue script e o os dados do arquivo teste.txt que usei para tentar fazer:

setwd("~/R")
teste<-read.table("teste.txt",header=T);teste

summary(teste)
               
library(sciplot)
               
lineplot.CI(teste$dia, teste$A, teste$trat, xlab="dia", ylab="A")

  dia trat replica     A
1    1    E      1L  0.01
2    1    E      2L  0.01
3    1    E      3L  0.01
4    1    E      4L  0.01
5    1    C      1L  0.01
6    1    C      2L  0.01
7    1    C      3L  0.01
8    1    C      4L  0.01
9    7    E      1L  0.15
10   7    E      2L  0.25
11   7    E      3L  0.18
12   7    E      4L  0.08
13   7    C      1L  0.10
14   7    C      2L  0.29
15   7    C      3L  0.42
16   7    C      4L  0.33
17  14    E      1L  0.01
18  14    E      2L  0.02
19  14    E      3L  0.00
20  14    E      4L  0.00
21  14    C      1L  0.48
22  14    C      2L  0.34
23  14    C      3L  0.28
24  14    C      4L  0.17
25  21    E      1L  0.09
26  21    E      2L  0.06
27  21    E      3L  0.01
28  21    E      4L  0.07
29  21    C      1L 15.17
30  21    C      2L 16.22
31  21    C      3L 14.07
32  21    C      4L 17.07

Oi, Roberta!

Então o problema está acontecendo quando você tenta importar os dados, é isso?

Se sim, copie e cole aqui a primeira mensagem de erro que aparece, e com qual comando ela surge, para descobrirmos o que está errado, ok?

Abraços!

Deu certo , consegui.