An analysis of the electron localization properties in doped graphene is performed by doing a numerical multifractal analysis. By obtaining the singularity spectrum of a tight-binding model, it is found that the electron wave functions present a multifractal behavior. Such multifractality is preserved even for second neighbor interaction, which needs to be taken into account if a comparison is desired with experimental results. States close to the Dirac point have a wider multifractal character than those far from this point as the impurity concentration is increased. The analysis of the results allows to conclude that in the split-band limit, where impurities act as vacancies, the system can be well described by a chiral orthogonal symmetry class, with a singularity spectrum transition approaching freezing as disorder increases. This also suggests that in doped graphene, localization is in contrast with the conventional picture of Anderson localization in two dimensions, showing also that the common belief of the absence of quantum percolation in two dimensional systems needs to be revised.Ever since its discovery 1,2 , graphene has been considered as an ideal candidate to replace Silicon in electronics 3 , since this first truly two-dimensional crystal has the highest electrical and thermal conductivity known 4 . However, graphene per se is not a semiconductor. Several proposals have been made to solve this issue 5 . Experimentally, it has been found that doped graphene presents a metal to insulator transition 6 when doped with H, producing a kind of narrow band gap semiconductor. The increase in localization around the Dirac point was roughly predicted from an electron wavefunction frustration analysis in the graphene's underlying triangular lattice 7-9 . Such theoretical results were made under the supposition that Hydrogen bonds to the 2p z Carbon orbital, and thus impurities act as vacancies 10,11 . This case corresponds to the split-band limit. This approach has been useful to predict localization and the pseudogap size, i.e., the region in which the inverse participation ratio increases by one order of magnitude 7 , in very good agreement with experiments 6 , although vacancies and impurities are indeed different 12,13 . However, there is a theoretical nuance to the idea of having a metal-insulator transition in two dimensions (2D). According to the well known Abrahams's et. al. scaling analysis, in 1D and 2D all states are localized for any amount of disorder, excluding the possibility of a mobility edge 14 . The experimental and numerical analysis shows that not all states are localized, and there exists a kind of mobility edge associated with a pseudogap around the Fermi energy 9,15,16 . This means that there is a problem that has not been solved. There are two possibilities, either electron-electron interaction produces delocalization 17 , or somehow the Abrahams's et. al. analysis does not completely applies to this case. In graphene, electron-electron interaction is very weak 18 , so in principle, t...