finally, the code below gave me the required imputation. thanks all! I used multiple imputation with Fully  conditional specification (FCS)     proc mi data= data minimum = 0 0 0 0 MINMAXITER=100 nimpute=10 out= imuputed_data ; var  (variables); fcs plots=trace(mean std); fcs discrim( categorical variables) /classeffects=include) nbiter =100 ; fcs reg (continuous variables) nbiter =100 ; run;       mini function was used in some to prevent getting some negative minimum values. minmaxter is to re-imupute if the system gets a negative ( or any avoided number in minimum)      ... View more