Fátima, boa tarde!

Se você tiver a estimativa de erro armazenada em um vetor, pode utilizar lines() para plotar o limiar.

### <BEGIN>

popfaF  <- c(37374,  54783,  86049, 126048, 287384)

obfaF   <- c(6,15,27,55,239)

txobfaF <- (obfaF*365/popfaF)*100; txobfaF

idade   <- c("1","2","3","4","5")

plot(idade, txobfaF, type="b", lty=2, ylim=c(1,40), xaxt='n',

     main="Sexo Feminino", xlab="Idade", ylab="Taxa de Mortalidade")

set.seed(457)

erro <- runif(5,1,5) ### substitua por seu cálculo do erro!

lines(idade, txobfaF+erro, lty='dashed', col=2)

lines(idade, txobfaF-erro, lty='dashed', col=2)

### <END>

Éder Comunello <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.