Senhores, <br><br>preciso estimar um vetor de parâmetros (Beta) usando as saídas do svd() para comparar com a estimação pelo método dos mínimos quadrados. <br>O grande problema é que nao sei isolar Beta(3x1) na seguinte equação:<br>
<br>H'XBeta = H'Y = lambda'Beta,   onde lambda'(1x3)<br><br>Eu tenho H'(1x3), X(3x3) e Y(3x1), mas nao estou visualizando uma maneira de isolar beta, pois H'X é um vetor.<br><br><br clear="all">X <- matrix(c(3,0,-3,0,6,3,-3,3,6),ncol=3,nrow=3)<br>
dvs <- svd(X)<br>lambda <- (dvs$d)^2<br>Y <- c(5,2,3)<br>D <- matrix(c(lambdas[1,1],0,0,0,lambdas[2,1],0,0,0,lambdas[3,1]), ncol=3, nrow=3) <br>C <- dvs$u<br>Vlinha <- t(dvs$v)<br><br>sabe-se que pela Decomposição em valores singulares X=CDV' (neste caso em que X é simétrica e nao temos nenhum autovalor igual a zero)<br>
<br><br>Obrigado.<br><br><br><br><br>-- <br>Att,<br>Thiago Morais de Carvalho (61 - 8169-5094)<br><br>