Caro Wenceslau
Aproveitando a discussão para aprender sobre a validação de modelos em
geoestatisticos. Dado que estest modelos e tecnicas são novas para mim.
A krigagem é um estimador exato se o numero de pontos estimados for o mesmo dos
avaliados, não ?
CVes posta enviado agora pouco sobre m=variavcao de nmicro escala e erro
de medida. O interpoklacor só é exato se o nugget é considerado ser o
primeiro destes
Se eu estimar 10.000 valores para um universo de 1000 valores medidos.
Os valores estimados para os 1000 pontos avaliados serão os mesmos
valores ?
ver acima
Em adicao note que no ajuste de max verossimilhaça
o argumento " components" gasrante o retorno daq decomposição das
observacoes segundo o modelo ajustado
args(likfit)
function (geodata, coords = geodata$coords, data = geodata$data,
trend = "cte", ini.cov.pars, fix.nugget = FALSE, nugget = 0,
fix.kappa = TRUE, kappa = 0.5, fix.lambda = TRUE, lambda = 1,
fix.psiA = TRUE, psiA = 0, fix.psiR = TRUE, psiR = 1, cov.model,
realisations, lik.method = "ML", components = TRUE, nospatial = TRUE,
limits = pars.limits(), print.pars = FALSE, messages, ...)
Saudações geoestatisticas
Wenceslau Teixeira
Embrapa Solos
_________________________________________________
Prezado Helio, boa tarde!
Pensando melhor sobre a questão do resíduo, o que coloquei no email
anterior não faz muito sentido, já que a krigagem é um estimador exato e a
operação retornaria praticamente o próprio valor inicial (e o resíduo
tenderia ser zero/nulo). O que faria sentido seria separar o efeito da
tendência ('drift') e descontar esse valor dos dados originais, tendo aí
sim o resíduo a ser modelado (sem o efeito da tendência). Para o produto
final teria que somar o drift novamente para retornar os valores.
Acredito que não seja muito difícil de fazer no R, mas só tentando fazer
pra ver.
Fui dar uma conferida no seu post antigo sobre o erro na validação cruzada:
<
http://r-br.2285057.n4.nabble.com/R-br-erro-na-validacao-cruzada-td4659526.html#a4659536
Mas o link dos dados não está mais disponível...
<https://www.dropbox.com/s/tfml7zjpc7hmbc3/dados250.txt>
Se importaria em dispor esses dados novamente?
Atte.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
