Your code indicates that you want to use the Weights variable, but your example shows expected values that aren't affected by the Weights variable.  If your example is correct then you might be able to use something like: data have;   infile cards dlm='09'x;   input SAMPLEID CLUSTER $ WEIGHTS ITEM61 ITEM5 ITEM69 ITEM27 ITEM48         ITEM707 ITEM314 ITEM11 ITEM412 ITEM36 ITEM124 ITEM1415;   cards; 800023 no re 345 0 0 0 85 0 76 5 42 0 67 0 15 800037 no re 28 0 0 23 0 0 59 0 34 0 85 0 46 800454 no re 380 0 0 20 67 0 0 56 84 0 16 0 73 8004702 no re 478 0 0 0 51 0 0 48 28 0 46 0 49 800523 no re 367 0 0 48 0 0 0 56 34 0 58 0 57 800580 no re 899 0 0 0 645 0 0 78 0 0 67 0 113 800614 no re 257 0 0 146 0 0 34 56 0 0 101 0 345 ; proc summary data=have nway;   var item61--item1415;   output out=need1 (drop=_:) min=; run; data need1;   format type $40.;   set need1;   type='Min'; run; proc summary data=have nway;   var item61--item1415;   output out=need2 (drop=_:) max=; run; data need2;   format type $40.;   set need2;   type='Max'; run; data test;   set have;   array items(*) item61--item1415;   do _n_=1 to dim(items);     if items(_n_) eq 0 then call missing(items(_n_));   end; run; proc summary data=test nway;   var item61--item1415;   output out=need3 (drop=_:) n=; run; data need3;   format type $40.;   set need3;   type='Number of times the code assigned'; run; proc summary data=test nway;   var item61--item1415;   output out=need4 (drop=_:) mean=; run; data need4;   format type $40.;   set need4;   array items(*) item61--item1415;   do _n_=1 to dim(items);     if missing(items(_n_)) then items(_n_)=0;   end;   type='Mean'; run; data want;   set need3 need4 need1 need2; run; ... View more

Hello Reeza, Thank you for the prompt response. Number of respondents is around 11,000 and variables - around 2,000... Do you think this could be an issue..? The paper that I mentioned in my post says to use ACECLUS before the actual PROC CLUSTER.. So I was confused, whether I can do this with my data. ... View more