Performs SVD on a data matrix using the base svd() function in R, with a workaround to avoid LAPACK errors. If an SVD of the data matrix gives an error, an SVD of its transpose will be performed. Used internally when computing the JIVE decomposition. Credit to Art Owen: https://stat.ethz.ch/pipermail/r-help/2007-October/143508.html.
svdwrapper(x, nu, nv, verbose=F )a numeric matrix whos SVD decomposition is to be computed.
the number of left singular vectors to be computed.
the number of right singular vectors to be computed.
logical. Print error message if needed.
An svd object, as returned by svd(x,nu=nu,nv=nv).
# NOT RUN {
x<-matrix(rnorm(100),nrow=10,ncol=10)
SVD = svdwrapper(x,nu=1,nv=1)
# }
Run the code above in your browser using DataLab