data MFH;
   input STR CLU CHRON SYSBP X;
   wgt=1;
   label
      STR="Stratum ID"
      CLU="Cluster ID"
      CHRON="Cluster sample sum of Chronic morbidity (CHRON)"
      SYSBP="Cluster sample sum for Systolic blood pressure (SYSBP)"
      X="Number of sample elements in the clusters";
datalines;
   1  1 70 29056  204 
   1  2 74 29417  210
   2  1 12  3692   26
   2  2 14  4564   30
   3  1 15  7741   59
   3  2 16  8585   63
   4  1  9  6277   45
   4  2 14  5668   43
   5  1 10  2322   17
   5  2 16  3960   30
   6  1 10  3080   21
   6  2  6  3252   22
   7  1 10  3966   27
   7  2  4  3261   24
   8  1 12  4156   28
   8  2  6  2852   20
   9  1 15  6617   46
   9  2 23  6616   48
   10 1 37 10552   73
   10 2 25 11032   77
   11 1 11  8759   60
   11 2 25  9876   72
   12 1 33  9901   69
   12 2 24  6828   47
   13 1 31  8624   61
   13 2 27  9390   66
   14 1 22  6960   48
   14 2 20  7130   49
   15 1 18  6646   49
   15 2 22  7094   49
   16 1 24  9841   69
   16 2 37 11786   83
   17 1 19  6910   48
   17 2 23  6446   45
   18 1 25 10742   73
   18 2 29  9026   61
   19 1 36  9350   65
   19 2 34  8912   62
   20 1  9  3810   26
   20 2 22  7098   51
   21 1 18  6998   53
   21 2 34  9970   69
   22 1 29 11146   79
   22 2 41 13215   94
   23 1 22  6596   48
   23 2 18  6002   41
   24 1 15  3808   27
   24 2  7  3148   22
;
run;


data MFH;
   merge mfh survey.weights;
run;


proc surveymeans data=mfh varmethod=brr(fay=.29289322);
   var chron sysbp x;
   ratio chron sysbp/ x;
   repweights w1-w1000;
   weight wgt;
run;