Caro Elias, obrigado pela demonstração.

Estou com uma dificuldade, se puder dar uma ajuda, agradeço muito, desde já.

No CT do Carvalho e Vieira eles usamos resíduos ortonomais,  da seguinte forma

Ek= dados$data-kc$pred/sqrt(kc$kr)   pelo menos foi o que entendi, me corrija se entendi errado.

No link que citei  faz-se a regressão, que ao meu modo de entender seria isso:

reg=lm(dados$data~Ek)    , por favor, me corrija se entendi errado.

para depois calcular a correlação de filliben e comparar com o tabelado e assim verificar se o modelo atende a hipótese dos resíduos ortonomiais normalmente distribuídos, com:

filliben(rstandard(reg))

E por ultimo, uma dúvida de comando no R, se puder me dar uma luz, agradeço mais uma vez!

Vamos supor que meus dados possuam 50 pontos com coordenadas.
Minha Krigagem gerou 1000 pontos.

Evitando usar o predict e krige.var da validação, pergunto:

Teria uma função que separasse 50 pontos da krigagem com as mesmas coordenadas dos dados originais?

Não necessariamente entre os 1000 pontos da krigarem terão as mesmas coordenadas das coordenadas dos meus dadosoriginais,
Seria válido utilizar a coord da krig mais próxima possível da coord do dado original, se sim, qual o limite de distância?

Abraço

Hélio






2013/12/6 Elias T Krainski [via R-br] <ml-node+s2285057n4661205h23@n4.nabble.com>
apenas para ilustrar o poder desse tipo de teste:

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)
}

e = replicate(1000, filliben(rexp(100, 1)))
u = replicate(1000, filliben(runif(100)))
n = replicate(1000, filliben(rnorm(100)))

hist(n, xlim=range(e,u,n))
hist(u, add=TRUE, border='blue')
hist(e, add=TRUE, border='red')

_______________________________________________
R-br mailing list
[hidden email]
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.



If you reply to this email, your message will be added to the discussion below:
http://r-br.2285057.n4.nabble.com/R-br-filliben-tp4661203p4661205.html
To unsubscribe from R-br, click here.
NAML



--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho