Realiza predição de um GeoModel

A função faz a predição para o mesmo conjunto de dados e também para novos dados de um modelo geoestatístico multivariado da classe GeoModel.

predict.geomodel(
  object,
  newdata,
  nugget = TRUE,
  krig.type = c("ordinary", "universal"),
  fast = TRUE,
  confirm = TRUE,
  confirm_id = "confirm_big_pred",
  confirm_threshold = 5000,
  intervals = FALSE,
  conf = 0.95,
  distref = "centroid",
  ...
)

# S4 method for class 'GeoModel'
predict(
  object,
  newdata,
  nugget = TRUE,
  krig.type = c("ordinary", "universal"),
  fast = TRUE,
  confirm = TRUE,
  confirm_id = "confirm_big_pred",
  confirm_threshold = 5000,
  intervals = FALSE,
  conf = 0.95,
  distref = "centroid",
  ...
)

Arguments

object

Objeto da classe GeoModel. Veja fit.geomodel.

newdata

Objeto espacial da classe sf ou sfc.

nugget

Lógico opcional informando se o efeito pepita deve ser incluído na predição. Veja detalhes.

krig.type

Caractere informando o tipo de krigagem entre "ordinary" e "universal". Padrão é 'ordinary'.

fast

Lógico informando se a predição deve ser realizada usando C++. Padrão é TRUE.

confirm

Lógico informando se um prompt deve ser levantado no terminal, caso newdata extrapole confirm_threshold.

confirm_id

Caractere utilizado apenas para captura da div em ambiente shiny.

confirm_threshold

Inteiro informando o número de observações necessário para levantar o prompt de confirmação de predição.

intervals

Lógico informando se a predição deve conter intervalos de confiança. padrão é FALSE.

conf

Decimal contendo o nível de confiança adotado. O padrão é 0.95.

distref

Caractere informando 'centroid' ou 'edge' para cálculo de distâncias. Veja detalhes.

...

Demais argumentos. Por enquanto não repassados.

Value

Lista com as predições para todas as variáveis, bem como as variâncias.

Details

Por padrão, se newdata é ausente, então presume-se que os mesmos pontos serão estimados, e portanto o efeito pepita será incluído por padrão (nugget = TRUE). Se newdata for informado, então por padrão nugget = FALSE. O efeito pepita pode forçadamente ser incluído ou omitido ao ser declarado na função.

fast=TRUE implica na utilização de um algoritmo em C++ para cálculo da predição. Por padrão o algoritmo calcula apenas a diagonal da matriz de covariâncias para estimativa da variância da predição. O uso de fast=FALSE não é recomendado, e pode acarretar em um uso intenso de memória.

confirm e confirm_threshold são utilizados para levantar uma solicitação de confirmação de predição. Caso confirm=TRUE e o número de observações (considerando todas as variáveis) ultrapassarem confirm_threshold, então a solicitação será levantada no terminal (ou na aplicação shiny, se em ambiente shiny).

Se intervals=TRUE, então os intervalos de confiança da predição são retornados para o nível conf de confiança informado.

O argumento distref recebe "centroid" ou "edge", e é passado para get.dists(). Tal argumento regula a posição de onde as distâncias entre objetos vetoriais é calculada. Se distref="centroid", então as distâncias são calculadas a partir do centroide das geometrias. Se distref="edge", as distâncias são calculadas a partir da borda das geometrias. Note que naturalmente objetos do tipo ponto terão suas distâncias calculadas sempre a partir do centroide.

Author

Gabriel Agostini Orso