The neutrino masses are generated in grand unified theory (GUT) constrained supersymmetric model with R-parity violation. The neutrinos acquire masses via tree-level neutrino-neutralino mixing as well as via one-loop radiative corrections. The theoretical mass matrix is compared with the phenomenological one, which is reconstructed by using neutrino oscillation and neutrinoless double beta decay data. This procedure allows to obtain significantly stronger constraints on R-parity breaking parameters than those existing in the literature. The implication of normal and inverted neutrino mass hierarchy on the sneutrino expectation values, lepton-Higgs bilinear and trilinear R-parity breaking couplings is also discussed.