A idéia de registrar pensamentos, experiências, fatos, num diário é antiga. Mas fazer esse diário de acesso público é uma coisa que só o milagre da internet permitiu. É um desafio interessante. Um diálogo com alguém, qualquer um, todos. Mais provavelmente, ninguém ou muito poucos. Quem é que tem algo interessante a dizer?
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!
Assinar:
Postar comentários (Atom)
Um comentário:
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?
Postar um comentário