I have a variable "education" which is 3-level and ordinal and I have a binary variable "urban" which equals to '1' if the individual is in urban area or '0' if they are not. I also have sample weights in a variable "sampleWeights" to scale my data up to a full county level-these weight values are provided by the study investigators.

I wish to do a chisq2 test to see if the education distribution is different for urban vs non-urban and do this at the county level using the weights.

Now to do this without using 'sampleWeights' I can do this:
Code:
 
 tabulate education urban, chi2 
However I am wondering, how can I incorporate these weights which are 'pweights' ?