Больше комментария (или серии комментариев), чем ответа, который вы ищете, но график не помещается в комментарии.
Ваш подход выглядит обреченным - если не к провалу, то к чрезвычайной трудности.
Нет никакой гарантии, что какой-либо из ваших пределов квантильного бина будет соответствовать любому из пределов бина гистограммы.
Точно так же нет гарантии, что разница между соседними квантилями будет кратна любой ширине бина гистограммы, которую вы можете выбрать.Вы могли бы поддаться искушению обмануть это, раскрасив столбец в зависимости от того, какая квантильная корзина встречается чаще, но это будет игнорировать детали.Итак, предположим, что ваша гистограмма была для [100, 200), но некоторые значения в этом интервале принадлежат одному квантильному элементу, а другие - другому: что бы вы сделали?И что бы вы сделали, если бы 3 или более квантильных бина попали в гистограмму?
Указав несколько гистограмм без указания начальных значений или ширины бинов, вы развязываете анархию.Stata будет принимать отдельные решения для каждой гистограммы, основываясь частично на размерах выборки.Вот что говорит ваш код, но не то, что вы хотите.
Ваши гистограммы ничего не знают об используемых вами аналитических весах.
Помимо этого, ваш вопрос поднимает всевозможные ненужные головоломки.
Зачем производить aux и ничего с этим не делать?Это стандартное искусство SO, чтобы показать минимальный код, необходимый для объяснения вашей проблемы.
Вы говорите, что заинтересованы в децилях, но непоследовательно также работаете с 1 и 95% процентилями.
Почему у вас есть такие нерегулярные значения с очень разными весами, неясно, но несущественно для вашего непосредственного вопроса.Но все это заставляет меня думать, что вы не можете легко или эффективно получить гистограмму, такую как пример, из ваших данных.У вас всего 53 точки данных, и поэтому вес не имеет значения, если вы не можете иметь более 53 непустых корзин.
Как границы ячеек падают относительно данных, может быть показано непосредственно без гистограммы.
С вашими примерами (спасибо!) Я делаю это
xtile aux = dn3001 [aw=hw0010], nq(10)
quantile dn3001, ms(none) mla(aux) mlabpos(0) scheme(s1color) rlopts(lc(none))
Я бы использовал логарифмическую шкалу обычно, но отрицательноценности исключают это.
Здесь я не ограничиваюсь строгими проблемами программирования, но этот вопрос неизбежно поднимает проблемы, к которым я обращаюсь.