A determinação da combinação de doses de nutrientes que dá uma produção máxima é um problema de interesse em ensaios de adubação. No presente trabalho, procurou-se ilustrar as análises de verossimilhança restrita (REML) e bayesiana (GS) de um delineamento em blocos incompletos parcialmente balanceados. Para essa ilustração, foram simulados dados de produção de uma cultura de citros, no primeiro ano de idade, com a combinação de doses de dois nutrientes. Para o ajuste do modelo e o cálculo de suas propriedades estatísticas, foi feita a inferência de verossimilhança restrita (REML) e a inferência bayesiana, usando a amostragem Gibbs, ambas aplicadas a um modelo linear misto. Os resultados demonstram que as duas metodologias são boas ferramentas de ajuste de modelos mistos e devem ser recomendadas para ajustes de superfícies de resposta.