Senhores, 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. O grande problema é que nao sei isolar Beta(3x1) na seguinte equação: H'XBeta = H'Y = lambda'Beta, onde lambda'(1x3) Eu tenho H'(1x3), X(3x3) e Y(3x1), mas nao estou visualizando uma maneira de isolar beta, pois H'X é um vetor. X <- matrix(c(3,0,-3,0,6,3,-3,3,6),ncol=3,nrow=3) dvs <- svd(X) lambda <- (dvs$d)^2 Y <- c(5,2,3) D <- matrix(c(lambdas[1,1],0,0,0,lambdas[2,1],0,0,0,lambdas[3,1]), ncol=3, nrow=3) C <- dvs$u Vlinha <- t(dvs$v) 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) Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)