bestfitmod is a model fitting function internally called in the code bdynsysdysymod. It calls the functions polyfitreg. It creates the testing and training datasubsets for model fitting, passes them on to polyfitreg and receives back the Sum Square of Errors which are then stored in matrix for each indicator. It returns internally the Sum Square of Error matrices to bdynsysdysymod. It requires the package caTools.