The stellar distribution of a small sample of edge-on spiral galaxies is examined in B, V, R, and I by fitting model distributions to the light profiles, both perpendicular to the galaxy planes and along the major axes. We have developed a method to compare the fits for the models obtained for different galaxies by introducing a global goodness-of-fit parameter P-ij. In general, the best fitting vertical model is more peaked than expected for an isothermal sheet distribution, i.e. it is either an exponential light distribution or a sech(z)-model. The most likely explanation for the deviations from an exponential or a sech(z) distribution at small z is dust extinction, whereas the excess light detected at large z for a few early-type spirals...