1. A method for treating anemia in a subject using erythropoietic stimulating agents (Epo), comprising: obtaining a target hemoglobin value (Hb target ), a current hemoglobin value (Hb k ), a rate of change of hemoglobin value (ΔHb k ), and an Epo dose received value (Epo k ) related to an amount of Epo being received by the subject; determining an initial recommended Epo dose (Epo k+1,i *) for each of i dosing regimens based on the target hemoglobin value (Hb target ), the current hemoglobin value (Hb k ), the rate of change of hemoglobin value (ΔHb k ), and the previous Epo dose prescribed value (Epo k ); determining a next recommended Epo dose (Epo k+1 ) based on the initial recommended Epo doses for the i dosing regimens weighted by degrees of matching computed using fuzzy membership functions; determining, using a dose-response model, a Hb k+1 value and a predicted hemoglobin response value (ΔHb k+1 ) for a subject with a respective Epo sensitivity profile based on the Hb k value, the ΔHb k value, and a proposed Epo k+1 value; determining, using an optimization algorithm, the Epo k+1,i * value by determining the proposed Epo k+1 value using the Hb target value, the Epo k value, the Hb k+1 value and the ΔHb k+1 value, and then iteratively re-determining the proposed Epo k+1 value until the proposed Epo k+1 values converge to the Epo k+1,i * value; and administering the next recommended Epo dose to the subject; wherein the dose-response model is substantially according to the function: [ dHb dt d ⁢ ⁢ Δ ⁢ ⁢ Hb dt ] = [ - 2 T i 1 T i 2 1 0 ] ⁡ [ Hb Δ ⁢ ⁢ Hb ] + [ K i T i 2 0 ] ⁢ Epo where: K i —erythropoietic response (1,000 Units Epo per 1 g/dL Hb change) T i —time constant Epo—dose (1,000 Units) Hb—hemoglobin level (g/dL) ΔHb—hemoglobin rate of change (g/dL per day); and wherein the optimization algorithm is substantially according to the function: OFV = ∑ k p = 1 H p ⁢ ( Hb target - Hb k p ) 2 + λ i ⁢ ∑ k c = 1 H c ⁢ Δ ⁢ ⁢ Epo k c 2 where: k p , k c —time steps (weeks) H p —prediction horizon (weeks) H c —control horizon (weeks) λ i —dose change suppression (non-dimensional) Hb kp —hemoglobin at step k p ΔEpo kc —change in dose from step k c-1 to k c ; Subject to constraints: 0≦Epo≦90,000 units per week 0≦Hb≦20 g/dL −0.5<ΔHb<0.5 g/dL per week where: ΔHb=Hb kp −Hb kp-1 .