One way that you can obtain domain statistics including percentiles like medians from the PROC SURVEYxxx procedures is to create new survey weights from the original survey weights, and run these procedures using these new survey weights. For observations within a specified level of a domain, the new survey weights equal the original survey weights, but for observations at another level of that domain, the new survey weights equal a small positive number close to zero (for example, 0.000000001). For each level of a domain, create similar new survey weights. Then, run the PROC SURVEYxxx procedures across all the levels of a domain using the corresponding new survey weights. If you do this often, you can "macro-ize" the creation of these new survey weights and merge them to the original data set. The justification for this method is found in the following reference: Graubard BI, Korn EL. Survey inference for subpopulations. American Journal of Epidemiology 1996;144:102-106.
... View more