New constraints on the lepton number violating (LNV) parameters are derived from the analysis of the neutrinoless double beta (0νββ) decay in the hypothesis that this process would occur through the exchange of heavy neutrinos and/or SUSY particles. For derivation, we use new values of both phase space factors (PSFs) and nuclear matrix elements (NMEs) calculated with numerical codes developed recently, as well as the most recent experimental lifetimes. The NMEs are computed with a shell model (ShM) code for 48Ca, 76Ge, and 82Se nuclei, while at present similar ShM results are available only for the first nucleus. We compare our results with similar ones from literature, obtained with ShM, QRPA, and IBM-2 methods, and conclude that more results are needed for a relevant analysis on the validity of NMEs associated with these decay mechanisms.