(r1) InitialGenMedOpt < ORUA

Tags: view all tags
<!--
      * Set ALLOWTOPICCHANGE = %USERSWEB%.TWikiAdminGroup,%USERSWEB%.AmeliaWhite,%USERSWEB%.MichaelOSullivan,%USERSWEB%.CameronWalker
-->

---+!! Optimising Registrar Rosters in General Medicine: Initial Project (Amelia White)

%TOC{depth="2"}%
---++ Problem Description

General medicine (GM) regristrars at Auckland Hospital are currently rostered to 'on call' shifts by a 12 week rotating roster. This roster has grown organically over the years and is subject to union rules. The GM registrars on call admit patients to their GM ward from the emergency department and the Acute and Planning Unit (APU). There are 4 GM wards, with 3 registrars attatched to each ward. The types of shifts that are allocated to registrars include 'long day' admitting shifts (0800-2230), and normal admitting shifts (0800-1600). If a registrar is not allocated to an admitting shift then they complete a 'Normal Day' (0800-1600).

The current roster requires development as currently there are high fluctuations in patient numbers across registrars. These flucations result in some registrars being overworked on some days while other registrars have a low number of patients. This affects the standard of care for patients and the quality of life for registrars. Also of importance to consider is 'continuation of care', that is where the same registrar is looking after a patient for the duration of the patient's stay.
---++ Objective

The purpose of this project is to investigate how the current roster for General Medicine registrar. The main goal is to implement a working model of an "Idealized roster". The timeline of this project is 10 weeks and as such is not expected to produce a full working roster but to provide the basis for further investigation, possibly as a 4th year project.

The desired outcomes of the summer project are:
   * Develop a detailed understanding of how the roster works and what is required 
   * Identify a significant subset of the important constraints to include in the optimization of the roster governing measures of patient care and staff quality-of-life. 
   * Develop a code shell for roster generation. 
   * Produce example of an improved roster 

---++ Progress

Initial meeting held 25th November with Andrew Peterson, Tim Denison, Rupert Handy, Nicolas Szecket, Mike O'Sullivan and Amelia White.

From the initial meeting it seems that the main issues is variation in patient numbers across the teams which results in registrars being overworked and can compromise continuation of care for patients. This variation means that on some days a registrar may have to see upwards of 15 patients whereas on other days they may have 2 or 3. This means that on a given day one sub team may be extremely busy whereas another team may be quiet. The important aspect at this stage to address appears to be "smoothing" out the number of patients admitted across the week for each team to reduce variability.

---+++ Current Roster Notes
 The current roster is based from historical preferences but it is not known how it was developed. There are 4 GM wards and 4 teams: White, Gold, Red and Black, each consisting of 3 sub-teams. Each team is allocated to a ward.

A and B call each take approximately 1/3 of the new patients each day, C and D call take the remaining 1/3. As of December 13th 2010 this will change so each admitting team will take 1/4 of the incoming patients each. Patients admitted by the night cover registrar are split evenly across all wards.

Patients stay on average 3.5 days, however this is a long tailed distribution with around 50% of the patients being discharged within 24 hours. 40% of the time there are not enough GM beds.

Workload of registrars tends to peak later in the day and thus their shifts are not currently allocated around the peak work demand periods.

Admissions are 10-20% lower on weekends and there is a cap of 30 maximum admissions to the GM wards.

More Post-acute shifts would be beneficial in the weekend if possible.

*Current Roster Rules*

   * Shifts A B C D call must be allocated to one registrar Mon-Friday 
   * Night shift and three days off must be allocated Sunday-Thursday 
   * A relief registrar covers the shifts of the registrar on night duty. They can only work one long day (A or B call) in any given 7 day period and do not work weekends. 
   * Friday A call must do P shift (post-acute) on Saturday and Sunday 
   * Saturday A call must do P shift on Sunday 
   * Friday C and D call's patients are covered by the 'P' call on the weekends 
   * Registrars must have at least every second weekend off 

*Shift Rules*

   * Shift length must be less than 15 hours 
   * Minimum of an 8 hour break between shifts 
   * A registrar can do a maximum of one shift per day 
   * Registrars are paid from 0730 to 1730, Monday to Friday, regardless of whether they work those hours or not 
   * A registrar must be replaced by a RMO of equal or greater standing 

---+++ Preliminary Formulation

*Objective:* Produce a schedule

*Constraints:*
   * Every week day must have shifts 'A B C D' allocated to one and only one sub team 
   * Every weekend day must have shifts 'A' and 2 'P' allocated to one and only one sub team 
   * Sunday through to Thursday must have 'N' night shift allocated to one and only one sub team. 
   * Each sub team can be allocated at most one shift 

*Addition constraints to add from meeting held 03/12/2010*

   * Friday 'A' call does a post acute ('P') shift on Saturday and Sunday. -- added 21/12/2010 
   * Saturday A call does a post acute ('P') on Sunday 
---++ Hospital Visit

To gain further understanding to exactly how things currently run at Auckland hospital, I shadowed several of the staff in General Medicine. I shadowed Dr Janar Sathanathan for several hours around the APU. Patients can be admitted directly to the APU from their GP to avoid double handing by the ED. One thing that I found surprising was the amount of time paper work took compared to the time actually needed to be spent with the patient.

Dr Anthony Jordan Showed me around the APU and Emergency Department (ED) and discussed several of the current issues. One issue brought up is that not all team positions are actually filled and there is currently no-one allocated to the relief roster. This means that that if a ward is short of a registrar, that ward has to stretch to cover the absent registrars work-load. He also mentioned that it may be useful to talk to Sarah Tracey at ARRMOS who currently develop the rosters based on templates in Excel. Apparently usually the first time the rosters come out, there is usually a problem with it not meeting some of the rules such as working two weekends in a row.

The Doctors that are on call admit from the ED and APU. In addition to this they can be allocated the Emergency pager which means they have to respond immediately to any critical cases. They can also be allocated the GP phone, which is the line GP's call in on if they are sending a patient to hospital. It appears having this phone can be highly disruptive to the doctors schedule.

I followed around a locum registrar in the ED for a couple of hours in the evening. The late afternoon to early evening is usually the busiest but the night I was there it was quite quiet. I discovered the next morning that some of the patients he admitted went to the team I was following the next morning. This meant that that no-one on the team that they were admitted to had the next day had seen the patient before morning rounds. In the ED each patient is assessed based on priority and can either be dispatched or, if they are expected to only be in for a short time, admitted to the APU or finally they can be admitted to the wards.

The morning following the evening I was in ED, I attended the GM morning meeting. This meeting occurs at 8am every weekday morning. The purpose is to review patients that have come in over-night and which wards and teams they have been allocated to. Ideally patients are allocated to the ward belonging to the team they have been assigned too. This unfortunately doesn't always happen and part of the meeting is spent 'finding' 'lost' patients. I shadowed Dr Rupert Handy, we first visited the new patients who had been kept in the APU overnight and then saw the new patients who had been allocated to the wards.
---+++ Example of patient allocation:

The below tables detail patient allocation before and after patients were moved between teams at the 8am meeting. The reason patients have to be moved between teams is when they have been sent to the wrong ward from ED or the APU.

<table width="80%" cellpadding="2" border="1"> <tbody><tr><th colspan="5">Before re-allocation</th><td rowspan="14"> </td><th colspan="5">After re-allocation</th> </tr><tr><td>Team</td> <td>Call</td> <td>New</td> <td>Total</td> <td>Total Ward</td> <td>Team</td> <td>Call</td> <td>New</td> <td>Total</td><td>Total Ward</td></tr> <tr> <td>B1</td> <td> </td> <td>5</td> <td>11</td> <td rowspan="3">29</td> <td>B1</td> <td> </td> <td>8</td> <td>14</td> <td rowspan="3">32</td> </tr> <tr> <td>B2</td> <td> </td> <td>1</td> <td>10</td> <td>B2</td> <td> </td> <td>1</td> <td>10</td> </tr> <tr> <td>B3</td> <td>Yes</td> <td>1</td> <td>8</td> <td>B3</td> <td>Yes</td> <td>1</td> <td>8</td> </tr> <tr> <td>W1</td> <td>Yes</td> <td> </td> <td>4</td> <td rowspan="3">24</td> <td>W1</td> <td>Yes</td> <td> </td> <td>4</td> <td rowspan="3">24</td> </tr> <tr> <td>w2</td> <td> </td> <td>1</td> <td>10</td> <td>w2</td> <td> </td> <td>1</td> <td>10</td> </tr> <tr> <td>w3</td> <td> </td> <td>5</td> <td>10</td> <td>w3</td> <td> </td> <td>5</td> <td>10</td> </tr> <tr> <td>R1</td> <td> </td> <td>11</td> <td>16</td> <td rowspan="3">33</td> <td>R1</td> <td> </td> <td>8</td> <td>13</td> <td rowspan="3">31</td> </tr> <tr> <td>R2</td> <td> </td> <td> </td> <td>10</td> <td>R2</td> <td> </td> <td> </td> <td>11</td> </tr> <tr> <td>R3</td> <td>Yes</td> <td> </td> <td>7</td> <td>R3</td> <td>Yes</td> <td> </td> <td>7</td> </tr> <tr> <td>G1</td> <td> </td> <td> </td> <td>5</td> <td rowspan="3">24</td> <td>G1</td> <td> </td> <td> </td> <td>5</td> <td rowspan="3">23</td> </tr> <tr> <td>G2</td> <td> </td> <td> </td> <td>4</td> <td>G2</td> <td> </td> <td> </td> <td>4</td> </tr> <tr> <td>G3</td> <td> </td> <td>5</td> <td>15</td> <td>G3</td> <td> </td> <td>4</td> <td>14</td> </tr> </tbody></table>

---++ Initial Roster Generation

An initial mathematical model was developed in Python. This model simply allocates one shift to every registrar and ensures that the correct number of admitting shifts are allocated each day. Python outputs this as a color-coded schedule into HTML format. This initial model has been set up to generate a schedule for one individual registrar for the 12 weeks. This roster can then be allocated to different teams so that it roatates. For example Team 2 starts the roster in week 2, team 3 in week 3 ect. This incorporates very basic constraints and just allocates the shifts so that the required shifts will be allocated each day.
---+++ Initial Rotating Roster

<img alt="example_roster.bmp" src="%ATTACHURL%/example_roster.bmp" />

The above image shows part of the generated schedule. This schedule currently is not practical, for example it does not match up night shifts with days off. This model has been set up to generate a schedule for each individual registrar and to develop a schedule for one registrar which is rotated through by each of the registrars. The advantage of the rotating schedule is that it is a smaller problem to solve and thus should solve faster (hopefully). The rotating roster also means that the work load on all registrars is equal as over the 12 weeks they all do the same roster just in different weeks
---+++ Ward Numbers

It is useful to obtain an estimate of number of patients per ward, each day based on the roster. Code was developed to do this based on the assumptions listed below.

*Assumptions:*

   * 28 patients admitted Monday-Friday 
   * 1/4 of these patients go to each on call team 
   * 7 patients admitted Saturday and Sunday (This needs to be modified once we have a better number of patient admissions and discharge rates on the weekend) 
   * 1/3 of patients are discharged by each team each day *To include in future:* 

   * Need better data regarding patient admissions - awaiting ethics approval. Are numbers significantly different across the week? 
   * Work out legend and title code so graphs are a bit prettier. 
   * Compare out estimates for the current schedule with actual ward numbers to determine how accurate our estimating method is. <img alt="current_scedule_predicted_nums.png" src="%ATTACHURL%/current_scedule_predicted_nums.png" /> 

*Above:* Estimate of patient numbers per ward based on first 3 weeks of current roster.

*Below:* Estimate of patient numbers per ward based on initial rotating roster.

<img alt="patient_loading.png" src="%ATTACHURL%/patient_loading.png" />

---++ Next Steps

   * Look at historical data regarding patient admission and verify split of admissions across teams. 
   * Shadow a medical team to gain further understanding as to how the roster affects patient care. 
   * Mathematically formulate further constraints and incorporate in model 
   * Incorporate patient loading on weekends 

<font color="#630000">Future things to think about</font>

   * Scheduling public holidays 
   * Annual leave 
   * 'Blank slate' -- Currently Junior Doctor rules are under negotiation and thus it is possible that a roster may be able to be developed from scratch, focusing on patient needs. Stephen Child spoke about this at the meeting 03/12/2010. 
   * If patients are returning patients then they get realloctaed to the previous team they were seen by? 
---++ The to-do list

This section to keep track of tasks and progress as it is made - area for my general thought dump.

<table width="80%" cellspacing="2" cellpadding="2" border="1"> <tbody><tr><th>Task</th> <th>Description</th> <th> Date Completed</th> </tr> <tr><td>Comment Code</td> <td>This really needs to be kept up-to-date and requires full descriptions of each functions. This will ensure that it is easy for people to continur working on in the future</td><td> </td></tr> <tr><td>List of what is currently included in schedule and what to include next</td><td> define where we are at now for and what to include next in the formulation</td><td> </td> </tr><tr><td>Patient number estimation</td><td>

   * Compare our estimate from a previous roster with data from that roster 
   * Get more data on admission and discharges in weekend and include in code 
</td><td> </td></tr> <tr><td>

Analyse previous data

</td><td> </td><td> </td></tr> </tbody></table>

---++++ What's in schedule and what isn't 15/12/2010 Updated 21/12/2010

*Whats in:*
   * Shifts A B C D allocated to one registrar on weekdays 
   * 2 P shifts and one A shift allocated on weekends 
   * N shifts allocated sun-Thurs 
   * Correct order of shifts for weekend 

*What to add next:*

   * Rules regarding whoever does A on Friday does P on the weekend --- added 21/12/2010 
   * Get rid of night shift but add relief roster constraints - just add to first week? --- added 21/12/2010 
   * Only one on call shift allocated to each ward? -- 21/12/2010 - may not be needed 
   * look at just minimizing number of patients -- have run down to 32 patients won't solve for lower than 32 <font color="#630000"> </font> <font color="#630000"> </font> <font color="#630000"> </font> 
---++ <font color="#630000"> Rosters produced by additional constraints </font>

<font color="#630000"> The roster below was produced by adding an additional constraint that a team can do no more than 2 'long days' in any 7 day period </font>

<font color="#630000"> <img alt="roster_max_35.bmp" src="%ATTACHURL%/roster_max_35.bmp" /> </font>

<font color="#630000"> </font>

<font color="#630000"> </font>

<img alt="max_35_nor_more_than_2_long.png" src="%ATTACHURL%/max_35_nor_more_than_2_long.png" />
---++ Holiday 23rd December to 12th January

Wonderful holiday worked on tan. result great tan!
---++ 12th Jan onwards:
---+++ What to do next from meeting with Mike and Cam 12/01/2011:

The stage we are at now is that the LP is producing a roster with the constraints of ensuring the union rules are met and that there is a maximum number of 35 patients per ward. This appears to be better than the current roster in regards to patient numbers.

<img width="1133" alt="" src="%ATTACHURL%/comparing_old_with_new_schedule_graphs.JPG" height="525" />

%BLUE%<u><br /></u>%ENDCOLOR%

Analyse data from ADHB. To do this I need to write some VBA to match up the number of patients admitted to each ward and compare this with our method of predicting the number of patients in each ward. We then want to compare our roster with the current roster using the number of patients from the data using our current assumptions. Following on from this we want to compare our assumptions regarding how the patients are allocated and amount discharged.

Goal of completing these by end of January. Following on into February we will review what has been done so far and read union documents to check for additional rules. Also should probably look at removing fractional number of patients in calculation or may need to look at impact of this assumption.
---++ Data Analysis

In the period 14/06/2010-12/12/2010 for the data given there were 3100 patients admitted and 6070 discharged. This is clearly inconsistent and so I used the admitting numbers which seemed o.k as it was an average of 17 patients admitted per day. When I visited the hospital the number of patients was between 18-2 and the wards have around 20 bed this seemed a reasonable figure. I used the discharge assumption that 1/3 of patients were discharged on Monday to Friday and 1/5 on Saturday and Sunday. These rates can be easily adjusted when we obtain more reliable data. Each ward starts the period with 18 patients in where previously there were between 18-23. This enables a lower maximum number of patients in the formulation and thus should spread the patient load more evenly. The maximum number of patient from the formulation is 19. Interestingly this formulation one 'on call' shift to each ward each day as was with the previous roster for June to December 2010. I'm not sure why this was changed for the roster from 13th December onwards.

[[%ATTACHURL%/data_analysis.pdf][data_analysis.pdf]]

Feedback from ADHB regarding the data that was supplied didn't include the number transfered into GM. Rupert supplied a summary of data from 2010 broken of patient admissions / discharges which included transfers. This data balances and hopefully we will be able to source similar data broken into a day by day format.

%BLUE%<u> [[%ATTACHURL%/ADT23_ArtN_Summary_data.pdf][ADT23_ArtN_Summary_data.pdf]]</u>%ENDCOLOR%
---++ Objective function

From the graphs of patient numbers from the current rosters there is still quite a big gap between patient numbers. There is the option of doing a two step optimization process, where first minimize the maximum number of patients, set this minimum as a constraint and then re-run maximizing the minimum number of patients in the ward. The better option is likely to be minimize the difference between these two numbers.

Using the data given averages were calculate for each day of the week. Admitting seems constant through the week with about 10 patients per day. This is not including internal transfer. Discharges are higher on Fridays and lower on the weekends. Average patient stay is 1 day 7 hours.
| Day | Admitting | Discharge |
| Monday | 10.21 | 39.31 |
| Tuesday | 10.60 | 39.08 |
| Wednesday | 9.50 | 38.19 |
| Thursday | 10.62 | 36.50 |
| Friday | 10.02 | 44.81 |
| Saturday | 4.73 | 21.19 |
| Sunday | 3.94 | 14.38 |

Trying to solve with objective seems smallest gap it can get is around 25 patients which seems quite big, although this is a smaller gap than when solving before. It seems to be definitely problem with only one on-call person being on weekends or could be assumptions regarding patient data. The current assumption is 28 admitted on Saturday and Sunday. This seems way too high when compared with the data above. I modified the input data to 20 admissions on weekdays and 10 on Saturday and 8 on Sunday. I have doubled these numbers from the admitting numbers above to model the transfers as well.

I tried to solve with the changes to the data and it is struggling. With running to a 0.05 and 0.1 tolerance between integer solution and best solution, Pulp crashes. I removed this tolerance and have tried to solve it again. After 1429.35 seconds (about 23 minutes) the integer solution is 10.6798 and the best possible is 7.8805, this is about what it was getting to before it crashed.

Image of patient numbers and roster: 

[[%ATTACHURL%/MinGap1.JPG][MinGap1.JPG]]

[[%ATTACHURL%/solution_table_rotate_full.html][solution_table_rotate_full.html]] 

I added a constraint to say that the difference between the maximum and minimum number of patients must be >=10 to try and get it to solve as it seems that the best gap must be about 10.

I am also considering that we are only 'smoothing' out patient numbers across wards and not across registrars in the wards. Adding these in as variables I suspect will increase the solve time but it could be something to investigate.

---++ References
   * [[https://matplotlib.sourceforge.net]]

-- Main.MichaelOSullivan - 20 Jan 2012
Attachments
Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
pdf	ADT23_ArtN_Summary_data.pdf	r1	manage	24.1 K	2011-01-21 - 01:16	AmeliaWhite	2010 summary data by team
rtf	ADT23_ArtN_Summary_data.rtf	r2 r1	manage	24.1 K	2011-01-21 - 01:14	AmeliaWhite	Summary data from 2010 broken down by team
jpg	MinGap1.JPG	r1	manage	129.9 K	2011-01-26 - 22:14	AmeliaWhite	Graph from minimising the gap between maximum and minimum number of patients. Stopped solver at 1736.60 seconds
docx	Patient_number_estimate.docx	r1	manage	191.5 K	2010-12-15 - 04:47	AmeliaWhite	Graphs of estimate of patient numbers to show to ADHB at meeting
Topic revision: r1 - 2012-01-20 - MichaelOSullivan
~~Edit~~
~~Attach~~