A finite-size scaling analysis of wave functions near the metal-insulator transition ͑MIT͒ point has been developed, and applied to the MIT in a two-dimensional disordered electron system in the presence of spinorbit interaction. The present method has the following advantages: ͑i͒ Quantities characterizing the critical behavior, such as the critical disorder W c or the localization exponent , are multiply calculated from independent scaling analyses of spatial parts with different intensities in wave functions. ͑ii͒ These quantities and the multifractality of the critical wave function are determined simultaneously. ͑iii͒ It is not necessary to treat many samples with different sizes. ͑iv͒ Much computing time is saved, and the scaling analysis can be done up to very large sizes. Using this method, we obtained W c ϭ5.86Ϯ0.04 and ϭ2.41Ϯ0.24 for a model of a two-dimensional symplectic system. ͓S0163-1829͑98͒08039-4͔