O bootstrapping funciona bem para acessar a incerteza na estimativa média, no entanto, lembro-me de ler em algum lugar que o bootstrap não faz um bom trabalho ao avaliar a incerteza nas estimativas quantílicas (particularmente a mediana).
Não lembro onde li isso e não consegui encontrar muita coisa com uma rápida pesquisa no Google. Considerações sobre isso e quaisquer referências serão muito apreciadas.
sqreg
comando (regressão simultânea-quantil) no Stata estima os erros padrão. Mas isso não prova nada, eu sei.