[R-br] erro na validação cruzada com pacote geoR (Éder Comunello)

[Paulo Justiniano]

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
> R-br em listas.c3sl.ufpr.br
> 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.
>