proc sql; create table sumproduct as select -4.608753929 as value1 ,0.28179071 as value2 ,col1,col2,col3 , calculated value1 + (calculated value2 *col1) as col4 , EXP( calculated col4)/(1+EXP(calculated col4))as col5 , sum(calculated col5*col3)/sum(col3) as sum1, sum(col2*col3)/sum(col3) as sum2 from have; quit; /*need a do-loop for value1?*/ data sump; set sumproduct; if sum1=sum2 then value1=-4.56266639; run; ... View more

%let value1 = -4.608753929; %let value2 = 0.28179071; data want; retain sumproduct1 sumproduct2 sum_col3 0; set have; col4 =&value1. + (&value2.*col1); col5=EXP(col4)/(1+EXP(col4)); sumproduct1=sumproduct1+col5*col3; col5_col3=col5*col3; sumproduct2=sumproduct2+col2*col3; col2_col3=col2*col3; sum_col3=sum_col3+col3; sumproduct1_=sumproduct1/sum_col3; sumproduct2_=sumproduct2/sum_col3; xx=col4-col5; run; sumproduct(col5, col3)/sum(col3) = sumproduct(col2,col3)/sum(col3)?? ... View more

/*reupload the code*/ data test1; set enrolled; if drugA=1 then do; if (end_dt + 7) >= fst_dt > fill_dt then a_mis=0; *no drugA misuse; if (end_dt + 7) < fst_dt then a_mis=1; *drugA misuse, because we expect it to be out of the system by then; end; if drugB=1 then do; if (end_dt + 30) >= fst_dt > fill_dt then b_mis=0; * no drugB misuse; if (end_dt + 30) < fst_dt then b_mis=1; *misuse of drugB; end; run; /*get unique mis-use record per patid per fst_dt*/ data a; set test1; where a_mis=1; keep patid fst_dt a_mis; proc sort nodupkey; by patid fst_dt ; run; data b; set test1; where b_mis=1; keep patid fst_dt b_mis; proc sort nodupkey; by patid fst_dt ; run; ... View more

Sorry, I don't totally understand your question, so I guess that you want to get unique mis-use records per patid per fst_dt? if not, could you mind to explain it clearly, filter condition or what you want? /*get unique mis-use records per patid per fst_dt*/ data a; set test1; where a_mis=1; keep patid fst_dt a_mis; proc sort nodupkey; by patid fst_dt ; run; data b; set test1; where b_mis=1; keep patid fst_dt b_mis; proc sort nodupkey; by patid fst_dt ; run; ... View more

data test1; set have; array a admin_date1-admin_date6; array b admin_date1-admin_date6 ; array x $ type1-type6; array y $ type1-type6; array v(5,6) $10 a1b1-a1b6 a2b1-a2b6 a3b1-a3b6 a4b1-a4b6 a5b1-a5b6; do i=1 to 5; do j=i+1 to 6; /* case 1: If both types = P, the difference between the dates must be 17 or more. case 2: All other combinations must be 24 or more days apart*/ if (x(i)='P' and y(j)='P' and b(j)-a(i)>=17) or ((x(i)^='P' or y(j)^='P') and b(j)-a(i)>=24) then do; v(i,j)='Y'; output; end; end; end; run; data want; set test1; format admin_date_1-admin_date_6 mmddyy10.; array a admin_date1-admin_date6; array b admin_date_1-admin_date_6 ; array x $ type1-type6; array y $ type_1-type_6; array v(5,6) $10 a1b1-a1b6 a2b1-a2b6 a3b1-a3b6 a4b1-a4b6 a5b1-a5b6; do n=1 to 5; do m=n+1 to 6; if v(n,m)='Y' then do; b(n)=a(n); y(m)=x(m); end; end; end; keep person i j admin_date_1-admin_date_6 type_1-type_6; run; ... View more

https://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_phreg_sect016.htm ... View more

ods output ParameterEstimates= Est; proc phreg data=have; model AVAL*y(1)=sexn AGEGR1N/ALPHA=0.05 rl; run;   hi , try the option -- ALPHA=0.05 rl ... View more