In this work, we use a stochastic-Chebyshev expansion approach to simulated the linear spectroscopic signatures of exciton-polariton under varying types of molecular frequency disorder: (I) Gaussian, (II) rectangular, and (III) Lorentzian. We find that the transmission spectra exhibit drastically varied characteristics between different electronic disorder types as well as, for a fixed type of disorder, between one-, two-, and many-molecule systems. Specifically, for the Gaussian and rectangular disorders, we found that the rate of collapse of the upper and lower polariton spectral features increases with an increasing number of molecules from N = 1 to 2. Further, the spectral features at weak disorders are drastically changing with increasing numbers of coupled molecules in both the Gaussian and rectangular disorders. At all numbers of coupled molecules, N, both the Gaussian and rectangular disorder types exhibit an increased Rabi splitting at weak disorder prior to its collapse, which agrees with previous analytical work in the large-N limit .