
Complementando informações, Disponibilizei mais dois artigos e o artigo do próprio Filliben, onde tem a tabela de percentuais do coeficiente de correlação. https://www.dropbox.com/sh/jbd2pmvzm41hs6b/Hk2tm7fr5v Hélio Em 6 de dezembro de 2013 08:36, Hélio Gallo Rocha <heliogallorocha@gmail.com
escreveu:
Caros listeiros, bom dia.
Estou com uma
Li um comunicado técnico da EMBRAPA utilizando o Teste de Filliben, onde analisa os resíduos ortonormais para validação modelos geoestatísticos. neste link: https://www.dropbox.com/s/e9ww5h5w9dz8wed/CT60_FILLIBEN.pdf
Fiz uma busca na lista e não achei nenhuma referencia.
Achei neste link: http://genepi.qimr.edu.au/staff/davidD/R/filliben.R
esta função:
filliben <- function(x) { y<-sort(x) n<-length(y) p<-(1:n-0.3175)/(n+0.365) p[1]<-1-0.5^(1/n) p[n]<-0.5^(1/n) p<-qnorm(p) cor(y,p) }
Para obter o resultado do teste:
filliben(rstandard(reglin_esc))
onde, segundo o link reglin_esc seria produto de uma regressão, pensei entre os dados originais e os estimados pela krigagem, até ai tudo bem.
No Comunicado técnico ele calcula o resíduo ortonomal.
Alguém poderia dar uma ajuda?
Agradeço antecipadamente,
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho