In order to better understand the differential code biases (DCBs) of global navigation satellite system, the IGGDCB method is extended to estimate the intra- and inter-frequency biases of the global positioning system (GPS), GLONASS, BeiDou navigation satellite system (BDS), and Galileo based on observations collected by the multi-GNSS experiment (MGEX) of the international GNSS service (IGS). In the approach of IGGDCB, the local ionospheric total electronic content is modeled with generalized triangular series (GTS) function rather than using a global ionosphere model or a priori ionospheric information. The DCB estimated by the IGGDCB method is compared with the DCB products from the Center for Orbit Determination in Europe (CODE) and German Aerospace Center (DLR), as well as the broadcast timing group delay (TGD) parameters over a 2-year span (2013 and 2014). The results indicate that GPS and GLONASS intra-frequency biases obtained in this work show the same precision levels as those estimated by DLR (about 0.1 and 0.2–0.4 ns for the two constellations, respectively, with respect to the products of CODE). The precision levels of IGGDCB-based inter-frequency biases estimated over the 24-month period are about 0.29 ns for GPS, 0.56 ns for GLONASS, 0.36 ns for BDS, and 0.24 ns for Galileo, respectively. Here, the accuracies of GPS and GLONASS biases are assessed relative to the products of CODE, while those of BDS and Galileo are compared with the estimates of DLR. In addition, the monthly stability indices of IGGDCB-based DCBs are 0.11 (GPS), 0.18 (GLONASS), 0.17 (BDS), and 0.14 (Galileo) ns for the individual constellation.