2
Por que não consigo obter um SVD válido de X por decomposição de autovalor de XX 'e X'X?
Estou tentando fazer SVD manualmente: m<-matrix(c(1,0,1,2,1,1,1,0,0),byrow=TRUE,nrow=3) U=eigen(m%*%t(m))$vector V=eigen(t(m)%*%m)$vector D=sqrt(diag(eigen(m%*%t(m))$values)) U1=svd(m)$u V1=svd(m)$v D1=diag(svd(m)$d) U1%*%D1%*%t(V1) U%*%D%*%t(V) Mas a última linha não retorna m. Por quê? Parece ter algo a ver com os sinais desses autovetores ... Ou entendi mal o procedimento?
9
r
svd
eigenvalues