Consider N independent stochastic processes (Xi (t), t. [ 0, T]), i = 1,..., N, defined by a stochastic differential equation with random effects where the drift term depends linearly on a random vector Phi i and the diffusion coefficient depends on another linear random effect Psi i. For these effects, we consider a joint parametric distribution. We propose and study two approximate likelihoods for estimating the parameters of this joint distribution based on discrete observations of the processes on a fixed time interval. Consistent and v N-asymptotically Gaussian estimators are obtained when both the number of individuals and the number of observations per individual tend to infinity. The estimation methods are investigated on simulated data and show good performances.