BookmarkSubscribeRSS Feed
asterina
Calcite | Level 5

Hi all,

I wonder if you can please help me on something that has got me stuck over the whole of last week!

I have visit information for a large number of patients. Every patient has a varying number of visits and visit date. The task is to classify the current medication treatment they have been prescribed by their doctor by looking ahead over the next 6 months and considering what else they may have got prescribed in the future (of next 180 days). For simplicity of describing this example, think of it whether we classify this as a concurrent medication or not - depending if they have prescribed anything else in the future also.

So I am really stuck how to do this - there is a large number of patients, a large number of records for each - some patients visited nearly every day, others once every 8 months. The information goes up to a few years for each. So it's tricky to set any hardcoded rules for them, or indeed to just look at next few observations.

Here is an example of what the data looks like. So for this example, we want to know if Patient1 was prescribed anything else over next 180 days other than drug x, which would then allow us to create a new variable indicating that "drug x" was given as a concurrent medication. ny example codes on how to solve this would be massively appreciated.

Visit_dateMax_dateTrt_at_visit
Patient101/12/0030/05/01Drug_x
Patient131/12/0029/06/01Drug_x
Patient124/03/0120/09/01Drug_e
Patient122/07/0118/01/02Drug_x
Patient102/12/0031/05/01Drug_x
Patient101/01/0130/06/01Drug_x
Patient125/03/0121/09/01Drug_x
Patient123/07/0119/01/02Drug_x
Patient120/11/0119/05/02Drug_a
Patient120/12/0118/06/02Drug_x
Patient113/03/0209/09/02Drug_x
Patient111/07/0207/01/03Drug_x
Patient108/11/0207/05/03Drug_x
Patient108/12/0206/06/03Drug_x
Patient101/03/0328/08/03Drug_x
Patient229/06/0326/12/03Drug_x
Patient227/10/0324/04/04Drug_e
Patient226/11/0324/05/04Drug_x
Patient217/02/0415/08/04Drug_x
Patient216/06/0413/12/04Drug_x
Patient214/10/0412/04/05Drug_e
Patient213/11/0412/05/05Drug_x
Patient204/02/0503/08/05Drug_x
Patient204/06/0501/12/05Drug_x
Patient202/10/0531/03/06Drug_x
Patient201/11/0530/04/06Drug_x
Patient223/01/0622/07/06Drug_x
Patient223/05/0619/11/06Drug_x
4 REPLIES 4
ballardw
Super User

Can you provide an example of the desired output for Patient1 based on the example data?

What role if any does the Max_date variable have?

asterina
Calcite | Level 5

thanks! Max_date is just the visit date window (visit date +180 days). I would like a new variable which will recode the treatment variable based on its current value - e.g. Single medication or conmed

DBailey
Lapis Lazuli | Level 10

proc sql;

create table want as

select

     t1.patient,

     t1.visit_date,

     t1.Trt_at_visit,

     sum(case when t1.trt_at_visit = t2.trt_at_visit then 1 else 0 end) as Future_Rx_Current_Med,

     sum(case when t1.trt_at_visit <> t2.trt_at_visit then 1 else 0 end) as Future_Rx_Different_Med

from

     have t1

     left outer join have t2

          on t1.patient = t2.patient

          and t2.visit_date > t1.visit_date

          and t1.visit_date < t1.visit_date + 180

group by

     t1.patient,

     t1.visit_date,

     t1.trt_at_visit;

quit;

Haikuo
Onyx | Level 15

Was not aware of you have done it yesterday ago, was still waiting on your new post. anyway, here is my suggested solution:

data have;

     infile cards dlm='09'x truncover;

     input patient:$10.   (Visit_date     Max_date) (:ddmmyy10.)     Trt_at_visit$;

     format visit_date Max_date ddmmyy10.;

     cards;

Patient1   01/12/00   30/05/01   Drug_x

Patient1   31/12/00   29/06/01   Drug_x

Patient1   24/03/01   20/09/01   Drug_e

Patient1   22/07/01   18/01/02   Drug_x

Patient1   02/12/00   31/05/01   Drug_x

Patient1   01/01/01   30/06/01   Drug_x

Patient1   25/03/01   21/09/01   Drug_x

Patient1   23/07/01   19/01/02   Drug_x

Patient1   20/11/01   19/05/02   Drug_a

Patient1   20/12/01   18/06/02   Drug_x

Patient1   13/03/02   09/09/02   Drug_x

Patient1   11/07/02   07/01/03   Drug_x

Patient1   08/11/02   07/05/03   Drug_x

Patient1   08/12/02   06/06/03   Drug_x

Patient1   01/03/03   28/08/03   Drug_x

Patient2   29/06/03   26/12/03   Drug_x

Patient2   27/10/03   24/04/04   Drug_e

Patient2   26/11/03   24/05/04   Drug_x

Patient2   17/02/04   15/08/04   Drug_x

Patient2   16/06/04   13/12/04   Drug_x

Patient2   14/10/04   12/04/05   Drug_e

Patient2   13/11/04   12/05/05   Drug_x

Patient2   04/02/05   03/08/05   Drug_x

Patient2   04/06/05   01/12/05   Drug_x

Patient2   02/10/05   31/03/06   Drug_x

Patient2   01/11/05   30/04/06   Drug_x

Patient2   23/01/06   22/07/06   Drug_x

Patient2   23/05/06   19/11/06   Drug_x

;

proc sql;

     create table want as

           select *, exists(select * from have where   patient=a.patient

                and a.visit_date <= visit_date <= a.visit_date+180

                and Trt_at_visit ne a.Trt_at_visit

                ) as Concurrent_flag

     from have a

     ;

quit;

Good luck,

Haikuo

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 4 replies
  • 725 views
  • 0 likes
  • 4 in conversation