A method for determining the likelihood of breast cancer recurrence in a lymph node-negative, invasive breast cancer subject comprising: (a) determining the expression levels of the RNA transcripts of Grb7, Her2, ER, PR, BC12, CEGP1, SURV, Ki.67, MYBL2, CCNB1, STK15, CTSL2, STMY3, GSTM1, BAG1 and CD68, or corresponding substitute genes, in a biological sample containing tumor cells obtained from said subject; (b) creating the following gene subsets: (i) growth factor subset: GRB7 and HER2 (ii) estrogen receptor subset: ER, PR, BCl2 and CEGP1 (iii) proliferation subset: SURV, Ki.67, MYBL2, CCNB1, and STK15; and (iv) invasion subset: CTSL2, and STMY3; wherein a gene within any of subsets (i)-(iv) can be substituted by a substitute gene which co-expresses with said gene in said tumor with a Pearson correlation coefficient of > 0.40 by determining the range of values for the gene within any of subsets (i)-(iv) and the substitute gene in the patient population and applying a linear transform to map the values of the substitute gene into the corresponding values for the gene within any of subsets (i)-(iv); and (c) calculating the recurrence score for said subject by the following equation: RS = 0.23 to 0.70 × GRB7axisthresh − 0.17 to 0.51 × ERaxis + 0.52 to 1.56 × prolifaxisthresh + 0.07 to 0.21 × invasionaxis + 0.03 to 0.15 × CD 68- 0.04 to 0.25 × GSTM 1 − 0.05 to 0.22 × BAG 1, wherein GRB7 axis = (0.45 to 1.35) x GRB7 + (0.05 to 0.15) x HER2; if GRB7 axis < -2, then GRB7 axis thresh = -2, and if GRB7 axis ≥ -2, then GRB7 axis thresh = GRB7 axis; ER axis = (ER + PR + Bcl2 + CEGP1)/4; prolifaxis = (SURV + Ki.67 + MYBL2 + CCNB1 + Stk15)/5; if prolifaxis < -3.5, then prolifaxisthresh = -3.5, if prolifaxis ≥ - 3.5, then prolifaxishresh = prolifaxis; and invasionaxis = (CTSL2 + STMY3)/2, wherein the individual contributions of the genes in, and are weighted by a factor of 0.5 to 1.5, and wherein a higher RS represents an increased likelihood of breast cancer recurrence. (a) determining the expression levels of the RNA transcripts of Grb7, Her2, ER, PR, BC12, CEGP1, SURV, Ki.67, MYBL2, CCNB1, STK15, CTSL2, STMY3, GSTM1, BAG1 and CD68, or corresponding substitute genes, in a biological sample containing tumor cells obtained from said subject; (b) creating the following gene subsets: (i) growth factor subset: GRB7 and HER2 (ii) estrogen receptor subset: ER, PR, BCl2 and CEGP1 (iii) proliferation subset: SURV, Ki.67, MYBL2, CCNB1, and STK15; and (iv) invasion subset: CTSL2, and STMY3; wherein a gene within any of subsets (i)-(iv) can be substituted by a substitute gene which co-expresses with said gene in said tumor with a Pearson correlation coefficient of > 0.40 by determining the range of values for the gene within any of subsets (i)-(iv) and the substitute gene in the patient population and applying a linear transform to map the values of the substitute gene into the corresponding values for the gene within any of subsets (i)-(iv); and (i) growth factor subset: GRB7 and HER2 (ii) estrogen receptor subset: ER, PR, BCl2 and CEGP1 (iii) proliferation subset: SURV, Ki.67, MYBL2, CCNB1, and STK15; and (iv) invasion subset: CTSL2, and STMY3; (c) calculating the recurrence score for said subject by the following equation: RS = 0.23 to 0.70 × GRB7axisthresh − 0.17 to 0.51 × ERaxis + 0.52 to 1.56 × prolifaxisthresh + 0.07 to 0.21 × invasionaxis + 0.03 to 0.15 × CD 68- 0.04 to 0.25 × GSTM 1 − 0.05 to 0.22 × BAG 1, wherein GRB7 axis = (0.45 to 1.35) x GRB7 + (0.05 to 0.15) x HER2; if GRB7 axis < -2, then GRB7 axis thresh = -2, and if GRB7 axis ≥ -2, then GRB7 axis thresh = GRB7 axis; ER axis = (ER + PR + Bcl2 + CEGP1)/4; prolifaxis = (SURV + Ki.67 + MYBL2 + CCNB1 + Stk15)/5; if prolifaxis < -3.5, then prolifaxisthresh = -3.5, if prolifaxis ≥ - 3.5, then prolifaxishresh = prolifaxis; and invasionaxis = (CTSL2 + STMY3)/2, wherein the individual contributions of the genes in, and are weighted by a factor of 0.5 to 1.5, and wherein a higher RS represents an increased likelihood of breast cancer recurrence. GRB7 axis = (0.45 to 1.35) x GRB7 + (0.05 to 0.15) x HER2; if GRB7 axis < -2, then GRB7 axis thresh = -2, and if GRB7 axis ≥ -2, then GRB7 axis thresh = GRB7 axis; ER axis = (ER + PR + Bcl2 + CEGP1)/4; prolifaxis = (SURV + Ki.67 + MYBL2 + CCNB1 + Stk15)/5; if prolifaxis < -3.5, then prolifaxisthresh = -3.5, if prolifaxis ≥ - 3.5, then prolifaxishresh = prolifaxis; and invasionaxis = (CTSL2 + STMY3)/2,