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 C-1 and C-2 vary from cross-link to entanglement-dominated networks. For the latter, we observed 40% and 60% contributions of 2C(1) and 2C(2) 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.