dear All  when the 2 codes are run I don't find the same predictive probability that is supposed to be used as propensity score. Any idea? Here is the code: data School; Music= 'Yes'; do j=1 to 60; if (ranuni(1312) > 0.4) then Gender='Female'; else Gender='Male'; Absence = ranexp(99); if (Absence < 0.5) then GPA= 4 + rannor(99)/3.5; else GPA= 4 - abs(rannor(99)/3.5); if (Gender='Female') then GPA= GPA + 0.02; id1= ranuni(99); id2= ranuni(99); output; end; Music= 'No'; do j=1 to 100; Absence= ranexp(99); if (ranuni(99) > 0.45) then Gender='Female'; else Gender='Male'; if (Absence < 0.5) then GPA= 4 + rannor(99)/4.2; else GPA= 4 - abs(rannor(99)/4.2); id1= ranuni(99); id2= ranuni(99); output; end; do j=1 to 40; Absence= 2 + ranexp(99); if (ranuni(99) > 0.5) then Gender='Female'; else Gender='Male'; GPA= 3.4 - abs(rannor(99)/4.2); id1= ranuni(99); id2= ranuni(99); output; end; run; proc sort data=school; by id1; run; data school; set school; StudentID= _n_; Absence= int(Absence*100+0.5) / 100; GPA= int(GPA*100+0.5) / 100; run; proc sort data=school; by id2; run; data Grades; set school; keep StudentID GPA; run; data School; set school; drop j id1 id2 GPA; run; proc sort data=School out=School1; by StudentID; run; proc sort data=Grades out=Grades1; by StudentID; run; data SchoolGrades; merge School1 Grades1; by StudentID; run;   ods graphics on; proc psmatch data=School region=treated; class Music Gender; psmodel Music(Treated='Yes')= Gender Absence; match distance=lps method=greedy(k=1) exact=Gender caliper=0.5; assess lps var=(Gender Absence) / stddev=pooled(allobs=no) stdbinvar=no plots(nodetails)=all weight=none; output out(obs=match)=OutEx4 matchid=_MatchID; run; ****** logistic ****; proc logistic data=School descending; class Music Gender; model Music(Event='Yes')= Gender Absence / link=cloglog rsquare; output out = ps_los pred = ps xbeta=logit_ps; run; ... View more