Poderia enviar o arquivo correspondente para poder rodar?
 
 
Prezados colegas, bom dia!
 

 
Então, é o seguinte, no site http://berkeleyearth.org/ eu tenho dados de temperatura média de toda a terra. Eu preciso extrair de acordo com os meus pontos geográficos, a temperatura média anual para cada ponto. Só que eu tenho os específicos que eu preciso. Comecei o script de forma que eu pudesse extrair para um ponto específico e para todos os anos, para depois arrumar o script para fazer para todos os meus pontos e anos de interesse.
 
Segue o script utilizado!
 
library(raster)
 
b <- brick("D:/Documentos/Desktop/analises_cap2/dados_climaticos/precipt.nc", varname <- "precip")
print(b)
b
 
idx <- getZ(b)
coords <- matrix(c(-55, -10), ncol = 2) 
vals <- extract(b, coords, df = T)
df <- data.frame(idx, t(vals)[-1,])
 
rownames(df) <- NULL
names(df) <- c('date','value')
head(df)
 
Seguindo o script no local onde deveria apresentar as variáveis extraídas , aparece NA. Alguém já fez este tipo de extração? Peço ajuda. 
 
Desde já, muito obrigada.
 
Ana Paula