Due to substantial effect of classification of resource models on future mine planning, one should come with an accurate method of estimation to guarantee that the minimum error is acquired in the estimation process. The known world class Cu-Mo deposit, Sarcheshmeh Porphyry deposit (central Iran) selected as the study area. The Hypogene zone of the deposit was chosen as the space in which estimation processes should be done. The mean value of Molybdenum and Copper extracted from the top part of this zone, where sampling operations have been done on a dense grid. The correlation coefficient of 0.45 allowed going through the process of interpolation. It was shown that taking account Cu as an auxiliary variable the interpolation process, the estimation had been improved. Simple Cokriging interpolation technique is applied and it was proved that using Cu, with mean value of 0.61 percent, as secondary variable will decrease the estimation variance of Mo interpolation which has the mean value of 0.022 percent. The chief influence of this reduction appeared when the resource should be classified. Only 1% decrease was obtained when Cu used as secondary variable, but in an industrial aspect it can be of great importance as a high number of voxels in “Indicated” class changed into “Measured” one. This led to 133 Mt more Mo-ore that were added to the previous “Measured” class blocks. Also, the transition zones where the changes in class of cells have occurred are identified; these zones are mainly the places where Mo has fewer samples than Cu.