The optical properties of organic chromophores, in the class of nonlinear optical molecules, have improved dramatically. These chromophores have large hyperpolarizability and large permanent dipole moments, in excess of 10 D. When these molecules are mixed with a host polymer at high chromophore density and processed using an aligning electric poling field, the resulting material does not have the electro-optic (EO) response promised by the individual molecules. We have used equilibrium statistical mechanics, and Monte Carlo (MC) methods in a lattice-based model, to investigate how the molecular properties translate into the EO response of the material. We compare our MC results with a previously developed analytic theory and find good agreement when certain assumptions about the analytic theory are made. By use of the MC methods to compute order parameters for chromophores at high density on a lattice in the presence of an external uniform poling field, we find that greater order and hence greater EO performance can be obtained by adjusting the geometry of the lattice. This suggests that modifications to existing chromophores could lead to improvements in the EO properties of the materials. We compare with other lattice and mean field models. In our lattice model, we find no ferroelectric phase transition.