We have derived the luminosity function for subdwarfs on the basis of the proper motion data in LHS Catalogue, utilizing the reduced proper motion diagram for the selection of sub dwarfs and the hybrid method combining the mean absolute magnitude method and V/$V_m$ method to estimate the distance and density of subdwarfs. The luminosity function found here is almost flat, showing a very slow increase up to $M_V;=;9$ or $M_B;=;10$, and the overall halo density is larger than those derived by Schmidt (1975), Chiu (1980), Reid (1984), Lee (1985), and Dawson (1986), but smaller than that by Eggen (1983). Comparison with 1/100 of disk stellar luminosity function implies that no conclusive dip in the halo luminosity function is found.