forked from freesurfer/freesurfer
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathestimateIOBTable.m
83 lines (78 loc) · 5.75 KB
/
estimateIOBTable.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
activeInsulin_3Hr = ...
... % 0 1 2 3 4 5 6 7 8 9 - ending minute digit
[ 10000, 10000, 10000, 10000, 10000, 10000, 9996, 9984, 9965, 9944, % 0 - 9
9920, 9895, 9876, 9847, 9816, 9783, 9760, 9723, 9685, 9645, % 10 - 19
9603, 9574, 9529, 9483, 9435, 9402, 9351, 9299, 9245, 9190, % 20 - 29
9152, 9094, 9035, 8974, 8933, 8870, 8805, 8740, 8673, 8628, % 30 - 39
8559, 8488, 8417, 8369, 8296, 8221, 8146, 8070, 8018, 7940, % 40 - 49
7861, 7782, 7728, 7647, 7565, 7482, 7398, 7342, 7257, 7172, % 50 - 59
7086, 7028, 6941, 6854, 6766, 6677, 6618, 6529, 6440, 6350, % 60 - 69
6290, 6200, 6109, 6019, 5928, 5868, 5777, 5687, 5596, 5536, % 70 - 79
5445, 5355, 5265, 5175, 5115, 5025, 4936, 4847, 4788, 4700, % 80 - 89
4613, 4526, 4439, 4382, 4296, 4211, 4127, 4071, 3988, 3906, % 90 - 99
3824, 3743, 3690, 3611, 3532, 3455, 3404, 3328, 3253, 3179, % 100 - 109
3106, 3058, 2987, 2917, 2848, 2803, 2736, 2670, 2605, 2541, % 110 - 119
2499, 2438, 2377, 2318, 2279, 2222, 2165, 2110, 2057, 2021, % 120 - 129
1970, 1919, 1869, 1837, 1789, 1743, 1697, 1653, 1624, 1582, % 130 - 139
1540, 1500, 1473, 1435, 1397, 1360, 1324, 1301, 1266, 1232, % 140 - 149
1200, 1178, 1147, 1116, 1086, 1057, 1037, 1009, 981, 954, % 150 - 159
936, 910, 884, 858, 833, 817, 792, 768, 744, 728, % 160 - 169
705, 682, 659, 636, 621, 598, 575, 553, 538, 516, % 170 - 179
493, 0, 0, 0, 0, 0, 0, 0, 0, 0 % 180 - 189
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, % 190 - 199
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, % 200 - 209
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, % 210 - 219
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, % 220 - 229
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, % 230 - 239
0, 0, 0, 0, 0, 0, 0, 0, 0, 0]; % 240 - 249
activeInsulin_3Hr = activeInsulin_3Hr /100;
activeInsulin_4Hr = ...
... % 0 1 2 3 4 5 6 7 8 9 - ending minute digit
[ 10000, 10000, 10000, 10000, 10000, 10000, 10000, 10000, 9996, 9984, % 0 - 9
9972, 9958, 9944, 9928, 9912, 9895, 9876, 9857, 9837, 9805, % 10 - 19
9783, 9760, 9736, 9711, 9685, 9658, 9631, 9603, 9574, 9545, % 20 - 29
9514, 9467, 9435, 9402, 9368, 9334, 9299, 9263, 9227, 9190, % 30 - 39
9152, 9113, 9074, 9015, 8974, 8933, 8891, 8848, 8805, 8762, % 40 - 49
8718, 8673, 8628, 8582, 8535, 8465, 8417, 8369, 8320, 8271, % 50 - 59
8221, 8171, 8121, 8070, 8018, 7966, 7914, 7835, 7782, 7728, % 60 - 69
7674, 7619, 7565, 7509, 7454, 7398, 7342, 7286, 7229, 7143, % 70 - 79
7086, 7028, 6970, 6912, 6854, 6795, 6736, 6677, 6618, 6559, % 80 - 89
6499, 6410, 6350, 6290, 6230, 6170, 6109, 6049, 5989, 5928, % 90 - 99
5868, 5807, 5747, 5656, 5596, 5536, 5475, 5415, 5355, 5295, % 100 - 109
5235, 5175, 5115, 5055, 4996, 4907, 4847, 4788, 4730, 4671, % 110 - 119
4613, 4555, 4497, 4439, 4382, 4324, 4268, 4183, 4127, 4071, % 120 - 129
4015, 3960, 3906, 3851, 3797, 3743, 3690, 3637, 3584, 3506, % 130 - 139
3455, 3404, 3353, 3303, 3253, 3203, 3154, 3106, 3058, 3010, % 140 - 149
2963, 2894, 2848, 2803, 2758, 2713, 2670, 2626, 2584, 2541, % 150 - 159
2499, 2458, 2417, 2357, 2318, 2279, 2241, 2203, 2165, 2129, % 160 - 169
2092, 2057, 2021, 1987, 1953, 1902, 1869, 1837, 1805, 1774, % 170 - 179
1743, 1712, 1682, 1653, 1624, 1596, 1568, 1527, 1500, 1473, % 180 - 189
1447, 1422, 1397, 1372, 1348, 1324, 1301, 1277, 1255, 1221, % 190 - 199
1200, 1178, 1157, 1136, 1116, 1096, 1076, 1057, 1037, 1018, % 200 - 209
1000, 972, 954, 936, 919, 901, 884, 867, 850, 833, % 210 - 219
817, 800, 784, 760, 744, 728, 713, 697, 682, 666, % 220 - 229
651, 636, 621, 605, 590, 568, 553, 538, 523, 508, % 230 - 239
% 0 1 2 3 4 5 6 7 8 9 - ending minute digit
493, 0, 0, 0, 0, 0, 0, 0, 0, 0]; % 240 - 249
activeInsulin_4Hr = activeInsulin_4Hr /100 ;
[rows,cols] = size(activeInsulin_3Hr) ;
ind = 1 ;
for row=1:rows
for col=1:cols
iob_table3(ind) = activeInsulin_3Hr(row,col) ;
ind = ind + 1 ;
end
end
[rows,cols] = size(activeInsulin_4Hr) ;
ind = 1 ;
for row=1:rows
for col=1:cols
iob_table4(ind) = activeInsulin_4Hr(row,col) ;
ind = ind + 1 ;
end
end
iob_table = iob_table3 ;
iob_change_table(1) = 0 ;
for ind=2:length(iob_table)
iob_change_table(ind) = iob_table(ind-1)-iob_table(ind) ;
end