Uno de los metodos utilizados para evaluar estabilidad fenotipica es el propuesto por Shukla, el cual calcula la varianza de los genotipos dentro de la interaccion genotipo por ambiente, para lo cual se hace uso de la estimacion de los componentes de varianza dentro del analisis de varianza combinado. El acceso al calculo de metodologias como REML e ML permitieron trabajar con datos que presentan algun grado de desbalance, sin embargo no solucionan de una manera adecuada el problema de la estimacion de componentes de varianza negativos, los cuales son asumidos como cero y redistribuidos en los demas componentes positivos. El uso de la metodologia bayesiana en la estimacion de componentes de varianza resuelve satisfactoriamente este problema sin afectar los demas componentes. En este trabajo, se utilizaron datos de produccion comercial de papa de 10 pruebas regionales realizadas en la region andina colombiana y se utilizo la metodologia bayesiana en la solucion del modelo mixto para la estimacion de la varianza de Shukla con base en una distribucion a priori no informativa de Jeffreys. Fueron obtenidas muestras de la distribucion a posteriori conjunta mediante el algoritmo Independence Chain, con un tamano de muestreo de 1,16x105 y un burn-in de 500. Los resultados muestran que en la estimacion REML de componentes de varianza tres genotipos presentan componentes de varianza estimados como cero. Las estimativas bayesianas son 89,35; 377,18 y 101,12; y los respectivos intervalos de credibilidad al 95% son: (2,13 – 371,70), (35,26 – 1363,67) y (2,33 – 434,53). Finalmente con estas estimativas no se afecto la estimacion de los demas componentes de varianza.