I have the following dataset below. I have been trying to subset the data 1) have patients (Patient_No) with any prior positive outcome and any subsequent negative outcome, 2) have patients with at least 1 prior positive outcome, at least 1 subsequent negative outcome, at least 90 days between the first low positive or negative outcome and most recent negative outcome, most recent outcome should be negative and there are no subsequent positive outcomes. I have tried using Proc sql first and last to select the patients but the code is getting too long and messy.
Data Diagnose;
Input @1 Patient_No $2.
@3 Date MMDDYY10.
@14 Visit_No $2.
@16 Outcome $12.;
Format Date MMDDYY10.;
Datalines;
1 10/21/2000 1 Positive
1 10/25/2000 2 Positive
1 11/01/2000 3 Negative
1 05/28/2001 4 Negative
2 11/22/2000 1 Positive
2 11/29/2000 2 Positive
2 12/28/2000 3 Positive
2 06/28/2001 4 Low positive
2 10/29/2001 5 Negative
3 12/12/2000 1 Positive
3 12/29/2000 2 Positive
3 02/21/2001 3 Positive
3 07/12/2001 4 Negative
3 08/29/2001 5 Positive
;
I hope this implements all of your rules:
data sub1 sub2;
do until(last.patient_no);
set diagnose; by patient_no;
select (Outcome);
when ("Positive") do;
if missing(firstPos) then firstPos = date;
end;
when ("Negative") do;
if missing(firstNeg) then firstNeg = date;
recentNeg = date;
end;
when("Low positive") do;
if missing(firstLow) then firstLow = date;
end;
otherwise;
end;
end;
/* Subset 1 : any prior positive outcome and any subsequent
negative outcome */
if firstPos < firstNeg
then output sub1;
/* Subset 2 */
if
/* at least 1 prior positive outcome */
not missing(firstPos) and
/* at least 1 subsequent negative outcome */
firstPos < firstNeg and
/* at least 90 days between the first low positive
or negative outcome and most recent negative outcome */
intck("day", min(firstLow, firstNeg), recentNeg) >= 90 and
/* most recent outcome should be negative and there are no
subsequent positive outcomes */
Outcome = "Negative"
then output sub2;
keep patient_No; /* Applies to both output datasets */
run;
Please post what you've tried, otherwise we may make suggestions that will not work for you.
I hope this implements all of your rules:
data sub1 sub2;
do until(last.patient_no);
set diagnose; by patient_no;
select (Outcome);
when ("Positive") do;
if missing(firstPos) then firstPos = date;
end;
when ("Negative") do;
if missing(firstNeg) then firstNeg = date;
recentNeg = date;
end;
when("Low positive") do;
if missing(firstLow) then firstLow = date;
end;
otherwise;
end;
end;
/* Subset 1 : any prior positive outcome and any subsequent
negative outcome */
if firstPos < firstNeg
then output sub1;
/* Subset 2 */
if
/* at least 1 prior positive outcome */
not missing(firstPos) and
/* at least 1 subsequent negative outcome */
firstPos < firstNeg and
/* at least 90 days between the first low positive
or negative outcome and most recent negative outcome */
intck("day", min(firstLow, firstNeg), recentNeg) >= 90 and
/* most recent outcome should be negative and there are no
subsequent positive outcomes */
Outcome = "Negative"
then output sub2;
keep patient_No; /* Applies to both output datasets */
run;
@PGStats, the code works fine, but I noticed that if the first visit is a negative, the code fails to meet the rules. For example, in the data set below. Patient 4 meets the rules, but has a negative first visit.
Data Diagnose;
Input @1 Patient_No $2.
@3 Date MMDDYY10.
@14 Visit_No $2.
@16 Outcome $12.;
Format Date MMDDYY10.;
Datalines;
1 10/21/2000 1 Positive
1 10/25/2000 2 Positive
1 11/01/2000 3 Negative
1 05/28/2001 4 Negative
2 11/22/2000 1 Positive
2 11/29/2000 2 Positive
2 12/28/2000 3 Positive
2 06/28/2001 4 Low positive
2 10/29/2001 5 Negative
3 12/12/2000 1 Positive
3 12/29/2000 2 Positive
3 02/21/2001 3 Positive
3 07/12/2001 4 Negative
3 08/29/2001 5 Positive
4 10/21/2000 1 Negative
4 11/25/2001 2 Positive
4 12/01/2001 3 Positive
4 06/28/2002 4 Negative
4 10/26/2002 5 Negative
;
April 27 – 30 | Gaylord Texan | Grapevine, Texas
Walk in ready to learn. Walk out ready to deliver. This is the data and AI conference you can't afford to miss.
Register now and lock in 2025 pricing—just $495!
Still thinking about your presentation idea? The submission deadline has been extended to Friday, Nov. 14, at 11:59 p.m. ET.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.
Ready to level-up your skills? Choose your own adventure.