segunda-feira, 23 de março de 2009

Dica do Stata: gráficos com intervalo de confiança



O Stata faz gráficos super bonitos, mas nem sempre é fácil chegar onde a gente quer. Em epidemiologia, a gente toda hora quer fazer gráficos de pontos ou de barras com intervalos de confiança. Aqui vão comandos básicos para gerar gráficos desse tipo.

Este gráfico é gerado pelo comando

twoway rcap icinf icsup escol || scatter coef escol, legend(off) xscale(range(.7 4.3))
xlabel(, valuelabels) ytitle(Mean development score)


Em itálico estão os nomes das variáveis: icinf e icsup são os limites de cada intervalo de confiança.

Este outro gráfico é gerado pelo comando

twoway (bar p10pct escore, barwidth(.5)) (rcap p10pctsup p10pctinf escore),
legend(order(1 "% in P10" 2 "95% CI" )) ylabel(0 (10) 60) ytitle(Percent)

Idem acima! Brincando com o comando e danda uma olhada na ajuda do Stata vc vai conseguir entender o que cada instrução faz! Boa diversão!

Um comentário:

Unknown disse...

No help do Stata tem um exemplo que ele calcula os valores máximos e mínimos pela fórmula: by month: egen lo = min(volume). Como foi calculado o icinf e o icsup, por favor?