<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Se em <br></span></div><div><span> A b = x</span></div><div><span>A tem dimensão 3x3, então vc quer encontrar 6 valore. Para isso, vc tem 3 valores em b e 3 valores em x. Isso é possível? No usual, quando vc quer saber os 3 valores de b e conhece A e x, vc tem bem mais valores conhecidos que desconhecidos. </span></div><div><br></div><div>### Exemplo <br></div><div>### considerando o contexto de LM<br>n <- 5<br>X <- cbind(1, runif(5))<br>b <- c(10, 5)<br>y <- rnorm(n, X%*%b, 0.5)<br></div><div><br></div><div>### estima b<br>(XX <- crossprod(X))<br>(Xy <- crossprod(X, y))<br>(bhat <- solve(XX, Xy))<br><br>### com bhat e Xy podemos obter XX?<br>### considere que<br>### XX %*% bhat = Xy<br>### para compatibilidade de operações seja<br>### XX %*% (bhat %*% bhat') =
Xy %*% bhat'<br><br>require(Matrix)<br>bb <- forceSymmetric(tcrossprod(bhat))<br>solve(bb, tcrossprod(Xy, bhat)) ## nao é igual a XX</div><div><br></div><div>pode ser que este 'sistema' possua mais de uma solução...<br><br></div><div>Att.<br><span></span></div><div>Elias T. Krainski<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; margin-top: 5px; padding-left: 5px;"> <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"> <font face="Arial" size="2"> <hr size="1"> <b><span style="font-weight: bold;">De:</span></b> FHRB Toledo <fernandohtoledo@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> R-Br <r-br@listas.c3sl.ufpr.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Domingo, 11 de Dezembro de 2011 12:07<br> <b><span style="font-weight:
bold;">Assunto:</span></b> [R-br] MQM<br> </font> <br><meta http-equiv="x-dns-prefetch-control" content="off"><div id="yiv651420565">Senhores,<br><br>Tenho uma dúvida que não necessáriamente se relaciona com o programa, mas que em muito tem a ver com teoria estatística e algebra de matrizes, espero que aqui ache pelo menos uma sugestão para a resposta ou uma recomendação de literatura.<br>
<br>O problema é dado o sistema de equações na forma matricial: Xb = y, pelo MQM, dada a matriz X e o vetor y, obtenho o vetor de soluções por X^-1y, respeitando todas as pressuposições.<br><br>Se caso conheço ao invés de X e y, conheço os vetores b e y, existe alguma forma de conhecer a matriz X, podendo apenas assegurar que ela é uma matriz quadrada simétrica e positiva definida?<br>
<br>att,<br>FH<br>
</div><meta http-equiv="x-dns-prefetch-control" content="on"><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </blockquote></div> </div></body></html>