Olá pessoal da lista!

Estou com dificuldade de implementar na função krige.bayes() o atributo trend. Primeiramente, por uma abordagem frequentista escolhi um modelo de trend, quando tento usar o mesmo na função krige.bayes() ocorre alguns erros com questões matriciais. Alguem tem alguma dica para corrigir isso?

Outra questão cé omo obtendo o DIC para os modelos que testarei na função krige.bayes(). Por exemplo, modelo com trend e componentes espaciais, modelo só com trend e modelo só com os componentes espacias.

A figura abaixo mostra um exemplo de um artigo que fez essas comparações utilizando binom.krige.bayes(). Fonte: Mapping, Bayesian Geostatistical Analysis and Spatial Prediction of Lymphatic Filariasis Prevalence in Africa

Wagner Wolff, PhD
"Luiz de Queiroz" College of Agriculture,
University of São Paulo
Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil
Phone: +55 19 982385582
http://orcid.org/0000-0003-3426-308X
https://github.com/wwolff7
http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1