x <- cbind(1, 3, 2, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 3, 1, 0, 1,0, 0, 0, 2, 0, 0,1, 0, 0, 0, 1, 0,1, 0, 0, 0, 0, 0,1, 1, 1, 0, 1, 0,0, 0, 0, 0, 1, 0,0, 0, 0, 0, 0, 1,0, 0, 1, 0, 0, 0,0, 1, 2, 0, 2, 1,1, 1, 0, 0, 0, 2,1, 1, 1, 1, 1, 0) x <- matrix(x, ncol=6) svd <- svd(x) // löschen der singulärwerte, zB svd$d[3] <- 0 svd$d[4] <- 0 svd$d[5] <- 0 svd$d[6] <- 0 // neue term2doc matrix in xx speichern: xx <- svd$u %*% diag(svd$d) %*% t(svd$v)