Greetings R users,
My interest in the Q2cum score comes my endeavor to replicate SIMCAP PLS-DA
analysis in R. I use the exact same dataset. After doing the analysis in
R, I can get the exact same R2Ycum. However, the Q2cum is significantly
off. Adding the Q2cum of the 1st and 2nd component comes close but that
seems unconvincing and I don't understand why the 3rd component Q2cum would
be so far off. Below is the code and output:
pls.reg.new=plsreg2(newserumvars.sc, newtimematrix, nc=3)
pls.reg.new$expvar #this provides the exact R2Y I get in SIMCAP
pls.reg.new$Q2cum
The Q2cum values should be close to
component 1 = .2149
component 2 = .4578
*component 3 = .6627* #value of concern
All of the R output values are close for the respective class except for
the final Q2cum value:
Q2cum.B Q2cum.FR8 Q2cum.S45 Q2cum
t1 -0.0546 0.3992 0.2973 0.2140
t2 0.2863 0.3824 0.7039 0.4573
t3 0.5839 0.7430 0.6980 *0.3735 *#*not close*
Done.
I have read that the plspm package used LOO cross validation. SIMCAP uses
the same method, however it leaves '7' out instead. Perhaps this might
be
an issue that results in this difference?
Any insight or assistance would be most appreciated. I certainly
appreciate any time you take to answer my question.
Regards,
Charles Determan
[[alternative HTML version deleted]]
