[R-br] filliben

[Hélio Gallo Rocha]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131206/3f328445/attachment.html>