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
_______________________________________________ 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.
