n <- 6 ; p <- 4
X <- matrix(rnorm(n * p), ncol = p)
y <- rnorm(n)
nlv <- 3
fm <- pcasvd(X, nlv = nlv)
xfit(fm, X)
xfit(fm, X, nlv = 1)
xfit(fm, X, nlv = 0)
X - xfit(fm, X)
xresid(fm, X)
X - xfit(fm, X, nlv = 1)
xresid(fm, X, nlv = 1)
Run the code above in your browser using DataLab