We built randomly cross-linked model poly(dimethylsiloxane) (PDMS) networks and used molecular dynamics methods to obtain stress-strain curves. The Mooney-Rivlin (MR) analysis was used to estimate the shear moduli. We applied primitive path analysis (PPA) and its variation, phantom primitive path analysis (3PA), to estimate the entanglement and the cross-link moduli, respectively. The MR moduli estimates are in good agreement with the sum of the entanglement and the cross-link moduli, and we observe that the stress-strain data collapse to a universal form when reduced with the PPA and 3PA moduli. We studied how the MR parameters C1 and C2 vary from cross-link to entanglement-dominated networks. For the latter, we observed 40% and 60% contributions of 2C1 and 2C2 to the shear modulus, respectively. Finally, we fitted several theoretical models to our stress-strain data. While all fits are good, the estimates for the entanglement and cross-link moduli vary significantly when compared to our independently obtained PPA and 3PA benchmarks.