Realizou-se um estudo sobre a variabilidade espacial de dados de resistência do solo à penetração (RSP), nas camadas de 0,0-0,1 m, 0,1-0,2 m e 0,2-0,3 m de profundidade, utilizando métodos estatísticos em forma univariada, isto é, utilizando a geoestatística tradicional, construindo os mapas temáticos por krigagem ordinária para cada camada em estudo. Foi analisada a RSP na camada de 0,2-0,3 m de profundidade por meio de um modelo espacial linear (SLM), em que se consideraram as camadas de 0,0-0,1 m e 0,1-0,2 m como covariáveis, obtendo um modelo de estimação e um mapa temático por krigagem universal. Os mapas temáticos da RSP da camada de 0,2-0,3 m de profundidade, construídos por ambos os métodos, foram comparados por meio de medidas de acurácia obtidas a partir da construção da matriz de erros e da matriz de confusão. Verificou-se semelhança entre os mapas temáticos. Todos os mapas mostraram que a RSP é mais alta na região norte.
geoestatística; máxima verossimilhança; matriz de erro
