Incorporating tumor stage or the PAM50 subtypes separately as blocking variables into each of the generalized linear models described above, revealed that tumor stage did not confound the observed interaction differences of ER and MYC between IBC and nIBC. Here, ESR1 is linked to neoplasm.