Using the theoretical ambiguities inherent in the seesaw mechanism, we derive the new analytic expressions for both quadratic and linear seesaw formulae for neutrino masses at low energies, with either up-type quark masses or charged lepton masses. This is possible through full radiative corrections arising out of the renormalizations of the Yukawa couplings, the coefficients of the neutrino-mass-operator in the standard model with two-Higgs doublets, and also the QCD-QED rescaling factors below the top-quark mass scale, at one-loop level. We also investigate numerically the unification of top-b-τ Yukawa couplings at the scale M₁=0.59×10⁸ GeV for a fixed value of tan β=58.77, and then evaluate the seesaw neutrino masses which are too large in magnitude to be compatible with the presently available solar and atmospheric neutrino oscillation data. However, if we consider a higher but arbitrary value of M₁=0.59×10¹¹ GeV, the predictions from linear seesaw formulae with charged lepton masses, can accommodate simultaneousely both solar atmospheric neutrino oscillation data.