[R-br] Problema com matrizes e vetores

[Thiago Touya]

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)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110928/689c2cdc/attachment.html>