comparação entre: ginv(x) com solve( t(x)%*% x ) %*% t(x)

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 --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

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@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@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.
participantes (2)
-
Cesar Rabak
-
Cleber N.Borges