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