[R-br] comparação entre: ginv(x) com solve( t(x)%*% x ) %*% t(x)

[Cesar Rabak]

A resposta a essa pergunta se centra na definição de "igualdade":

1) SE você deseja a igualdade com acurácia numérica exata, a resposta é
"nunca";

2) SE você aceita uma "faixa" de "precisão" dos resultados, a resposta é
"na maioria dos casos".

Para mais detalhes procure (tanto nos arquivos desta lista, que estarão em
português) como na Net sobre a precisão numérica em computação e a norma
IEEE 754.

HTH
--
Cesar Rabak


2018-03-16 15:39 GMT-03:00 Cleber N.Borges via R-br <
r-br em listas.c3sl.ufpr.br>:

> olá a todos
>
> Ainda que a dúvida seja mais conceitual do que do R propriamente.... peço
> ajuda aqui
>
> Minha pergunta é:
>
> Até que ponto eu posso usar essa igualdade?? Ou mesmo se posso considerar
> isso igual...
>
> #############################################
>
> library( MASS ) # para usar a função ginv
>
> x <- matrix( runif(100*100), 100, 100 )
>
> all.equal( ginv(x),  solve( t(x)%*% x ) %*% t(x) ) #  TRUE
> #############################################
>
> O motivo de minha dúvida e vontade de usar isso nos modelos lineares com
> variáveis linearmente dependentes
>
> do tipo de misturas químicas, onde a soma das proporções logicamente
> sempre tem que ser 100%.
>
> Por hora, exponho somente essa dúvida (entre tantas :-) )
>
> Agradeço antecipadamente por qualquer indicação nesse tema de 'ginv'
>
> Muito Obrigado
>
> Cleber Borges
>
>
>
>
>
>
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient> Livre
> de vírus. www.avast.com
> <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=emailclient>.
> <#m_6061262705293632801_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180317/b86176e9/attachment.html>