Plotar intervalo de confiança
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar? pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n') Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"
Fátima, O post que compartilho não tem quase nada a ver com o que você pergunta não fosse o fato de que nele você pode ver uma bela demonstração do uso do latticeExtra::segplot() Acho que é um ótimo ponto de partida para você fazer os gráficos com os intervalos de confiança! https://ridiculas.wordpress.com/tag/segplot/ Todo o crédito ao Walmes! att, FH 2013/11/26 Fátima Lima Paula <fatima.lima.paula@gmail.com>
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar?
pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n')
Obrigada Fátima
-- "Minha felicidade depende da qualidade dos meus pensamentos"
_______________________________________________ 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.
Desculpe Fernando. Mas isso seria partir do mais complicado para o simplérrimo. Preciso de algo simples. Obrigada Fátima Em 26 de novembro de 2013 13:43, FHRB Toledo <fernandohtoledo@gmail.com>escreveu:
Fátima,
O post que compartilho não tem quase nada a ver com o que você pergunta não fosse o fato de que nele você pode ver uma bela demonstração do uso do latticeExtra::segplot()
Acho que é um ótimo ponto de partida para você fazer os gráficos com os intervalos de confiança!
https://ridiculas.wordpress.com/tag/segplot/
Todo o crédito ao Walmes!
att, FH
2013/11/26 Fátima Lima Paula <fatima.lima.paula@gmail.com>
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar?
pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n')
Obrigada Fátima
-- "Minha felicidade depende da qualidade dos meus pensamentos"
_______________________________________________ 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.
_______________________________________________ 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.
-- "Minha felicidade depende da qualidade dos meus pensamentos"
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 <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Valeu, Éder. Estou fazendo isso. Mas o erro ficou tão pequeno que não dá nem para ver. Vou checar e ver se está correto. Obrigada. Em 26 de novembro de 2013 16:25, Éder Comunello <comunello.eder@gmail.com>escreveu:
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 <c <comunello.eder@gmail.com>omunello.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.
-- "Minha felicidade depende da qualidade dos meus pensamentos"
Obrigada pela atenção, Éder. O erro foi meu mesmo. Já está tudo certo. Abs Fátima Em 26 de novembro de 2013 21:26, Éder Comunello <comunello.eder@gmail.com>escreveu:
Fátima, boa noite!
Para fins de apresentação, talvez você possa alterar o intervalo de confiança utilizado pra produzir um intervalo maior e ficar visível no gráfico.
Éder Comunello <c <comunello.eder@gmail.com>omunello.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.
-- "Minha felicidade depende da qualidade dos meus pensamentos"
Fátima, Veja a função plotCI do pacote gplots. HTH -- Cesar Rabak 2013/11/26 Fátima Lima Paula <fatima.lima.paula@gmail.com>
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar?
pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n')
Obrigada Fátima
-- "Minha felicidade depende da qualidade dos meus pensamentos"
_______________________________________________ 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.
Boa noite Membros, Retomando a questão de plotar o intervalo de confiança com plotCI, gostaria de fazer como no exemplo da função: # plot means and data(state) tmp <- split(state.area, state.region) means <- sapply(tmp, mean) stdev <- sqrt(sapply(tmp, var)) n <- sapply(tmp,length) ciw <- qt(0.975, n) * stdev / sqrt(n) # plain plotCI(x=means, uiw=ciw) Só que eu gostaria de ligar a parte superior e inferior do intervalo de confiança das quatro médias com uma linha contínua, dando o aspecto de envelope e ainda se possível colorir a área interna, para tanto, pergunto alguem conhece alguma função em algum pacote que faça isto? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== Em 27/11/2013 11:25, Cesar Rabak escreveu:
Fátima,
Veja a função plotCI do pacote gplots.m envelo
HTH
-- Cesar Rabak
2013/11/26 Fátima Lima Paula <fatima.lima.paula@gmail.com <mailto:fatima.lima.paula@gmail.com>>
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar?
pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n')
Obrigada Fátima
-- "Minha felicidade depende da qualidade dos meus pensamentos"
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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.
_______________________________________________ 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.
Obrigada, Cesar. Vou ver. Em 27 de novembro de 2013 20:15, ASANTOS <alexandresantosbr@yahoo.com.br>escreveu:
Boa noite Membros,
Retomando a questão de plotar o intervalo de confiança com plotCI, gostaria de fazer como no exemplo da função:
# plot means and data(state) tmp <- split(state.area, state.region) means <- sapply(tmp, mean) stdev <- sqrt(sapply(tmp, var)) n <- sapply(tmp,length) ciw <- qt(0.975, n) * stdev / sqrt(n)
# plain plotCI(x=means, uiw=ciw)
Só que eu gostaria de ligar a parte superior e inferior do intervalo de confiança das quatro médias com uma linha contínua, dando o aspecto de envelope e ainda se possível colorir a área interna, para tanto, pergunto alguem conhece alguma função em algum pacote que faça isto?
Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
Em 27/11/2013 11:25, Cesar Rabak escreveu:
Fátima,
Veja a função plotCI do pacote gplots.m envelo
HTH
-- Cesar Rabak
2013/11/26 Fátima Lima Paula <fatima.lima.paula@gmail.com>
Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar?
pop=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",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n')
Obrigada Fátima
-- "Minha felicidade depende da qualidade dos meus pensamentos"
_______________________________________________ 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.
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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.
_______________________________________________ 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.
-- "Minha felicidade depende da qualidade dos meus pensamentos"
Alexandre, bom dia! Adotando a mesma ideia utilizada no caso da Fátima, dá pra fazer algo bem simples, mas que pode servir. ### <BEGIN> require(gplots) # plot means and data(state) ### pra apagar o que não vai usar: rm(list=ls(patt='^state')[!(ls() %in% c('state.area', 'state.region'))]) tmp <- split(state.area, state.region) means <- sapply(tmp, mean) stdev <- sqrt(sapply(tmp, var)) n <- sapply(tmp,length) ciw <- qt(0.975, n) * stdev / sqrt(n) # plain plotCI(x=means, uiw=ciw) ### rbind(means,stdev,n,ciw) ### sumário ### envelope simples plotCI(x=means, uiw=ciw) eixoX <- 1:4 lines(eixoX, means+ciw, lty='dashed', col=2) lines(eixoX, means-ciw, lty='dashed', col=2) ### envelope simples hachurado plotCI(x=means, uiw=ciw) polX <- c(eixoX,rev(eixoX)) polY <- c(means+ciw,rev(means-ciw)) #colors()[grep("light",colors())] ### pesquisando cores #polygon(polX, polY, lty='dashed', col='lightgray', bor=NA) ### sem borda ao redor do polígono polygon(polX, polY, lty='dashed', col='lightgray', bor=2) ### com borda vermelha plotCI(x=means, uiw=ciw, add=T) ### <END> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Um porém... Acredito que o envelope é inapropiado no caso do exemplo dado, pois assume um comportamento contínuo entre os "tratamentos" enquanto eles não tem necessariamente alguma relação. -- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
Obrigado Éder, Ficou muito bom o intervalo de confiança colorido, faz um tempão que tento fazer isto. No meu caso são distâncias para um mesmo tratamento, usei o exemplo do gplot apenas para colocar o problema, mesmo não havendo esta continuidade, Redobrados agradecimentos, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== Em 28/11/2013 11:45, Éder Comunello escreveu:
Um porém... Acredito que o envelope é inapropiado no caso do exemplo dado, pois assume um comportamento contínuo entre os "tratamentos" enquanto eles não tem necessariamente alguma relação.
-- Éder Comunello <c <mailto:comunello.eder@gmail.com>omunello.eder@gmail.com <mailto:omunello.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.
Estou aqui fazendo um exercício pessoal. require(gplots) plotmeans(salario~escolaridade,p=0.95,col="red",connect =F,n.label=F,ci.label=F,mean.labels=F, xlab="Escolaridade",ylab="Média") Existe a possibilidade ao invés de representar as média por apenas um ponto representar as média por uma barra e no topo desta barra o colocar o IC? André Oliveira Souza Em Quinta-feira, 28 de Novembro de 2013 16:50, ASANTOS <alexandresantosbr@yahoo.com.br> escreveu: Obrigado Éder, Ficou muito bom o intervalo de confiança colorido, faz um tempão que tento fazer isto. No meu caso são distâncias para um mesmo tratamento, usei o exemplo do gplot apenas para colocar o problema, mesmo não havendo esta continuidade, Redobrados agradecimentos, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== Em 28/11/2013 11:45, Éder Comunello escreveu: Um porém... Acredito que o envelope é inapropiado no caso do exemplo dado, pois assume um comportamento contínuo entre os "tratamentos" enquanto eles não tem necessariamente alguma relação.
--
É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. _______________________________________________ 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.
Andre, bom dia! Tem um exemplo no help da {gplots} que acredito ser o que você procura. S ### <code r> require(gplots) # Example with confidence intervals and grid hh <- t(VADeaths)[, 5:1] mybarcol <- "gray20" ci.l <- hh * 0.85 ci.u <- hh * 1.15 mp <- barplot2(hh, beside = TRUE, col = c("lightblue", "mistyrose", "lightcyan", "lavender"), legend = colnames(VADeaths), ylim = c(0, 100), main = "Death Rates in Virginia", font.main = 4, sub = "Faked 95 percent error bars", col.sub = mybarcol, cex.names = 1.5, plot.ci = TRUE, ci.l = ci.l, ci.u = ci.u, plot.grid = TRUE) mtext(side = 1, at = colMeans(mp), line = -2, text = paste("Mean", formatC(colMeans(hh))), col = "red") box() ### </code> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 30 de novembro de 2013 19:57, Andre Oliveira <andreolsouza@yahoo.com.br>escreveu:
Estou aqui fazendo um exercício pessoal.
require(gplots) plotmeans(salario~escolaridade,p=0.95,col="red",connect =F,n.label=F,ci.label=F,mean.labels=F, xlab="Escolaridade",ylab="Média")
Existe a possibilidade ao invés de representar as média por apenas um ponto representar as média por uma barra e no topo desta barra o colocar o IC?
André Oliveira Souza
Em Quinta-feira, 28 de Novembro de 2013 16:50, ASANTOS < alexandresantosbr@yahoo.com.br> escreveu: Obrigado Éder,
Ficou muito bom o intervalo de confiança colorido, faz um tempão que tento fazer isto. No meu caso são distâncias para um mesmo tratamento, usei o exemplo do gplot apenas para colocar o problema, mesmo não havendo esta continuidade,
Redobrados agradecimentos,
Alexandre
-- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
Em 28/11/2013 11:45, Éder Comunello escreveu:
Um porém... Acredito que o envelope é inapropiado no caso do exemplo dado, pois assume um comportamento contínuo entre os "tratamentos" enquanto eles não tem necessariamente alguma relação.
-- Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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.
_______________________________________________ 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.
_______________________________________________ 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.
participantes (6)
-
Andre Oliveira -
ASANTOS -
Cesar Rabak -
FHRB Toledo -
Fátima Lima Paula -
Éder Comunello