To solve the model Y = XB, I use the function:
object <- pls(X, Y, ncomp = 2)
being X the independent matrix with 50 continuous variables, and Y a unique dependent variables.
I would like to include fixed effects (categorical) in the PLS model, apart from the X matrix ( Y = diet + XB), and I don't know how.
I would be very grateful if you could help me Thanks in advance,