Sampling delay times from empirical distributions by hour of the day/week

Question

question

Maryam H2 asked Dec 4, '23 Maryam H2 commented Jan 27, '24

Sampling delay times from empirical distributions by hour of the day/week

How can I read the delay time from an empirical distribution for each hour of the week? The data is organized with varying sample sizes for each hour, ranging from 0 (Sunday 12 AM) to 48 (Monday 12 AM). Each sample corresponds to a delay time, and we want to sample delay times from this empirical distribution for each hour of the week.

It's possible to separate days of the week but then is it possible/easier to read the data for each weekday by the hour from an empirical distribution, again at each hour we have an unequal number of samples and we want to sample from this distribution.

The data is in the GlobalTable1 in the model attached here.

Thanks

samplingbyhouroftheweek.fsm

Software Version:

FlexSim 23.2.0

flexsim healthcare by hour emprical sampling

samplingbyhouroftheweek.fsm (54.3 KiB)

5 |100000

Attachments: Up to 12 attachments (including images) can be used with a maximum of 23.8 MiB each and 47.7 MiB total.

Answer 1 · 2023-12-04T19:28:29Z

Jeanette F answered Dec 4, '23 Maryam H2 commented Jan 27, '24

Hello @Maryam H2,

An Empirical Distribution for an hour of a day will need to be its own table, so from your data we have created a Table for each as a Cell in a table.

The Table's first column needs to be the weight distribution and the second column is the value. We assumed an even distribution and calculated it so. Here is the Table for Hour of Week 0;

You can then reference that Cell of the Table for the Empirical Distribution for that hour.

empirical(Table("GlobalTable2").cell(Model.dateTime.hour, 2))

samplingbyhouroftheweek_1.fsm

1701717681496.png (5.9 KiB)

1701718029508.png (12.6 KiB)

samplingbyhouroftheweek-1.fsm (60.2 KiB)

· 10

5 |100000

Attachments: Up to 12 attachments (including images) can be used with a maximum of 23.8 MiB each and 47.7 MiB total.

Maryam H2 commented · Dec 05, 2023 at 02:42 PM

Hi @Jeanette F

Thanks for the answer, I have two questions:

1- Sometimes when I use the script to check the value samples by the code below, it returns something that is not in that hour's table. Do you know why that is?

For example, to check the value of samples for the first hour.

empirical(Table("GlobalTable2").cell(1, 2)

2- In the other script, where you created a table and cloned it to a new global table, what is the $1 does here?

..FROM GlobalTable1 WHERE [Hour of the Day] = $1"

0 ·

Jason Lightfoot ♦♦ Maryam H2 commented · Dec 05, 2023 at 03:00 PM

1. Should be Model.dateTime.hour+1 since it's using the number of hours passed rather than the current hour the model is in.

2. That's part of a loop where the $1 is a parameter for the hour of the day, which changes for each cycle through the loop.. I'm surprised by this - I thought you were going to use hour of the week.

0 ·

Maryam H2 Jason Lightfoot ♦♦ commented · Dec 06, 2023 at 06:50 PM

got it, thanks!

In my model, I seperated days to have their own table, so for each week of the day I have a Global Table and that's why I used Hour of the Day column.

0 ·

Show more comments

question