En la actualidad el modelamiento geológico carece de una cuantificación de la incertidumbre en la aparición de las unidades geológicas (UG) en el espacio. Junto a ello, la estimación de leyes minerales es efectuada independientemente entre las unidades geológicas modeladas. De este modo, el control ejercido por las unidades geológicas en la estimación de leyes se reduce a la pertenencia de las muestras a cierta unidad geológica y al volumen de la unidad modelada, sin explotar cabalmente la interrelación existente entre UG y leyes.Este trabajo presenta un método para establecer los recursos minerales e incertidumbre asociada, vinculando las leyes minerales y unidades geológicas por medio de la simulación geoestadística conjunta (cosimulación) de variables categóricas y continuas. El método propuesto utiliza el modelo multigaussiano para variables continuas (leyes) y gaussiano truncado para variables categóricas (UG), los cuales son combinados para caracterizar la interdependencia espacial de las unidades geológicas y leyes. El método además permite la incorporación de conocimiento geológico adicional (no medido directamente) en la estimación de leyes.El método propuesto y los enfoques tradicionales existentes son aplicados en un caso de estudio en el sector Sur Sur del yacimiento Río Blanco Los Bronces, utilizando las leyes de cobre asociadas a mineralización primaria y la litología registradas en datos de sondajes. Los modelos realizados son comparados y contrastados con datos de producción.Los resultados indican que en el enfoque propuesto, se expresa claramente que la ocurrencia de cierta unidad litológica influye en las leyes simuladas y, de modo inverso, el valor de las leyes influye en la posible aparición de cierta unidad. Por otro lado, la incorporación de conocimiento geológico adicional ejerce control sobre las leyes simuladas. De este modo, el modelo permite transferir la certeza en la aparición de cierta unidad litológica a las leyes asociadas a ésta.Debido a la hipótesis de independencia asumida entre las leyes de unidades contiguas, los modelos tradicionales generan simulaciones con cambios abruptos al pasar de una unidad litológica a otra, mientras que el modelo propuesto presenta transiciones graduales, más representativas de lo observado en las muestras. Los resultados indican que el modelo propuesto presenta un mejor desempeño que los enfoques tradicionales en términos de predictibilidad en leyes, el cual se incrementa al incorporar el conocimiento geológico adicional. Finalmente el modelo planteado se presenta como una alternativa interesante a los modelos tradicionales, en particular en depósitos con mineralización diseminada y/o donde tengan lugar transiciones graduales entre las leyes de unidades en contacto.