We present results for the QCD equation of state with nonzero chemical potential using the Taylor expansion method with terms up to sixth order in the expansion. Our calculations are performed on asqtad 2+1 quark flavor lattices at $N_t=4$.