Beginning from the density functional theory, the linear response theory within the framework of time dependent density functional theory has been applied to study the electron energy loss spectra of ZnSe. The ground state is obtained using first-principles Full-Potential Linearized Augmented Plane Wave method. After settling the crystal structure, the electronic properties such as band structure and density of states are calculated. The electron energy loss spectra are found considering exchange-correlation kernels namely adiabatic local density approximation, long range contribution and the bootstrap. For the sake of completeness and comparison, results using random phase approximation are also incorporated. The inclusion of local field effects in all the kernels improve the agreement of optical and electron energy loss spectra with experimental data. Moreover, agreement in the plasmon peak positions, found using different kernels considering local field effects, with experimental data has also improved.