Prezado Elias,

Fornecendo dados junto com seu código, permitindo torná-lo reproduzível (CMR), você obtém ajuda mais facilmente.

Suponho que o 17 seja em função da variável que você quer representar no grid que está exportando (predição). Estou supondo que 17 seja o índice da última coluna adicionada (cbind).

Um objeto SpatialPixelsDataFrame permite multivariáveis, mas pode ser que o formato escolhido exija a seleção de uma única variável pra compor um grid.

================================================
Éder Comunello

Em 27 de setembro de 2013 07:43, ELIAS MENDES COSTA MENDES escreveu:

Obrigado Éder, era isso mesmo, agora que rodei um comando para tirar todos os NAs deu certo.
Já aproveitando o email talvez pudesse me ajudar mais uma vez. No script que estou usando o comando que uso para salvar o mapa predito em um formato que
abra no SagaGis é: writeGDAL(grid_V[,17], "V_rna37.sdat", drivername="SAGA", mvFlag=-99999) Você sabe explicar o porque do 17?

Os comandos anteriores a esse foram.

# fazer o mapa para V%, ou aplicar o modelo para a area de estudo
grid_V <- cbind(grid1, predito) # usar "predito" ou "test1"
names(grid_V)
gridded(grid_V) <- ~X+Y # transformar "grid_V" de DataFrame para SpatialPixelsDataFrame
proj4string(grid_V) <- CRS(utm) # definir projecao

names(grid_V)
spplot(grid_V, "predito") # plotar no R apenas a coluna "predito"