Hi again @hchen8, so far so good! 
Then it could be due to:
- Too many missing or zero values.
nearZeroVar()can be used to handle this, either externally or internally in the pls function. - Too many components (empty residual matrices)
- Variables that are highly correlated/nearly the same (multi collinearity)
Let me know if any of these solved the problem.
- Christopher