Multichannel near-infrared spectroscopy (NIRS) is a functional neuroimaging modality that enables easy-to-use and noninvasive measurement of changes in blood oxygenation levels. We developed a clinically-applicable method for estimating resting state functional connectivity (RSFC) with NIRS using a partial correlation analysis to reduce the influence of extraneural components. Using a multi-distance probe arrangement NIRS, we measured resting state brain activity for 8 min in 17 healthy participants. Independent component analysis was used to extract shallow and deep signals from the original NIRS data. Pearson's correlation calculated from original signals was significantly higher than that calculated from deep signals, while partial correlation calculated from original signals was comparable to that calculated from deep (cerebral-tissue) signals alone. To further test the validity of our method, we also measured 8 min of resting state brain activity using a whole-head NIRS arrangement consisting of 17 cortical regions in 80 healthy participants. Significant RSFC between neighboring, interhemispheric homologous, and some distant ipsilateral brain region pairs was revealed. Additionally, females exhibited higher RSFC between interhemispheric occipital region-pairs, in addition to higher connectivity between some ipsilateral pairs in the left hemisphere, when compared to males. The combined results of the two component experiments indicate that partial correlation analysis is effective in reducing the influence of extracerebral signals, and that NIRS is able to detect well-described resting state networks and sex-related differences in RSFC.