Tentei achar seu email mas nao consegui, então vai por aqui mesmo. Se puder me ajudar, eu agradeço.
no seu post http://lmonasterio.blogspot.com/2011/04/entao-voce-acha-que-e-classe-media.html, vc tabula alguns percentis da renda domiciliar per capita, mas não consegui reproduzir. Voce deu algum tratamento aos dados (pesos, talvez)?
Eu apenas gerei as variáveis e calculei os percentis.
como vc pode ver, a diferença maior foi no percentil 99 (no seu tinha dado uns 4500, aqui já deu 1E+12). Posso ter feito alguma m, mas não identifiquei.
Ah... nao me lembro... teria que voltar ao codigo no spss. (sim, eu fiz no spss.... mea culpa) ... Mas, rapaz, "1e+12" é muita coisa! Vc nao deve ter limpado os dados, neh? Eh 1 quackilhao de reais...
Leonardo,
ResponderExcluirTentei achar seu email mas nao consegui, então vai por aqui mesmo. Se puder me ajudar, eu agradeço.
no seu post http://lmonasterio.blogspot.com/2011/04/entao-voce-acha-que-e-classe-media.html, vc tabula alguns percentis da renda domiciliar per capita, mas não consegui reproduzir. Voce deu algum tratamento aos dados (pesos, talvez)?
Eu apenas gerei as variáveis e calculei os percentis.
Usei os pesos amostrais da pnad, mas não a complexidade da amostra. Que versão da pnad vc usou?
ResponderExcluirDeu
A propósito, deu muito diferente?
Ps . Para me contactar, use o Google profiles.
Leonardo,
ResponderExcluirNão deu mto diferente não, vou colar aqui os códigos do R (usei o pacote survey).
...depois de ter lido os dados com aquele pacote do IBGE:
pnad <- svydesign(id = ~1, weights = ~pesos, data = dadospes)
svyquantile(~renddompcap, pnad, c(0,.05,.10,.25,.50,.75,.95,.99,.995,.999), ci=TRUE, na.rm=T)
E o resultado:
calculei os percentis 0 - 5 - 10 - 25 - 50 - 75 - 95 - 99
renda dom p capita (para cada percentil):
0 68 108 200 391 740 3487.771 1e+12
como vc pode ver, a diferença maior foi no percentil 99 (no seu tinha dado uns 4500, aqui já deu 1E+12). Posso ter feito alguma m, mas não identifiquei.
Ah, imagino o que pode ser: você trabalhou com os dados de domicílios, eu to usando o de pessoas. é isso?
ResponderExcluirAh... nao me lembro... teria que voltar ao codigo no spss. (sim, eu fiz no spss.... mea culpa)
ResponderExcluir...
Mas, rapaz,
"1e+12" é muita coisa! Vc nao deve ter limpado os dados, neh? Eh 1 quackilhao de reais...
é, tinha uma observação absurda ali. Vou deixar aqui os códigos caso te interesse. Obrigado pelos comentários!
ResponderExcluirdadospes <- dadospes[dadospes$renddompcap < 1e+7,]
pnad <- svydesign(id = ~1, weights = ~pesos, data = dadospes)
svyquantile(~renddompcap, pnad, c(0,.05,.10,.25,.50,.75,.95,.99,.995,.999), ci=TRUE, na.rm=T)
resultado (um pouco diferente):
$quantiles
0 0.05 0.1 0.25 0.5 0.75 0.95 0.99 0995 0999
renddompcap 0 98 144 257 465 820 2337.918 5300 7000 12600
Vamos produzir tabletes enquanto perdemos a indústria textil e siderúrgica. Que ótima troca não?...
ResponderExcluir