[R-br] filliben

[Hélio Gallo Rocha]

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 em 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] <http://user/SendEmail.jtp?type=node&node=4661205&i=0>
> 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<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
> .
> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>



-- 
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/20131207/1cedb6ed/attachment.html>