The purpose of this study was to investigate the application of the Monte Carlo technique to the calculation and analysis of output factors for electron beams used in radiotherapy. The code EGS4/BEAM was used to obtain phase-space files for 6, 12 and 20 MeV clinical electron beams from a scattering-foil linac (Varian Clinac 2100C) for a clinically representative range of applicator and square or rectangular insert combinations. The source-to-surface distance used was 100 cm. The field sizes ranged from 1 x 1 cm2 to 20 x 20 cm2. These phase-space files were analysed to study the intrinsic beam characteristics and used as source input for relative dose and output factor computations in homogeneous water phantoms using the code EGS4/DOSXYZ. The calculated relative central-axis depth-dose and transverse dose profiles at various depths of clinical interest agreed with the corresponding measured dose profiles to within 2% of the maximum dose. Calculated output factors for the fields studied agreed with measured output factors to about 2%. This demonstrated that for the Varian Clinac 2100C linear accelerator, electron beam dose calculations in homogeneous water phantoms can be performed accurately at the 2% level using Monte Carlo simulations.