产生于某任意时刻W前的TR天W内。因此，在血红蛋白动态的稳定状 态下，某任意时刻下存在于血液中的血红蛋白的总量（tota化g)可W用W下的式（1)表示。
[0032] totalHg=TRXG (1)
[003引另一方面，血液中的血红蛋白浓度（CHg)，可W作为上述的血红蛋白总量除朗t定 为体重的8% (男性的情况）或7% (女性的情况）的血液量度V)得到的值用下式（2)表 /J、- 0
[0034]CHg=totalHg/BV (2)
[0035] 而若向上述式似代入前述式（1)则可W得到W下的式（3a)。
[0036]CHg=TRXG/BV (3a)
[0037] 另外，若改写该式（3a)则可W得到W下的式（3b)。
[0038]G=CHgXBV/TR (3b)
[0039] 因此，若与所产生的血红蛋白量相比，所去除的血红蛋白量多或者少的情况下，即 血红蛋白动态为非稳定状态的情况下，血红蛋白浓度降低或者升高。此时，通常所去除的血 红蛋白量增减的结果，血红蛋白浓度不会降低或者升高，所产生的血红蛋白量增减的结果， 血红蛋白浓度升高或者降低。
[0040] 因此，若将每一天的血红蛋白产生量定义为血红蛋白产生速度，则血红蛋白动态 为非稳定状态、血红蛋白浓度变动的情况下，两时刻期间的生物体内的总血红蛋白量之差 为两时刻期间的血红蛋白产生速度之差累积而产生的。因此，若将两时刻期间的生物体内 的总血红蛋白量之差除W两时刻期间的天数则得到两时刻期间的血红蛋白产生速度之差， 该关系用W下的式（4)表示。需要说明的是，W下的式（4)中，AG表示两时刻的血红蛋白 产生速度之差，CHgl表示某任意时刻的血红蛋白浓度，CHg2表示在该时刻之前测得的血红 蛋白浓度，而BV表示血液量、TT表示该两时刻期间的天数。
[0041]ΔG=(CHglXBV-CHg2XBV)/TT (4)
[0042] 而上述式（4)中，血红蛋白产生速度日益增加时，ΔG〉0,血红蛋白产生速度日益 减少时，ag<0。也就是说，与两时刻的血红蛋白浓度相等的情况（血红蛋白动态处于稳 定状态的情况）相比，两时刻的血红蛋白浓度不同的情况（血红蛋白动态处于非稳定状 态的情况）下，血红蛋白产生速度仅增加AG。因此，非稳定状态中的血红蛋白产生速度 (Greal)，通过利用前述式（3b)求出的G加上AG而成的式巧a)表示。
[0043]Greal=G+ΔG巧a)
[0044] 另外，若向该式巧a)代入前述式（3b)和式（4)则可得到求出非稳定状态下的血 红蛋白产生速度的式巧b)。需要说明的是，为了使得成为基准的时刻清楚，式（3b)的血红 蛋白浓度改写为CHg2。进而，在此TR表示红细胞寿命、TT表示该两时刻期间的天数。
[0045]Greal=CHg2XBV/TR+(CHglXBV-CHg2XBV)/TT 巧b)
[0046] 另一方面，向体内给予ESA的结果发现，所产生的血清ESA浓度如图1所示那样指 数函数地减少。即，ESA浓度的衰减曲线通过式（6)示出。需要说明的是，在此，b为常数， t为给予ESA后的经过时间。另外，CesA(t)表示时间t时的ESA浓度、CesA(O)表示刚给予 ESA之后的时刻、即t= 0时的ESA浓度。
[0047] Cesa(t) =Cesa(0)Xexp(_bXt) 巧）
[004引其中，Cesa(O) =Desa/Vesa、Desa为ESA的给药量、Vesa为ESA的分布体积。常数b根据ESA的种类使用由已经报告的数据求出的值，例如达贝泊汀α的情况下，求出b= 0. 40943,另外，依伯汀α、依伯汀β的情况下求出b= 1. 2053。
[0049] 对于ESA的给药频率而言，通常依伯汀α、依伯汀e，1周3次、1月12~15次、 每隔2~3天给药，达贝泊汀α，1周1次、1月4~5次、每隔7天给药。而若各ESAW运 种时间间隔给药则任意时刻的ESA给药时，此前的ESA给药伴随的血清ESA浓度降低至大 致零。因此若例如由前次的血红蛋白测定时刻直至当前血红蛋白测定时刻为止期间给予 η次的ESA，则由前次的血红蛋白测定时刻直至当前血红蛋白测定时刻为止期间的作为ESA 浓度的积分平均值的平均ESA浓度（MCesA)通过W下的式（7a)或（7b)算出。
[0050]
[005引其中，CesA的下标（1、2、…η)表示ESA的给药顺序，TD表示ESA的给药间隔。另 一方面，若Desa设为ESA的给药量、VESA设为ESA的分布体积则Cesa(O) =Desa/Vesa。若将其 代入式（7b)则得到W下的式（8a)。
[0054]
[005引其中，式（8a)中，DesA的下标（1、2、…η)表示ESA的给药顺序。另外，作为ESA的 分布体积的VesaW体重的约5%、具体而言对于男性而言体重的5. 2%、对于女性而言体重 的4.6%表示。
[0056] 接着，若改写该式（8a)则得到W下的式（8b)。
[0057] MCesa_ {D ESAl+〇Ea化· · · · · +DesaD }
[0058]X{l-exp(-bXTD) } /b/VesA/TT (8b)
[0059] 而若利用体外（invitro)的实验则培养红细胞的溶液中的血红蛋白产生速度 (G)例如如图2所示那样与ESA浓度（CeJ的对数值成比例（参照永野伸郎等、…肾和透析 (腎透析）"、60化)，1039-1046, 2006)，可W利用下述的式（9)表示。其中，在此，a为常 数，作为血红蛋白产生对于ESA浓度的灵敏度认识。
[0060] G=aXln(CESA) 巧）
[0061] 需要说明的是，上述的式（9)中，若ESA浓度接近零则ESA浓度的对数值向着负无 限大降低，因此，G也向着负无限大降低。因此，该式存在在极低的ESA浓度时并非有效的 缺点。
[0062] 特别是上述血红蛋白产生速度与ESA浓度的关系，在患者的血清ESA浓度的通常 的范围内，也能够近似于直线性比例关系。若血红蛋白产生速度与血清ESA浓度之间的关 系近似于直线性比例关系，则在血清ESA浓度的通常的范围内，G的计算值产生一些误差， 但是运具有即使是低的ESA浓度也是有效的优点。
[0063] 而式（9)在生物体内也成立，因此若将先前的平均ESA浓度适用于上述式（9)则 可W得到下式（10a)。
[0064]G=aXln(MCESA)(10a)
[0065] 进而若改写该式（10a)则可W导出提供作为血红蛋白产生对于ESA浓度的灵敏度 的a的式（1化）。
[0066]a=6/1π(Μ〇ε5α) (l〇b)
[0067] 需要说明的是，上述式（10b)中，a根据患者另外根据患者的状态变化，因此测定 血红蛋白浓度时，通过最新的血红蛋白浓度和前次的测定时的血红蛋白浓度、自前次的血 红蛋白测定时刻直至当前血红蛋白测定时刻为止期间的平均ESA浓度算出上述a。
[0068] 而对于透析患者而言，在每约1个月（每4~5周）开始透析时测定血红蛋白浓 度，而对于本发明而言，通过如此测得的最新的血红蛋白浓度、前次（1个月前）测定时的血 红蛋白浓度、和自前次测定时直至此次测定时为止期间的1个月给予的ESA量，确定在下次 的采血时可得到目标血红蛋白产生速度的ESA给药量，并将其给予患者。于是，为了确定在 下次采血时得到目标的血红蛋白产生速度的ESA给药量，使用下述方法。
[0069] 目标血红蛋白浓度稳定持续的状态下，血液中的血红蛋白全部在过去90天W均 匀速度产生，存积于血液中。因此目标血红蛋白产生速度，通过作为目标血红蛋白浓度与血 液量的乘积即目标总血红蛋白量除W作为红细胞的平均寿命的90天来得到。目P，若作为下 次采血时想要存在于生物体内的血红蛋白量的目标血红蛋白浓度设为targetCHg、目标血 红蛋白产生速度设为targetG、血液量设为BV、红细胞的平均寿命设为TR( = 90天），则目 标血红蛋白产生速度通过W下的式（11)表示。
[0070]targetG=targetCHgXBV/TR (11)
[0071] 另一方面，若能够达成目标血红蛋白浓度的平均ESA浓度设为targetMCESA、将前 述式（10a)适用于targetMCesA和targetG,则可得到W下的式（12a)。
[0072]targetG=aXIn(targetMCesA) (12a)
[0073] 若改写上述式（12a)则可得到由目标血红蛋白浓度算出实现其的平均ESA浓度的 式（1化）。
[0074] targetMCESA=exp(targetG/a) (12b)
[00巧]其中，a由当前时刻的血红蛋白产生速度和基于与此相比在之前进行的ESA的给 药的平均ESA浓度，使用前述式（10b)算出。
[0076] 接着若向式（12b)代入式（11)，则可得到由目标血红蛋白浓度求出达成其的平均 ESA浓度的式（13)。
[0077]targetMCESA=exp(targetCHgXBV/TR/a) (13)
[0078] 在此，通过式