Thanks! ... View more

Thanks.  ... View more

Thank you Reeza!   I guess my question is not just for PROC PSMATCH, but about "weighted mean" in general, including PROC MEANS.   Do we calculate the weighted mean of X when Weight = Y by calculating the mean of (X1*Y1, ... Xn*Yn)?    Thanks! ... View more

Thanks! Please take a look when you have a chance. ... View more

Thank you Reeza for your help. Here are my codes:    proc psmatch data=drugs region=allobs; class Drug Gender; psmodel Drug(Treated='Drug_X')= Gender Age BMI; psweight weight=atewgt nlargestwgt=6; assess lps var=(Gender Age BMI) / varinfo plots=(barchart boxplot(display=(lps BMI)) wgtcloud); id BMI; output out(obs=all)=OutEx1 weight=_ATEWgt_; run; data weight; set OutEx1; age_wt=age*_ATEWgt_; bmi_wt=bmi*_ATEWgt_; ps_wt=_ps_*_ATEWgt_; run; proc sort; by drug; proc means noprint; by drug; var age age_wt bmi bmi_wt _ps_ ps_wt _atewgt_; output out=mean mean=agem age_wtm bmim bmi_wtm _ps_m ps_wtm _atewgtm sum=ages age_wts bmis bmi_wts _ps_s ps_wts _atewgts; run; proc print;run; The PROC PSMATCH part (and parts prior that) was copied from Example 98.1.  None of the weighted means (*_wtm) match that in the tables 98.1.2 and 98.1.3, even though _atewgts match the weights in Tables 98.1.2 and 98.1.3.    Thanks! ... View more

I calculated "means" for age, bmi, _ps_ and the "weighted means" based on the products  _ATEWgt_ *(age, bmi, and _ps_) using PROC MEANS.    All the "means" match the means in the tables but none of my "weighted means", even though the sums of _ATEWgt_ match the ones in the tables for each drug.    Apparently, _ATEWgt_ was not supposed to be used as the multiplier to get the weighted means. Then, how to get the weighted means?    Thanks! ... View more